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ABSTRACT 


It is well known that if the ratio of specific heats, I’, in a star is a constant and smaller than #, the 
star is dynamically unstable. In this paper the case in which I is variable in the star is discussed in 
some detail, especially for the homogeneous and standard models. It is shown that if ! becomes smaller 
than ¢ in a part of the star, that part has to be very extensive to render the star unstable. For example, 
if ! = }in the central part of a star and is as small as 1 in the outer part, this has to extend to a depth 
at which the temperature is of the order of half the central temperature. Peculiar forms of I(r) could in- 
crease the instability considerably, but it is doubtful whether any of these could be of physical im- 
portance. 


1. General considerations.—It has been known for a long time! that a star with a con- 
stant ratio of specific heats, I’, smaller than 4 is dynamically unstable. If T is variable, 
this should apply to a certain mean value of I’. This becomes clearer if we consider a 
formula such as the one given by the author in a recent paper.® For a small radial 
adiabatic deformation, such that 


6 


o” is given by 


0 0 po 
(1) 


* Fellow of the Belgian-American Educational Foundation, at the Yerkes Observatory. 


1 Ritter, Anwendungen des mechanisches Warmetheorie auf kosmologische Probleme, 1879. The result 
derived there concerns a homogeneous star. For a more general result see, e.g., S. Rosseland, Pub. Oslo 
U. Obs., No. 1, p. 20; or S. Chandrasekhar, An Introduction to the Study of Stellar Structure, 1939, p. 52, 
and references given there. 


2 Cf. L. Biermann and T. G. Cowling, Zs. f. Ap., 19, 1, 1939 (first part of the paper, where a formula 
of type [1], except for the term in dI’/dro, is used). 


3P. Ledoux, Ap. J., 102, 56, 1945. 
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where the suffix zero refers to equilibrium values. Furthermore, 


Gm (r) dm R 


To 


= — 


represent, respectively, the gravitational potential energy and the moment of inertia 
with respect to the origin, and é is the solution of the differential equation 


@t, f4—p,1 dP] dé , 


satisfying the boundary conditions 


Ero =0 at ro =0 (3) 
and 


sp = 5*] = — (4) 


In equation (2), 1 = Gm(r)po/roPo, and T is the general adiabatic exponent defined by 
dQ = dU — Pdp/p* = 0, where the variation of the internal energy dU is expressed in 
terms of P and p. In a star, I is, in general, a function of the ratio of the pressure of 
radiation to the total pressure (1 — 8), of the degree of ionization, and of the number of 
degrees of freedom of the particles.‘ 

Equation (2) reduces to Eddington’s equation for small radial pulsations if I’ is con- 
stant. Formula (1) can also be written as 


R R R dé 


if we integrate the last term of the numerator by parts. 

If o?, as given by equation (1) or equation (5), is positive, a small radial oscillation will 
result, which will be damped out or not according to the condition of vibrational stabil- 
ity.® But if o? is negative, the star will be violently unstable, as the small disturbance 
considered will grow exponentially with the time. 

Cases have been considered in which, because of a very small abundance of hydrogen 
and helium, I becomes smaller than 4 in an appreciable part of the star.® The stability 
of these models has been investigated? by a general method, using the minimal property 
of the total available energy of the star in case of stability. However, the rigorous ap- 
plication of such a condition is rather difficult, and it seems worth while to re-examine 
the problem from the point of view of the sign of 0”. 

If T is a constant and if £ is the amplitude of the fundamental mode, then it readily 
follows from equation (1) that the star is unstable for this mode if I’ is smaller than 4. 

The case of the harmonics, however, is not so simple. If I’? denotes the critical value 
of I for the harmonic of order , then, on the basis of the general theory of equations 
of type (2), we can expect that I? will be smaller, the higher the order of the harmonic 
considered. We shall see that this is true in the examples which have so far been dis- 
cussed in some detail. Thus, in general, we may expect that if is a constant for the whole 


4R. H. Fowler and E. A. Guggenheim, M.N., 85, 961, 1925. 


5 Cf, e.g., T. G. Cowling, M.N., 94, 768, 1934; 96, 42, 1935; 98, 528, 1938; and P. Ledoux, Ap. J., 94, 
537, 1941. 


6 L. Biermann, Zs. f. Ap., 18, table on 356, 1939; also n. 2. 
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star, it is through the fundamental mode that the stability of the star will first be en- 
dangered when we go along a series of models with decreasing I’. 

Considering, next, the general case of variable I’, we shall first re-write equation (2) 
in the form 


where 
p=IP,r4 and g=[ Gr—4) Gm —3np, 


Now let I’, be the minimum value of I along r. Denoting by a suffix m the values cor- 
responding toT’,,, we have pm < pand, if dI'/dro is everywhere negative or zero, gm < q. 
Hence o2 > (02) m. 

If dI'/drg is positive in some regions of the star, we cannot say anything rigorous about 
the relative values of g and gm. But, in general, the second term in g is small compared 
to the first one, and we can expect that in most cases the preceding inequalities will 
hold. Thus if I’, is greater than I, the star will be stable toward the harmonics of 
order n or greater, and the fundamental mode remains the most dangerous one for the 
stability of the star. Since P’8 = 4, the star will be stable if I is greater than 4 all over the 
Star. 

However, there might be exceptions; for, writing equation (1) in the form 


R R 
f tdlo= — f GP—4)P ol dro, (1” 


let us consider again the fundamental mode. The pressure P is essentially positive and 
decreases from the center to the surface, where it can be taken as being zero. If (3f — 4) 
is everywhere positive (I' > 4), the star will certainly be stable if (3 — 4)Po decreases 
everywhere with increasing r. 

But if (3% — 4) Po, for instance, remains constant, then o? = 0, and the star is on the 
verge of instability, although I is greater than $ everywhere. Of course, this particular 
case implies an infinite I at the surface, and physically it has no meaning. 

In fact, if (3 — 4) is everywhere positive, we can admit for physical reasons only 
the range of I, 0 < 32 — 4 < 1; and (31 — 4)P» will be an oscillating function in r, 
each oscillation corresponding to the ionization of a new shell of electrons of some fairly 
abundant element. If (3 — 4)Po starts at the center with its maximum value (this 
is likely as long as the central temperature remains high), these oscillations will be super- 
posed on a generally decreasing curve. And if there should be m regions of increasing 
(31 — 4)Po, there will be (7 + 1) regions of decreasing (3 — 4)Po—and one of them 
just near the surface. Starting from the surface, we can associate them by pairs (de- 
creasing, increasing), and the decreasing one will have a greater weight than the increas- 
ing one, as the corresponding r and & will be greater. 

Therefore, on this account, the right-hand member of equation (1’) will tend to be 
positive; but, furthermore, we are still left with a region of decreasing (3I' — 4)Ppo at 
the center. Thus it is very unlikely that this case could lead to instability. 

If (31 — 4)Po starts by increasing at the center, then we shall have as many regions 
in which (31 — 4)Po increases as regions where it decreases. The same reasoning would 
apply, but in this case the argument loses some of its force, as we have no extra region 
of decreasing (3I' — 4)Po near the center. Thus cases in which an important ionization 
takes place just at the center might require a more careful analysis. 

However, for the time being, we shall assume that the fundamental mode is the most 
dangerous for the stability of the star and that it is stable if I is everywhere greater 
than #. 
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2. Approximate formula for o?.—Now we shall consider the case in which I becomes 
smaller than 4 in some part of the star. For the fundamental mode of pulsation one 
usually gets quite a good approximation? for o” by assuming &é to be a constant in formula 
(5). In that case, equation (5) becomes 


% 


2 = 
T (6) 
where 
IPodVo 
= (7) 
PodVo 


Thus in a first approximation the sign of o? and consequently the stability or instabil- 
ity of the star will be determined by the value of I’, averaged with respect to I "Pod Vo 


over the whole star. As the pressure P decreases rapidly with increasing r, the external 
layers, where I is most likely to become smaller than 4, will have to be very extensive 
to influence I’ appreciably. 

3. Applications of formula (6).—As examples we shall consider the cases of the homo- 
geneous and standard models, and for the sake of simplicity we shall further divide the 
star into two parts: a central region, where I; = %, and an outer part, where [’, = 1, an 
extreme case very favorable to instability. 

If r. is the radius of the sphere separating the two regions, we know that, for r, = R, 
o” is positive and the star is stable; for r. = 0, o” is negative and the star is unstable. As 
o” is an eigen-value of equation (2), which is of the Sturm-Liouville type, it will vary con- 
tinuously with the coefficients of equation (2), and when r, decreases from R to zero, 
o” will decrease continuously, passing through the value zero and becoming negative. 
Thus o? will have one zero in the interval 0 < r < R, which will separate the stable from 
the unstable configurations. 

In a first approximation we can determine the critical value r, corresponding to o? = 
0 from a formula of the type (6), which in this case can be written as 


R 
or: PydVot orf PydVot+4Q. (8) 
0 


In the case of the homogeneous model, 


% = — 1$ and 


and our condition becomes, after introducing the numerical values of I; and I’., 


This equation has one root in the interval 0 < r./R < 1. We find 


or, in terms of mass or temperature, 


—~0.26 and —~0.6 


if 7) is the central temperature (small masses). 
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For the standard model 29 = — 3GM?/2R, and in Emden’s variables the condition is 
6.8969 
GY.= +5 uistds = (10) 


We find that equation (10) leads to the following values: 


Thus in both cases we see that, even with a I’ as small as 1 in the outer part, this must 
be very extensive before the stability of the star is endangered; in fact, the outer part 
contains appreciably more than half the mass. 

However, our approximation = constant, on which equation (8) is based, is rather 
crude. If we go back to equation (5), we see that if ¢ increases with r, the mean value of 
I’ should, in fact, be smaller than the one defined by equation (7), and this would in- 
crease the instability of the star. But, on the other hand, the third term of the numerator 
would be positive and contribute to the stability. 

To gain a definite idea as to the magnitude of these effects a rigorous method of treat- 
ing the problem will be developed in the following section and applied to the two cases 
considered here. 

4. Rigorous treatment in the case of a discontinuity of !'.—We shall consider again a 
star composed of two parts, separated by a surface of discontinuity of I’ (sphere r,), and 
we shall distinguish by the suffixes 7 and e the values relating to the internal and the 
external parts, respectively. 

We have to solve equation (2); but in this case, apart from the boundary conditions 
(3) and (4), we shall have two more conditions at the surface of discontinuity of T': a 
kinematical condition which reduces to 


A (&;) =B (&.) Te (11) 
and a dynamical condition, 
(6P;) (6P.) 
or, explicitly, 
dé, 


(12) 


= 
dr 


= Br, (3 


where £; and &, already satisfy equations (3) and (4), respectively. 
The condition that the homogeneous system (11) and (12) admits solutions other than 
the trivial ones A = B = O is that its determinant vanish. We must have 


ri (set re (38 


As £; and &, are functions of 0”, equation (13) provides us with an equation to determine 
its value. However, in general, equation (2) does not admit of analytical solutions which 
are explicit in 

The direct procedure would then be for a given value of r, to choose a value of o?, 
compute by numerical integration the corresponding values of £; and &, and determine 
whether they satisfy condition (13). Repeat this until a correct value of o” has been 
found. If it is positive, move 7, toward the center and start all over again, and so on 


(13) 


dé. 
d 


* 
j 
= 0 
+r 
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until a critical value of r., r’, is found such that, for r, > r2, o? > 0, and, for r. < ri, 
a? < 0. This would be a very arduous task. 

However, as we are mainly interested in the critical value of r,, we can obtain it much 
more simply by putting o? = 0 in equation (2) and treating r, as a parameter, which can 
then be determined by equation (13). In this way only one numerical integration will be 
required. 

a) The homogeneous model.—Writing x = r/R and using relation (9), we can write 
equation (2) as? 

4—6x? dé 
dx? dx 


3a? 4 
and a=(3-5). 


+Jé=0, (14) 
where 


If we suppose that o? = 0, then 
J=—2a. (15) 


Equation (i4) has two regular singularities, one at r = 0 (x = 0) and one at r = 
R (x = 1). The roots of the indicial equation at x = 0 are 0 and —3. The general 
solution would then have the form 


AX, (x) + A’x 3 [Gi +d; (x) lg x], (16) 


where \;(x) is holomorphic and ¢; regular in the vicinity of x = 0. As the solution (16) 
must satisfy the boundary condition (3), we must take A’ = 0, and we are left with 


Adi (x). 


T. E. Sterne, in the paper we have already referred to, has shown that \,(«) is of the 
form 


di (x) = ana, (17) 
k=0 


where the coefficients are determined in accordance with the relations 


2k (2k +5) —J 


= (2k+2)(2k+5)° (18) 


a=1 and 


If I’ is a constant over the whole star, then this solution should also satisfy condition 
(4), which is possible only if J has one of the values 


Ji, = 2k (2kR+5) ; (19) 


for the series (17) will then terminate with the term a2,x** and will represent the eigen- 
function of order k; the corresponding eigen-value of o” can be deduced from equation 
(19) and the definition of J. For any other value of J than those given by equation (19) 
the series does not terminate. However, it is converging for 0 < x < 1 and diverging 
only at x = 1. 

Another point brought out clearly by Sterne’s analysis is that, while the fundamental 
mode becomes unstable for ! < 4, the higher harmonics continue to be stable. Indeed, 
as we may directly verify from equation (19), the first harmonic becomes unstable only 
if ! < 2 and the second one if T < 34. 


7T. E. Sterne, M.N., 97, 582, 1937. 
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At r = R, the indicial equation has the double root 0. We can again, from the bound- 
ary condition (4), show that the solution of equation (14) will be of the form & = BX,(x); 
and, expanding \,(*) as a series in powers of (1 — x), we obtain from equation (14) 


d(x) = (1—-)*, (20) 


k=0 


with the following recurrence formula for the coefficients: 
= 3 [(k — 1) +4) — J] — + (R41) — 1) + (21) 


If we take 
J = (k—1)(k+4) (22) 


and if we admit that this value of J also satisfies the relation 
(kR—1) (23) 


so that b, = 0, it is easy to prove that d,4: and all the succeeding coefficients also vanish 
and the series will terminate with the term d,-:(1 — x)*". 

However, as we can readily verify with the first few terms, the relation (23) will be 
satisfied only if k is odd. The series will accordingly terminate with an even power of 
(1 — x). If Iisa constant for the whole star, our solution (20) must also satisfy equation 
(3), and this is possible only if the series terminates, that is, for the values of J given by 
equation (22) for k odd. One verifies that these J’s are identical with those given by 
equation (19). Thus, starting from the surface, we recover the same eigen-values as 
those given by Sterne. This is, of course, what we should expect. For any value of J 
different from the ones given by equation (22), the series (20) converges for 0 < x < 1; 
but it diverges at x = 0. 

In our case the values of J are fixed by equation (15). With the values I; = , 
I’, = 1, which we have already adopted, J; = —1.2 and J, = +2, and the solutions 
£; and &, will not terminate. However, as our condition (13) has to be applied at a point 
0 < x < 1, we need not be concerned with the possible divergence of the series we have 
referred to above. 

Computing the numerical values of the coefficients a2, and 5; for this particular case 
and introducing the corresponding solutions (17) and (20) into equation (13), we obtain 
the following values of x, or r-/R: 


2=0.575, (x.)4=0.717, (x.)¢= 0.725 , = 0.727, 


which correspond, respectively, to approximations limited to terms of degree 2, 4, 6, 
and 8. These values converge so very rapidly that we may adopt the last value as 
precise enough for our purpose. 

Comparing this value 


Me _ Te 
i 0.727 or —= 0.384 or T= 0.47 


with those given by our first approximation, we see that the instability is somewhat 
greater than is disclosed by equation (6). This corresponds to the fact, illustrated in 
Figure 1, that & increases slightly from the center to the surface, d§/dr experiencing a 
discontinuity when the compressibility changes. 

b) The standard model.—The case of the standard model can be treated in the same 
way, except that we cannot obtain series converging all over the interior of the star. But, 
before discussing this problem, we shall first determine whether in this case also the gen- 
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eral proposition on the stability of the star toward its different modes is satisfied. The 
exact verification would require numerical integrations. However, we can avail our- 
selves of the method developed in a previous paper.* It was shown there that successive 
approximations for a” can be obtained by solving determinants of order 2, 4, 6, etc. Be- 
cause of the form given there to those determinants it is evident that in a given approxi- 
mation (determinant of rank 27) the condition for one or more of the roots o? to be zero 
is that the minor of rank 7 in the lower left corner be zero. This provides an equation 
ina = (3 — 4/T), which enables us to find successive approximations for the critical I”. 

In this way we obtain for the first harmonic the critical value T'} = 0.891 in a first ap- 
proximation and I} = 0.922 in a second approximation. For the second harmonic, the 
first approximation gives I’? = 0.537. Of course, this does not provide upper limits for 
the critical I’,; but, as the method converges fairly rapidly as far as the values of o? are 


res/s ret 


Qu o4 as or ae os 


Fic. 1.—Variation of the amplitude with the radius in the homogeneous model for a discontinuity 
of T as indicated. 


concerned, we can assume that the same is true of the I’, and that their order at least is 
correct: ¢ > I! > 2. In the case of I} it is even safe to assume that its real value 
will be smaller than 1. 

Returning to equation (2), we can re-write it in terms of Emden’s variables for the 
case o? = 0 in the form 


dz dzlz u dz uaz 


where a has the value a; = 0.6 for z < 2, and a, = —1 for z > z. With the help of 
well-known series for u, we can obtain series for £; and &, in the neighborhood of r = 0 
and r = R, satisfying conditions (3) and (4), respectively. With these we can start our 
numerical integrations, one at the center and one at the surface. 

The necessary calculations are rather light, as it turns out that the solution which we 
have to extend furthest is &. This is the simplest, too, as & and dé,/dz vary very slowly 
in a large part of the star. The simplest way to solve equation (13) in this case seems to 
be to compute its value at different points and then interpolate between its positive and 
negative values. This shows that the critical value of z, is of the order of 2.45. 

Figure 2 represents the corresponding solution £, which exhibits the same peculiari- 
ties as the ones discussed in the case of the homogeneous model. However, here & in- 


8 P. Ledoux and C. L. Pekeris, Ap. J., 94, 124, 1941. 
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creases more rapidly, as should have been expected from the greater central condensa- 
tion of the standard model. And, owing to this, the correct values, 


FF~0.66, and 720.47, 


deviate more from the first approximation than in the case of the homogeneous model. 

Of course, a varying I! may mean also a varying mean molecular weight, uz, which 
would result in a slightly different distribution of density. Also the small value of T 
adopted in the external layers would lead to convection, which would again affect the 
distribution of mass. However, the results just obtained for two models with widely dif- 
ferent central condensations show that small deviations, such as those referred to above, 


Fic. 2.—Variation of the amplitude with the radius in the standard model for a discontinuity of T as 
indicated. 


would affect very little the critical values, especially the critical temperature T,. Thus 
we can conclude that, even if I’, has its smallest possible value [, = 1 in the external 
layers, these layers have to extend to a depth at which the temperature is approximate- 
ly half the central temperature before the star becomes unstable. 

5. Extension of these results to cases in which T, ¥ 1.—Of course, I’, in practice will 
not reach such a small value as I’, = 1; and Figure 3 (full curve) shows how the critical 
ratio, r./R, as given by our approximate formula (10), varies with I’, for the standard 
model. 

We know the correct critical values for I’, = 1, and another integration in the case of 
I’, = 1.2 gives 


0.71, 
To 


On the other hand, as I, tends toward #4, r./R tends toward zero, and £ tends to become 
a constant, so that the correct critical values will approach more and more the approxi- 
mate ones and their variation with I’, cannot differ very much from the representation 
given by the dashed curve. For I’, = 1.3, the external layers have to reach a depth at 
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which the temperature is of the order of 88 per cent of the central temperature (r./R = 
0.13; m./M = 0.09) before the star becomes unstable. 

6. Extension to cases in which T; ~ 3.—On the other hand, until now we have al- 
ways supposed that I’; = %, which corresponds to a negligible radiation pressure, 
[(M/Mo)# < 1]. If the mass increases, I’; will decrease and tend toward 4 for very 
large masses. 

With I, = 1, our approximate equation (10) gives the results summarized in Table 1. 


ao i 


19 i 12 13 43 


Fic. 3.—Variation of the critical ratio r,/R with I, in the case of the standard model and I; = 5/3 


TABLE 1 
(M/Mo©)a* ry re/R | m/M | T./To (M/Mo)a* ry | | T-/To 
5/3 0.23 | 0.35] 0.7 1.379 | 0.41} 0.78} 0.38 
1.538 27 .46 62 1.366 43 
1.481 .29 1.350 | 0.49} 0.88; 0.30 
1.429 | 0.33} 0.60] 0.51 


These results are also plotted in Figure 4 (full curve). Of course, the real critical values 
will be different; but we know the correct value for '; = §, and another integration for 
I’; = 1.363 gives 


m 

T, 29, and ~0.89, 
Again the correct and approximate curves should approach one another as I; tends 
toward 4, so that we can again tentatively draw a dashed curve, as in Figure 4, to repre- 
sent the real variations of r,/R. 

Thus, even for large masses, the external layers, where I’, is as small as 1, must still 
be very extensive. For instance, if (M/Mo)n? = 400, we find from the dashed curve in 
Figure 4 that I’, must be equal to 1 as far as 


r.~0.48R , where 


before the star reaches the limit of stability. 
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7. Small Y's in the vicinity of the center.—In contrast to our discussion in the previous 
sections, we can contemplate the opposite situation and suppose that T; < 4and TI, > #. 
If we take I; = 1 and I, = #, the two terms in equation (10) simply change sign, and 
thus in a first approximation we still get the critical values 


r,/R&0.23 , m,/M=0.35, and T7,/T~0.7. 


As the integrand of the second term in equation (10) has a well-marked maximum 
around z = 1.4, the I’s corresponding to that region will have the greatest weight in 
the evaluation of I’, and therefore it is in that region that a small I’ could have the greater 
effect. 

If we take I = 1 in an interval between 2; and z, having z = 1.4 as its mid-point and 
I = 3 everywhere else, we find that the star is on the verge of instability in a first ap- 


473 14 15 1s 


Fic. 4.—Variation of the critical ratio r./R with T; in the case of the standard model and I, = 1 


proximation, when 2; & 0.9 and z: & 1.9. Thus for instability the region of ! = 1 must 
extend from r/R ~ 0.13 to r/R ~ 0.28; accordingly, it occupies a spherical shell of 
thickness R/7 and has a mass of about 0.4 M; the temperatures at its boundaries are, 
respectively, 88 and 61 per cent of the central temperature. 

8. Numerical applications—No detailed application has been carried out, but often 
the results given in the preceding sections enable one to decide whether a star is stable 
or not. In general, it appears that the conclusions of L. Biermann and T. G. Cowling? 
will be confirmed. 

As an example we shall treat the case of the sun. R. H. Fowler and E. A. Guggen- 
heim‘ have shown that, quite generally in the interior of a star, only two consecutive 
states of ionization of a given element need be considered at the same time. For a given 
temperature 7 and electron concentration N,, there will be approximately as many 
atoms in one state as in the other if the difference of ionization energy between them is 


of the order of y, as defined by 
2 (24m_kT) 


Under these circumstances we can also assume, following Biermann, that all the states 
with ionization potentials x < y — kT will be ionized and those with ionization poten- 
tial x < ~+ kT will be un-ionized. 
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For the sun the standard model gives a fairly good representation, and, for an abun- 
dance of hydrogen by weight X of the order of 4, it leads to a central temperature of 
2 X 107 degrees and a central density of the order of 75 gm/cm*. Under these conditions 
the ionization will already be well advanced, and in a first approximation one can write 


(= (25) 


With this value of V., equation (24) gives for y a value of the order of 4400 ev. at 
the center of the sun. Looking at a table of ionization potentials, one realizes immediately 
that no important ionization will set in here. Thus, near the center of the sun, the matter 
can be considered as behaving like a monatomic gas, and, as the radiation pressure can 
be neglected, we can take T = 8. 

At the point where the temperature drops to half, its central value y is still of the order 
of 3300 ev., and we can again take = 3. But we should require I to be equal to 1 here 
and in the rest of the star to bring it to the verge of instability, and this is obviously 
impossible. 

Even for a vanishing abundance of hydrogen, if we take the following mixture to 
represent the relative abundance by number of atoms of the other elements: O:Mg:Fe = 
8:3:1, the same considerations as before show that the sun would still be stable. 

These results can easily be extended to the other stars of the main sequence. However, 
for stars of very large masses the pressure of radiation becomes important; and, as that 
case was not considered by Biermann and Cowling, we shall treat an example. Let us 
take an extreme case—the Trumpler stars—and let us consider especially NGC 
6871.5, for which M ~ 400 Mo and R & 16.6 Ro. If we suppose that it is built on the 
standard model and composed of pure hydrogen, the central temperature, 79, is of the 
order of 10° degrees and the central density of the order of 6.6 gm/cm*. The ratio 6 of 
the gas pressure to the total pressure is 0.375. Under these conditions, since hydrogen is 
very easily ionized, we can assume that, in an extensive region starting from the center, 
the adiabatic exponent for the gas is } and the combined exponent I for the matter and 
radiation will be 1.4. From Figure 4 we see, then, that for instability I should be equal 
to 1 in a region extending from the surface to r > 0.43R, where T © 37 o. But at this 
temperature the hydrogen is still completely ionized, and the region where I is equal to 
1.4 will extend much farther toward the surface. Thus the star is stable. 

If the star does not contain any hydrogen, the central conditions remain more or less 
the same, as the increase of yu is practically compensated by the decrease of 8, which 
becomes 6 = 0.1. At the center the gas can again be considered as monatomic, and 
I’ = 1.35. With this value of I in the central part we see from Figure 4 that, for instabil- 
ity, I should be equal to 1 from the surface to a point r © 0.54R, where T & (})7o. 
But, for that temperature and the corresponding density, y is still of the order of 
20,000 ev., and I will retain its value 1.35 much further, and the star is still stable. 

However, it is known from other evidence’ that the standard model gives a very poor 
representation of the internal structure of the Trumpler stars and a homogeneous model 
probably would be better. The problem is then a little complicated by the fact that, 
for such a model, 8 varies with depth; but for pure hydrogen one can again determine 
that the star is stable. 

For vanishing hydrogen content, the central temperature is of the order of 3.2 X 107 
degrees, B ~ 0.06, and y ~ 28,000 ev., so that we can still take for I’ the value cor- 
responding to a mixture of monatomic gas and radiation. We find T ~ 1.343. At the 
point where T = $7) and r © 0.96R, y is still of the order of 12,500 ev., and for the 
mixture considered above we can still neglect ionization and T = 1.4. Computing T at 


® Chandrasekhar, of. cit., p. 313. 
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different points inside that region, we find that it remains very close to its central 
value, as far as r ~ 0.8R, so that the average taken with respect to f *PdV cannot be 


very different from it, say ! = 1.35. Using equation (6), we find that [ must be smaller 
than 1 in the external region to reach instability. However, we know that the instability 
is somewhat greater than revealed by equation (6); and I’, although greater than 1, 
will be fairly small in that external region, so that probably we could not have much 
greater masses built on such a model except if they contain an appreciable amount of 
hydrogen. 

In that respect it seems that the condition of vibrational stability is more restrictive. 
It has been shown,’ for example, that a star built on the standard model, for which 
Kramer’s law of opacity and Bethe’s law of generation of energy are valid, would be- 
come unstable when the quantity (M/Mo)j,? is somewhat greater than 100. Since in 
the case of NGC 6871.5, M/Mog is of the order of 400, we should have to take an ex- 
tremely large abundance of hydrogen to avoid instability. 

A change of model in the sense of a greater homogeneity would not help, since in that 
case the amplitude & would increase less rapidly from the center to the surface and the 
stabilizing terms which arise near the surface would have less weight in the integral ex- 
pressing the condition of vibrational stability. 

For stars of very small masses falling. in the region of low hydrogen content in the 
Hertzsprung-Russell diagram, it does not sein either that any instability of the kind 
considered here will appear. For instance, in the case considered by Biermann" in connec- 
tion with his theory of the nova phenomena, M ~ 0.5Mo, R = }Ro, the radiative 
equilibrium might become unstable from the surface to a point where the temperature is of 
the order of 10°-° degrees to 10’ degrees. But this is of the order of only one-tenth of the 
central temperature, which in this case is 7‘) ~ 10° degrees. From this point to the point 
where T = $7», I is certainly greater than 1. Thus, although the radiative equilibrium 
can be unstable fairly deep, the star remains dynamically stable, at least for a radial 
perturbation. 

Of course, as Biermann and Cowling? have shown, it is for large radius that dynamical 
instability appears most easily; and past a certain value of the radius (for a given mass) 
one can determine the minimum abundance of hydrogen necessary to keep the con- 
figuration stable. The method developed here could be used to obtain more precise values 
of the critical radius or of the minimum amount of hydrogen, but it would require de- 
tailed computation. 

9. Special cases.—This discussion has left out some special cases, but it is doubtful 
whether they have any physical interest. For instance, we limited ourselves to real 
values of o?. But in the case of the Roche’s model, T. E. Sterne’ has shown that all modes 
are unstable for [' < 4, all the corresponding o? then having an imaginary part. 

As the coefficients of equation (2) are essentially real, o? can have only an imaginary 
part if the same is true of & But we can study & in the vicinity of the singularities of 
equation (2) and see whether, for models having a physical meaning, the appearance of 
imaginary values is possible. For such models p, P, and I remain finite and different from 
zero for 0 < r < R. We can also suppose that dI'/dr remains finite in the same interval. 

Then equation (2) can have singularities only at r= 0 andr = R. Atr=0, y= 
Gm/(r)p/Pr tends toward zero as r? and (1/T)(dI'/dr) remains finite. The coefficient of 
dé/dr will therefore tend toward infinity as 4/r. In the same way the coefficient of & will 
tend toward infinity as 1/r. Thus r = 0 is a regular singularity, and the corresponding 
indicial equation, 0(@ — 1) + 40 = 0, has real roots only, and the solution will be 
essentially real near r = 0. At r = R, the singularities of equation (2) will depend on 


10 P. Ledoux, Ap. J., 94, 537, 1941. 
11 Op. cit., p. 344. 
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the behavior of the ratio p/P and the quantity dI'/dr. As long as we use equation (2), it 
seems natural in this connection to take the same point of view as the one adopted in the 
theories of the internal structure of stars, namely, that P, p, and 7 tend toward zero 
at r = Rand that Kramer’s law of opacity is valid. 

We may then verify” that p/P will tend toward infinity as 8/(1 — r/R) and that 
dI/dr, which ultimately varies as d8/dr, remains finite. Thus r = R is also a regular 
singularity, and the indicial equation, 6(@ — 1) + 8@ = 0, can have only real roots. 
Thus & will be real everywhere, and so will o. 

Finally, if a star approaches dynamical instability (« becomes small), there will be a 
point at which the usual method of perturbation used to obtain the condition of vibra- 
tional stability will cease to be applicable, since the perturbation o’ of o (due to the non- 
adiabatic processes) will no longer be small compared to o and the terms in o” will not 
be negligible compared to go’ or o?. 

It would probably be interesting to study the interaction of these two types of instabil- 
ity and try to obtain a more general criterion of stability. 


2S. Chandrasekhar, M.N., 96, 647, 1936; also J. Tuominen, Ann. Acad. Sci. Fenn., Ser. A, Vol. 48, 
No. 16, 1938. 
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ABSTRACT 


Line intensities measured on Mount Wilson coudé plates of 10 Lacertae provide data for an investiga- 
tion of the atmosphere of this star by Unséld’s method. From the Stark broadening of the hydrogen lines 
and the ionization equilibrium, the electron pressure, temperature, and surface gravity are estimated: 
log P. = 2.80; effective temperature JT = 29,600° K; and log g = 4.44. Within the large errors of the 
data, the relative abundances of helium, carbon, nitrogen, oxygen, neon, magnesium, and silicon turn 
out to be the same as in 7 Scorpii. The discrepancy in the computed neon to oxygen ratio between the 
planetary nebulae and the hot stars remains. Calculations of target areas for collisional excitation of for- 
bidden neon lines, on the one hand, and f-values for the permitted transitions observed in the stars, on 


the other, are necessary before this question can be clarified. 


In a discussion of the continuous spectrum it is shown that the need for a blanketing-effect correction 
of the type indicated by Unsdld is almost eliminated if the mean value of the stellar absorption coefficient 
is computed in the manner suggested by Chandrasekhar. 


INTRODUCTION 


In the course of a photometric investigation of the spectra of the brighter nuclei of 
the planetary nebulae it became apparent that a representative early-type absorption- 
line star could serve as a useful comparison object. The electron pressure, effective tem- 
perature, and composition of main-sequence stars of about the same temperature as the 
central stars of the planetaries might provide useful norms for the possibly more com- 
plex atmospheres of the latter objects. At least, comparisons between absorption-line O 
stars and planetary nuclei, on the one hand, and characteristic Wolf-Rayet stars and 
Wolf-Rayet type nuclei, on the other, may be valuable. 

Unséld’s study of 7 Scorpii' provides us with useful data of the type we are seeking, 
but it seemed worth while to apply similar methods to at least one other representative 
main-sequence star. Accordingly, for comparison with the planetary nuclei I have 
chosen 10 Lacertae, an O-type, main-sequence object, which falls close to the tempera- 


ture and level of ionization of representative planetary nuclei. 


THE OBSERVATIONAL MATERIAL 


The data upon which the present discussion is based consist of coudé plates, obtained 
with the 100-inch reflector at Mount Wilson Observatory. The suitable spectra are given 
in the accompanying tabulation. 


Ex- Wave-Length 
Plate Date Emulsion Disp. Camera Range 
posure 
(AA) 
Ce: S62 vices 1937 Sept. 7 120" | Imp. Ec. Sp. Prism “ye 4900-3500 
KOS) rer 1945 Aug, 25 130 East. 103a-0 Grating 114-inch 4140-3500 
4063... ..«: 1945 Nov. 22 212 East. 100a-0 Grating 114-inch 4140-3500 


* This investigation was made possible by a co-operative arrangement with the Mount Wilson Ob- 


Servatory. 


1 Zs. f. Ap. 21, 22, 1941. 
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Table 1 gives the wave lengths adopted from Struve’s measures? and my own, the 
identifications, and the quantity log W/d + 6, where W is the equivalent width and \ 
is the wave length. Log W/) is the physically significant quantity for curve-of-growth 
discussions. 

For lines to the violet of \ 4143 the tabulated intensities are based upon all three 
plates; for the longer wave lengths only the Ce 952 data are available. For many of the 
weaker lines the intensities are very uncertain, and errors of a factor of 2 in W are quite 
possible. Despite their much greater dispersion, the grating plates show no more lines 
than the prism plates. A comparison of the hydrogen lines on the prism and the grating 
plates shows that a correction of about 20 per cent for light scattered by the grating must 
be applied. The agreement between the corrected intensities derived from the grating 
plates with those found from the prism plates is satisfactory, considering the weakness of 
the lines. Many of the lines lie just at the limit of measurability on the coudé plates. 


TABLE 1 
LINE INTENSITIES IN THE SPECTRUM OF 10 LACERTAE 
Identi- log Identi- log Identi- log 

fication | W/A+6 fication | W/A+6 fication | W/A+6 
3587 Het 1.62 || 3998.66...| Nur, Su? | 0.99 || 4275.67... OU 0.54 
3613.54... Het 1.30 || 4003.47... Niu 1.32 || 4303.88... Ou 1.25 
$634.27... Het 1.61 || 4009.36... Het 1548 Om 0.20 
3664.0.... New 1.15 |} 4026.22... Het 2:37 4325.64...) Ou, Cm 0.08 
3694.24... New 1.32 || 4056.09... Cm 1.15 || 4345.60... Ou 1.25 
3695.4.... Om 0.74 || 4060.70... On 1-20) 1.19 
3699.15... Om 0.92 || 4062.9.... Ou 0.70 || 4349.42... Ou 1.42 
3703 .4.... Om 1.18 || 4067.93... C Ir 2231: 4558.27... OU 0.86 
3705.10... Het 1.70 || 4069.11...) Cu, Cm 1.32 || 4366.9.... Ou 0.86 
3707.22... O 1.09 || 4069.81...} On, 0.53 
3709.68... Nen 1.14 |} 4072.18... OU 1.32 || 4388.0.... Het 1.97 
3713.03...| New, Om 1.37 || 4074.00... Our 1.04 || 4392.0.... S11 1.02 
3715.05... O ul 1.20 || 4075.90... Ou 1.38 || 4437.48... Het 0.98 
3727 .33...| Neu, Om 1.52 || 4078.86... Ou 1.02 || 4465.5.... OU 0.87 
3750.25... Ou 0.76 || 4081.01... Our 0.90 || 4471.5.... Het 2.248 
3754.76...| Nim, Or} 1.54 |} 4083.89... Ou 0.80 || 4481.25... Mgtt 1.04 
3757.26... Our 1.54 || 4085.20... Ou 0.79 || 4510.94... Ni 1.2:* 
3759.95... Our 1.71 || 4088.95... SiIv 1.82 || 4514.84... Niu 1,53 
3762.36...| Om, Sitv 1.35 || 4097.32 Nia 0.64 
3766.46... New 1.10 || 4103.40 Niu 1.35 || 4925-5;... Niu 1.19 
SiIv 0.91 || 4116.11 1.62 || 4541.75... Heu 1.92 
3774.0.... Our 1.22 || 4119.25 Ou 1.25 || 4552.68... Si tr 0.54 
3776.98...| Neu, Om 1.10 |} 4120.82 Het 1.71 || 4631.18... SiIv 1.16 
3791.40...| Om, Simm] 1.44 || 4143.84 Hei 1.70 || 4634.28... Ni 0.90 
3803.04... Ou 0.80 || 4152.54... C Ir 1.19 || 4638.62... Ou 1.05 
3805.55... Het 4:02 Su? 1.01 || 4641.61...) Nim, Om 1 25 
3806.56... Si U1 0.95 |} 4156.51...] On, Cm 1.26 |! 4647.42... C Ir 1.67 
3819.75... Het 2.21 || 4163.00... C Ir 1.15 || 4650.24...) Cm, Om? | 1.75 
3867.42... Het 1.30 || 4169.04...} Her, Om 1.04 || 4651.46... C 1.30 
Ou 1:05: Om 0.78 || 4654.41... SiIv 1.30 
3911.99... Ou 0.75 || 4186.91 Cm 0.58 
3926.74... Het 1.07 || 4189.85... Om 1.02 || 4661.85... Ou 0.87 
3938.82... N ut 0.92 || 4195.77... Nim 1.12 || 4685.80... Het 2.22 
3954.37... ou 0.80 || 4200...... Heu, Nur} 1.49 || 4705.32... Ou 1.12 
3961.67... Om 1.45 || 4212.38... SiIv Het 1.76 
3964.88... Het 1.56 || 4253.80...] Om, S m1? 1.28 || 4921.93... Het 245 
3995.00... Nu 0.83 || 4267.19... Ci 1.28 


* The colon denotes a line intensity possibly vitiated by a plate defect. 


2 Ap. J., 90, 722, 1939; see also ibid., 74, 225, 1931; and 77, 321, 1933. 
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The discussions of curve of growth and abundances were carried out some time ago on 
the basis of the results from Ce 952. Since grating intensities are available for only a 
portion of the total wave-length range and since the prism intensity system agrees well 
with the corrected grating intensity system, I have not recomputed the abundances of 
helium and the heavier elements. Although a few improvements might be obtained, it 
would seem appropriate to hold new calculations in abeyance until additional homo- 
geneous data are available for the entire spectral range. On the other hand, since the 
confluence of the Balmer series plays an important role in the calculation of the hydro- 
gen abundance, all data have been utilized for this problem. 

In Table 2 the second column gives the central depth of the lines measured in terms 
of the continuous spectrum as unity, the third gives the logarithm of the equivalent 


TABLE 2 
THE HYDROGEN LINES IN 10 LACERTAE 
Wave Length log W log NoH Wave Length log W log 
0.41 0.39 13.99 0.40 0.34 15.47 


width in Angstrom units, and the computed quantity in the last column (see eq. [2]) en- 
ables us to estimate the population of the second level in hydrogen. The central depths 
of the hydrogen lines are not given beyond \ 3835, where the lines begin to overlap 
seriously, and the true position of the continuum is unknown. The material suggests 
that the true central intensities of the hydrogen lines are nearly constant, and I have 
adopted R, = 0.47 as the most probable value. 


INTERPRETATION OF THE LINE INTENSITIES 


For the estimates of electron pressure, temperature, and composition, I have followed 
very closely the method given by Unsdld in his r Scorpii paper. Except for O 11, where 
additional line strengths had to be computed, Unsdéld’s f-values have been used through- 
out. Both 7 Scorpii and 10 Lacertae are main-sequence stars of apparently similar at- 
mospheric structure, and the procedure should suffice to fix the relative compositions 
of the atmospheres of these two stars, even though the absolute values might be in error. 
For the purposes of the initial calculations, the assumed temperatures of 31,500° K was 
taken from Kuiper’s scale for an 09 star. Hence we need to correct Unsdld’s log f/(R-Awp) 
by only —0.07 in order to apply them to 10 Lacertae. Here f is the Ladenburg “f” or 
oscillator strength; the Doppler breadth Awp is given by 


./2kT 
where T is the kinetic temperature and M is the mass of the atom in question. 

The total number of hydrogen atoms in the second level above the photosphere, 
No»2H, may be estimated in three ways:? 

3 The ordinary curve-of-growth procedure is not applicable to H, because the effective damping con- 
stant, which is determined by the interatomic Stark effect, increases from line to line in the higher mem- 
bers of the series. Hence other procedures for estimating abundances must be employed. For a discussion 
of the problem of the hydrogen lines in stellar spectra see Unséld, Physik der Sternatmosphdren, p. 286; 
see also Merrill and Wilson, Ap. J., 80, 19, 1934. 
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1. On the assumption that the hydrogen lines are formed in an optically thin layer 
(No2Hf<1), the equivalent width, W,, for the 2—n transition is related to No2H by 
the expression 


2 
= 0.886 No. 2H. (2) 


Among the lower members of the series, this relation is certainly not correct, and the 
computed No 2H will be toosmall. Equation (2) approaches validity as we go to higher and 
higher series members. Near the limit of the series, another difficulty comes in. Because of 
the overlapping of the line wings, one tends to draw the continuum too low, and the esti- 
mates of W,, are consistently too small. A plot of No.H against m gives a curve that 
rises for a time and then levels off and finally dips down as the measured W,, becomes 
too small. If we extrapolate the rising branch of the curve to a large value of n, in the 
fashion indicated by Unsdld, we obtain the limiting value for an optically thin layer. In 
this manner I have obtained log No.2H = 15.72. 

2. One may employ the value of the Balmer discontinuity as measured by Barbier 
and Chalonge‘ in conjunction with the ‘thin-layer approximation” to estimate the 
number of hydrogen atoms in the second hydrogen level. Since log D = 0.04, we obtain 


1.38 X10-"No,.H = 0.09 , (3) 
whence log No2N = 15.81. 

3. Finally, we may treat the dip in the continuum at the Balmer limit by the inter- 
polation formula more commonly employed in the “‘finite-layer”’ treatment of absorption 
lines.5 

1 


Here R is the amount of the depression below the level of the undisturbed continuum, 
ky = 1.38 XK 10-'7No2H is the absorption coefficient at the Balmer limit, and R, = 0.47 
for 10 Lacertae. Since R = 0.09 from the date of Barbier and Chalonge, log No2H = 
15.90. 

As the most probable value, I have adopted log No2H = 15.80, a result very similar 
to that obtained for 7 Scorpii by Unsdld. 

One may compute the electron density from the equivalent widths of 8, Hy, and Hé 
or estimate it from the number of resolvable hydrogen lines on a good plate (Inglis- 
Teller method). 

If N is the number of charged particles effective in broadening a line of equivalent 
width W,, then 


we? (5) 


12 
N = 6.5 X10" 


where C,, is a constant for each line. If, following Unsdld, we suppose that N = 1.5N,, 


we find 
4861 4340 4101 


Pee 13.88 14.18 14.34 
The mean value is log V, = 14.13. 


4 Ann. d. ap., Vol. 3, No. 2, 1940. 
5 See Physik der Sternatmosphiaren, eq. (69.3), p. 266. 
6 Ap. J., 90, 439, 1939. 
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The number of resolvable lines of the Balmer series depends on the Stark broadening 
and therefore on the ionic density. If the density of the ions is large, as in the atmosphere 
of a dwarf, the Stark broadening will be large and the resolved Balmer lines will be fewer 
than in a star of low atmospheric density. If ,, is the principal quantum number of the 


10 LACERTAE 


last resolvable Balmer line, the Inglis-Teller formula gives 


log N = 23.26—7.5 log nm « 


Struve has listed the Balmer lines up to m = 14 in 10 Lacertae; thus we obtain log V, = 
14.6. We cannot be sure, however, on the basis of the present data, that some kind of 
turbulence does not enter to widen the Balmer lines; hence I have adopted log V, = 


14.10 for further calculations. 


From the neutral helium lines and the thin-layer approximation, we can obtain a 
lower limit to the population of the 2°P level by the methods indicated by Unsdéld. In 


this way it is found that 


log N (28P) =14.5+0.3. 


With respect to ionized helium the situation is more favorable for our purposes in 
10 Lacertae than in 7 Scorpii, since we can measure the intensities of two of the Pickering 
lines, \ 4541 and \ 4200. If we wish to employ the thin-layer approximation (eq. [2]), 
as Unséld suggests, we must reduce the measured profile R(A) by means of equation (4) 
to obtain ky. From the computed &) one may then compute the corrected width by means 


of 


Here W, is the equivalent width reduced to a thin layer. This procedure gives us a lower 
limit to the NH computed from equation (2). The numerical results are given in the 
accompanying tabulation. Notice that the NV, calculated from the 4686 line (reduced 


Wi=fkdn. 


Lower Level nN log W log W’ log Nn’A log N4H 

4686 —0.11 +0.06 13.14 13.00 


from n’ = 3 to n’ = 4 with the Boltzmann equation and 7 = 31,500° K) is ten times 


smaller than that computed from the Pickering lines! 


It seems not unlikely that the 4686 line may be affected by incipient emission, while 
the Pickering lines are little disturbed. In the spectra of certain O-type stars the \ 4686 
series is sometimes seen in emission, while the Pickering series appears in absorption. 
Hence I have adopted log NH = 14.00 from the two Pickering lines. 

The calculation of the abundances of the other ions—C m1, Nu, N m1, On, Nenu, 
etc.—has been carried out in the manner described by Unsdld, so far as the 10 Lacertae 
data will permit. In this hotter star the lines are fewer and weaker, and the accidental 
errors in the intensities are more serious. In particular, the Om, Net, and N m1 lines 
are much weaker in 10 Lacertae, but the O m1 lines are stronger and more numerous. I 
have plotted log W/(R.2AXp) against log f/(R-Awp) in the usual way and have obtained 
the log NH values by a comparison with the theoretical curve of growth. The material 
does not suffice to fix the damping constant a; hence, on the basis of Unséld’s work, I 


have adopted log a = —1. 
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The columns of Table 3 list, first, the observed ion, the excitation potential of the 
lower term (or configuration) of the transition involved, subtracted from the ionization 
potential, viz., x, — x,,s; the logarithm of the number of atoms in the level s in the rth 
stage of ionization, log N,,,H; and, finally, log NV,41, HP, calculated from the combined 
Boltzmann and Saha equations for a temperature of 31,500° K (5040/7 = 0.16) and 
also 29,500° K (5040/T = 0.17). The lines of O11 enable one to compute the number 
of O11 ions; the O 11 lines, in their turn, enable the calculation of the number of O 1v 
ions to be made; etc. In other words, the lines arising from an atom in the rth stage of 
ionization enable us to compute the number of ions in the (r + 1)th ionization stage. 


TABLE 3 


DATA FOR THE CALCULATION OF IONIC ABUNDANCES 


L0G 
No. oF = ADOPTED 
Ton Lines* | Ny,sH 
6=0.16 6=0.17 
{ 7.64 1 11.50 20.98 20.80 
| 7.92 3 12.80 21.60 21.46 21.20 
18.25 3 13.40 21.40 21.15 
11.09 1 12.30 21.31 
{ 8.00 1 13.00 22.00 21.72 
8 13.10 21.10 20.93 21.00 
lie'99 1 12.00 19.67 19.42 
6.17 2 13.40 22:50 22.22 
9.25 8 13.40 22.10 21.94 22.00 
12.12 6 13.50 22.20 22.00 
{17 64 4 14.20 21.70 30.42} 20.80 
17.64 4 13.60 20.60 20.42 
{ 6.20 1 12.30 21.90 32.08) 21.95 
13.82 5 13.90 22.30 22.08 , 
Pe 6.14 1 11.88 20.97 20.84 
14.39 1 12.2 20.8 20.6 
{ 8.67 12.30 20.70 20.58) 20.70 
{21.00 2 14.00 21.40 


* The third column gives the number of good lines available for estimating the abundance of the ion 
concerned. Lines possibly blended are given half-weight. 


If we use T = 31,500° K in connection with the Saha equation and the computed 
ratios, He N 1m, O1v/Om, Sitv/Si m1, we obtain log P, = 3.20 and 
log N, = 14.56. On the other hand, a temperature of 29,600° K gives log V, = 14.10, 
in closer agreement with our previous estimate from the Stark broadening of the hydro- 
gen lines. Quite possibly the true electron density and effective temperature fall between 
these limits. 

For the purposes of an illustrative calculation I have adopted log P, = 2.80, T, = 
29,600° K. Actually, the relative abundances of the constituents are affected only very 
little by the exact choice one makes here. For example, an electron pressure of 1,600 
dynes and 7, = 31,500° K give essentially the same relative abundances. Factors to re- 
duce the observed ionic abundances to the relative proportions of the actual elements re- 
sult in small corrections to the observed ionic densities (Table 4). 
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The close agreement of the proportions of helium, carbon, nitrogen, oxygen, neon, 
silicon, and magnesium in the atmospheres of 10 Lacertae and 7 Scopii is probably 
fortuitous in view of the uncertainties accruing from the low line intensities in the hotter 
star. The apparently greater abundance of hydrogen in 10 Lacertae may be real, but it 
is more tempting to attribute it to the observational uncertainties, even though the 
hydrogen abundance should be better established than those of the other elements 
whose relative proportions agree so well with the 7 Scorpii results. Nevertheless, the lines 
of the heavier elements all seem to be weaker than they should be if the proportion of 
hydrogen to the total mixture is the same as in 7 Scorpii. 

In 10 Lacertae as in 7 Scorpii, the neon abundance turns out to be about the same as 
the oxygen abundance, in marked contrast to the results obtained from the planetary 
nebulae,’ where the abundance of neon was found to be between 0.1 and 0.01 that of 
oxygen on the basis of the intensities of the forbidden lines of Ne m1 and Ne v and of 
plausible assumptions concerning the target areas for collisional excitation. Until the 


TABLE 4 


COMPOSITION OF THE ATMOSPHERE OF 10 LACERTAE 


4 
RELATIVE NUMBERS OF ATOMS 
Atom toc NHPe 
10 Lac t Sco 

|! 25.33 1600 1000 
24.26 134 180 
22.02 0.8 1.0 
21.97 0.7 
20.95 0.06° 0.064 
20.84 0.05 0.06 


target areas are accurately computed, however, we cannot decide whether this apparent 
discrepancy is real or not. It may be that 10 Lacertae and 7 Scorpii are “neon-rich” ob- 
jects. An example in point in Nova Persei, which showed unusually strong forbidden 
neon lines. Further analyses of coudé spectra of early-type stars might help here. 

The figures given in Table 4 probably represent a lower limit to the ratio of the heavier 
elements to hydrogen. The f-values for many of the transition arrays may be appreciably 
less than: the value of 1, which Unséld had assumed. Hence the abundances of carbon, 
nitrogen, oxygen, etc., relative to hydrogen may have to be increased by a factor pos- 
sibly of the order of 2. An improvement in the abundance ratio can be attained only 
when the relevant f-values have been correctly computed. But it is difficult to see how 
the oxygen group could contribute as much as 35 per cent of the total mass of the stellar 
atmosphere, as Biermann,’ in a recent paper, supposes is the situation for the mean com- 
position of the sun. On the basis of the present results, an oxygen-group concentration 
in the neighborhood of 3 per cent seems quite possible; but it seems unlikely that the f- 
values could be off by the factor of 20, which the 35 per cent oxygen-group concentration 
would require. 


7 Aller and Menzel, Ap. J., 102, 239, 1945. 
Zs. f. Ap., 22, 244, 1943. 
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REMARKS ON THE THEORY OF THE CONTINUOUS SPECTRUM 


The problem of the continuous spectrum of 10 Lacertae and 7 Scorpii is essentially the 
same. Unsdld has discussed this question in detail. Since the lines appear to be formed in 
accordance with the mechanism of local thermodynamic equlibrium or “absorption,” 
the effective optical depth of the photosphere is 


1 
T+ 7) 
where 
3 he 
8 XRT k 
Bo = 1 — By (8) 
and 
k 
y= (9) 


is the effective optical depth in the integrated radiation, while ¢ is the optical depth at 
the wave length in question. Radiative equlibrium theory requires that 


R,=bh, (10) 


where 1 — R, is the central intensity of a strong line. Now 1/(1 — R,) is the ratio of the 
flux of the total radiation of the star in the neighborhood of the line to the flux in the cen- 
ter of the strong absorption line, i.e., /,//,o. Burkhardt® has tabulated this quantity as 
a function of the parameters, a = hc/ART, and k/ky. At the line center, k,(a), the line 
absorption coefficient, is very large; and the emergent radiation corresponds to the 
boundary temperature, 7. Given R, and some estimate of the effective temperature, we 
can determine k/k, empirically, provided that our assumptions about the increase of 
temperature with optical depth and the constancy of k/k, with ¢ are fulfilled. In this 
way Unsdld obtained an empirical estimate of k/k, of the order of 2.5 at \ 4000. Since 
R. = to = 0.45 at \ 4000 in 7 Scorpii, the corresponding optical depth in the integrated 
radiation is 1.12. Unsdld then calculated & and k from the theoretical atomic absorption 
coefficients and found that the &/k, ratio should be less than 0.5 for \ 4000. From the 
mass of the material above the photosphere and the computed mean absorption coeff- 
cient, he found the mean optical depth to be tp) = 0.055. The twenty-fold discrepancy be- 
tween the two results was explained by supposing that the large absorption beyond the 
Lyman and helium series acted like the absorption lines in the solar spectrum to retard 
the outward flow of radiation and thereby to modify the temperature distribution in the 
underlying layers. That is, instead of the temperature’s being given by the usual expres- 
sion, 74 = T§(1 + 31/2), the rise in temperature with optical depth is much steeper. 
With this blanketing-effect correction and a mean optical depth of 1.00, Unséld was able 
to reconcile the observed features of the continuous spectrum. 

In a recent paper!’ Chandrasekhar has pointed out that there is no justification for the 
use of the Rosseland mean coefficient of absorption in stellar-atmosphere problems and 
that the proper procedure is to calculate a mean absorption coefficient according to 


(1) 


b= ft. (11) 


° Zs. f. Ap. 13, 56, 1936. 
10 4p. J., 101, 328, 1945. 
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The quantity (F{/F) has been tabulated. Furthermore, if ,/k is independent of the 
optical depth, the temperature distribution on the gray-body assumption continues to 
be valid in the second approximation. 

I have calculated the continuous absorption coefficient as a function of frequency for 
an atmosphere with an electron pressure of 630 dynes (log P. = 2.80) and an effective 
temperature of 29,600° K. Contributions to the continuous absorption coefficient arise 
from photoelectric absorption by hydrogen and by helium and from electron scattering. 
The hydrogenic contribution has been evaluated in the manner indicated by Unsdld, 
the helium absorption is taken from the work of Greenstein for the excited levels" and 
from Vinti’s calculations from the ground state.” The electron scattering is computed on 
the assumption that each atom in the atmosphere supplies one electron. The resultant 
ky, is very similar to that depicted by Unsdéld; slight differences arise from the lower 
helium abundance which I have employed and from the slightly different temperature 
and electron pressure. The adopted effective temperature 7, of 29,600° K corresponds to 
an optical depth of about 0.7; hence the weighting function (F{!/F) chosen is that tabu- 
lated by Chandrasekhar for ¢ = 0.7. Graphical integration yielded a value of k = 3.0 
per gram of the stellar material. The mean absorption coefficient is relatively insensitive 
to optical depth; a calculation carried out for ¢ = 0.10 with the same & gives k = 2.4, 
but I have retained k = 3.0 for the present calculations. The variation of k,/k in the 
observable spectral range is shown in the accompanying table. The straight mean value 


ky ky/k ky ky/k 
3.83 1.30 3650 1.22 0.42 


of R. is about 0.45. If we take the corresponding value of F,/F,,o = 1.82 and use T = 
29,000° K, we find an empirical k/k, of about 3 rather than the theoretical value of 2 at 
\ 4000 from the above table. Thus there still exists a discrepancy of the sign indicated 
by Unsdéld, but the amount is much reduced. 

The mean optical depth, éo, is given by 2 = NHMok or NH4, i.e., the amount of ma- 
terial above the photosphere multiplied by the absorption coefficient. The number of 
atoms above the photosphere of 10 Lacertae is 3.7 X 10”, and the absorption coefficient, 
a, per atom is 6.17 K 10~*4; hence fy) = 0.228 .We must correct, however, for the fact that 
we did not take the negative absorptions into account in computing NH. The mean 
value of (1 — e~*¢/*7e) for X 4000 is 0.70; hence we must correct HN to 5.3 X 10”, and 
we find a mean optical depth fp of 0.33. The theoretical /) computed from equations (7), 
(8), and (9) with k/k, = 2 for \ 4000 gives fo = 0.5 and /p = 1.0. Thus there is a dis- 
crepancy of a factor of 3 between the optical depths computed by the two methods and 
the deviation remains in the direction indicated by Unsdld. 

The present calculations, although approximate, do show that, when the mean absorp- 
tion coefficient is calculated in the fashion suggested by Chandrasekhar, the necessary 
correction for the blanketing effect is largely eliminated. The discrepancies remaining 
seem to be of the order of magnitude of those frequently encountered in astrophysical 
work; and, until detailed numerical integrations for model stellar atmospheres are car- 
ried out, we cannot decide whether these discordances arise from errors in our assump- 
tions (e.g., the constancy of k,/k with optical depth) or from the neglect of the blanket- 
ing effect. A complete theory of radiative transfer ought to be able to handle large varia- 


4p, J., 95, 299, 1942. 
2 Phys, Rev., 42, 632, 1932. 
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tions of k, with ) (such as are encountered in the hotter stars) without invoking addi- 
tional corrections for a blanketing effect, as least when the atmosphere is in local 
thermodynamic equilibrium. 

It is important to establish the relation between the ultraviolet energy distribution 
and the energy distribution in the observable range, in connection with the calculations 
of the temperatures of the planetary nuclei by the Zanstra-Menzel method and for the 
determination of bolometric corrections. 

All atmospheric models that assume uniformly stratified layers provide no mechan- 
ism for the release of extra amounts of energy at the stellar surfaces, either through 
“windows” in the far ultraviolet part of the spectrum or by high-speed particles. The 
model we have assumed requires that the atmosphere become very opaque in the far 
ultraviolet because of the ionization of He1 and, finally, He 1 atoms from the ground 
level. 


In conclusion I should like to express my gratitude to Drs. W. S. Adams and P. W. 
Merrill for permission to use their plates of 10 Lacertae. Dr. Adams and later Dr. 
Bowen kindly gave me the opportunity to work at the Mount Wilson Observatory on 
this problem. Thanks are due to Henry Norris Russell, Jesse Greenstein, Leo Goldberg, 
D. H. Menzel, S. Chandrasekhar, W. Baade, Olin Wilson, and Paul W. Merrill, who read 
the manuscript and made numerous valuable suggestions. 
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THE LIGHT-CURVES OF R AQUARII 


CECILIA PAYNE-GAPOSCHKIN AND CONSTANCE Boyp 
Harvard College Observatory 
Received July 24, 1946 


ABSTRACT 


The visual and photographic light-curves of R Aquarii are presented and analyzed by separating the 
“red component” from the “blue component” photometrically. The relative behavior of the two ‘‘com- 
ponents” is shown to be closely related to the spectroscopic phenomena observed at the Mount Wilson, 
Yerkes, and Harvard observatories. 


The complex variable star, R Aquarii, consists of a long-period variable (P = 383 
days), an extensive and peculiar bright-line nebula, and what appears to be a variable 
stellar component capable of producing high excitation. The purpose of the present paper 
is to analyze the light-curves of R Aquarii in order to examine the bearing of these data 
on the physical nature of the system. 

The variations of R Aquarii have been studied on plates of the Harvard collection. 
One of us (Miss Boyd) measured plates made before JD 27000; most of the plates were 
measured twice. In addition, Miss Boyd determined the magnitudes of comparison 
stars, selected from the visual sequence, on the International scale; these magnitudes are 
basic to the present investigation. Some years later R Aquarii was placed on the system- 


TABLE 1 
COMPARISON STARS FOR R AQUARII 
Designation| DM No. HD IPg _— Designation | DM No. HD IPg ee 
—16°6345. | | —16°6370 | 223315 | 9.31 G5 
45 6507 | 322850: | 7.52) 15 6345 | 222861 | 9.90 F5 
16 6373 | 223428 7.63 KO 10.49 F5 
15 6491 | 223006.| 8.76.) p.....:... 16: 11.35 F8 
—16 6363 | 223306 | 8.72 —16 6355 |........ 


atic program for the study of the brighter variable stars, carried out by C. Payne- 
Gaposchkin and S. Gaposchkin with the aid of a grant from the Milton Fund of Harvard 
University. The measures were made by S. Gaposchkin, and each plate was measured 
once. This series of observations overlaps the period covered by Miss Boyd and extends 
about 3500 days beyond; however, it employs only the small-scale patrol plates, whereas 
Miss Boyd, by studying all possible early plates, achieved more complete coverage during 
the early years. All photographic observations have been reduced with Miss Boyd’s 
magnitudes. They are not published in detail but are available to scientific workers at 
the Harvard Observatory. 

In addition to the photographic record, there exists a rather complete and continuous 
series of visual observations, made by the members of the A.A.V.S.O. and others. 
Mr. Leon Campbell has generously made available these data, combined into 10-day 
means. 

The light-curves are shown in Figure 1. Visual observations are shown by circles, 
photographic observations by dots: small dots denote weight 1; medium dots, weight 
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Fic. 1.—Visual (circles) and photographic (dots) light-curves of R Aquarii 
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2-9; large dots, weight 10 or more. Each estimate contributes weight 1 to the mean, 
except doubtful estimates, which contributed weight 3. 

The well-known disturbance of the light-curve between JD 24000 and JD 28000, as- 
sociated with remarkable spectroscopic phenomena, is shown by both visual and photo- 
graphic curves. The star was evidently bluer than usual during the disturbed interval. 

A complex light-curve may be analyzed on the following simple assumptions: (1) 
Two sources of light are present. (2) These sources, though variable in brightness, are of 
constant color. (3) In the present case it is assumed that the components resemble a long- 
period variable (color index 1.3 mag.) and a blue star (color index 0.0 mag.). The appar- 
ent color index of an unresolved pair of stars, each of known and constant color, is a 
function of their difference of magnitude only. 


TABLE 2 
PHOTOGRAPHIC LIGHT-CURVES OF RED AND BLUE COMPONENTS 


MAXIMA oF In- MIntMa OF IN- Maxima oF In- oF In- 
TEGRATED CURVE | TEGRATED CURVE TEGRATED CURVE | TEGRATED CURVE 
JD JD 

Red Blue Red Blue Red Blue Red Blue 

Comp. | Comp. | Comp. | Comp. Comp Comp Comp. | Comp. 


The light-curves of R Aquarii prove to be difficult to analyze, because the average 
color index has not remained constant, even during apparently undisturbed intervals. 
From JD 21000 to JD 23000 the star was bluer throughout its range than it was either 
before JD 19000 or after JD 28000, when its average color index was 1.3 mag. Possibly 
the abnormal blueness in the 2500 days that preceded the appearance of the ‘‘compan- 
ion” spectrum is real and physically significant. Peculiarities of the light-curve at the 
minima near JD 21400 and JD 21800 suggest that the blue component was even then 
visible; but no colors have been determined at those dates. 

After a number of trials it was decided to obtain a color index for the red component 
from the colors before JD 19000 and after JD 28000, these being the maximum observed 
values. The analysis was carried out from JD 23000 on, for individual maxima and 
minima, with the results shown in Table 2 and Figure 2. 

It is unnecessary to discuss the spectra in the disturbed interval concurrently with the 
light-curve, as a very complete spectroscopic discussion through that epoch has been 
made by Merrill.! In Figure 2 of his paper, Merrill analyzes the light-curve into two 


1 Ap. J., 81, 


312, 1934. 
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Fic. 2.—Brightness of red and blue components at maxima and minima of long-period variable. 


Circles and light lines, red component; dots and heavy lines, blue component. Spectroscopic changes are 


indicated in the figure; the references to the “Ofp” component are from Merrill’s paper; w = weak; 
Ss = strong. 
TABLE 3* 
SPECTROSCOPIC OBSERVATIONS OF R AQUARII 
JD or Year Companion Nebular Lines 
Type Ofp; decreasing in- | Outstanding 
tensity 
1920, Fury, Weak Moderate intensity 
9926, Strong; type Bep; Fem | Weak 
predominant 
Weak; type Ofp Moderate intensity 
Strong; same type as | Very weak 
Dec., 1926 
Weak; [Fem] and [Sm] | Very weak 
predominant 
20520-6190... Y1 [O 11] strong; [O 111] 4363 weak; intensities of Ni, No, 
show effect of occultation by T70; evidence of strati- 
fication 
y2 Return to [Fe m1] stage; [O 11] 


Aa Merrill, The S pectra of Long-Period Variable Stars, p. 84, 1940. 
V1: P. Swings and O. Struve, ibid., 91, 616, 1940. 
Y2: P. Swings and O. Struve, ibid., 95, 159, 1942. 


@ 


| 
| 
rs 
| 
a 
| 
i 
i 
t 
es 
| 
| 
t 
} 
| 
| 
| 
| 
| 
‘te 


LIGHT-CURVES 361 


components; and the close parallel between his results and ours is an indication that our 
analysis has been carried out on acceptable lines. 

Merrill’s analysis goes as far as JD 27700, and our curves extend 3000 days further. 
Swings and Struve have added to the spectroscopic information since the ‘‘companion”’ 
grew weaker. A summary of the more conspicuous spectroscopic changes is given in 
Table 3. 

The light-curves show that both supposed components of R Aquarii are indeed vari- 
able; the large range of the long-period variable is well known, and the blue star also 
seems to vary, by as much as 5 mag. A striking and significant fact is that the range of 


TABLE 4 


CoMPONENT | RED COMPONENT CoMPoNENT | RED COMPONENT 


JD JD 
Average | Average | Average | Average | Average | Average 
Pg. | Maxi- | Mini- eet Pg. | Maxi- | Mini- ace 
Mag. mum mum | Mag. mum mum 
16000-18600...| 11.82 | 7.86 | 12.20| 4.34) 25300-28000...) 9.0 9.81 | 12.05 | 2.24 
21000-23500. ..)....... 8.64 | 12.37 | 3.73 | 28000-31000...| 11.81 | 8.73 | 43: 28:1 4.35: 
| | | 


the long-period variable is suppressed while the blue component is bright. The average 
behavior of the components is summarized in Table 4. 

From Table 4 it may be seen that the early and late ranges are sensibly equal; that in 
the 2500 days before the disturbed period, both the photographic magnitude of the long- 
period light-curve at maximum and the range were below normal; and that during the 
disturbed interval the photographic maximum of the long-period light-curve was de- 
pressed by 2 mag. and its range was reduced to 2.25 mag. Since JD 28000 the maximum 
of the long-period curve has increased progressively. 

The present paper is concerned primarily with the analysis of the light-curves; the 
system is to be discussed in conjunction with other complex systems, such as Z Androm- 
edae and T Coronae Borealis. Apparently, R Aquarii presents a combination of the 
behavior of a long-period variable and a star subject to recurrent outbursts and allied to 
the recurrent novae. Whether the two sets of phenomena proceed from a single or a 
multiple source must be the subject of later discussion. 
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THE LIGHT-CURVES OF Z ANDROMEDAE AND AX PERSEI 


CECILIA PAYNE-GAPOSCHKIN 


Harvard College Observatory 
Received July 24, 1946 


ABSTRACT 


The visual and photographic light-curves are compared and analyzed for Z Andromedae. The behavior 
of the “‘red component” and the “blue component” is analyzed photometrically and is shown to be closely 
related to the remarkable spectroscopic variations of the star. The photographic light-curve of AX Persei 
is shown to bear a strong resemblance to that of Z Andromedae, and the similarity of the spectroscopic 
phenomena in the two systems suggests that they have a close physical resemblance. 


The two stars discussed in the present paper are excellent and similar examples of 
complex variable stars which appear to consist of a red component and a blue star 
that is a close relative of the recurrent novae. Both have been studied in the systematic 
program that was financed by a grant from the Milton Fund of Harvard University. 

Miss Virginia Brenton made 1102 observations of Z Andromedae; the star has already 
been discussed in general terms by Mrs. Greenstein! and by Dr. Richard Prager, who 


TABLE 1 
COMPARISON STARS FOR Z ANDROMEDAE 


BD | HD | Spectrum | I Pg BD | HD | Spectrum| I Pg 
4474177... | 221188 | AO | 7.75 || +48°4094......... | A 10.70 

| | | 
TABLE 2 
COMPARISON STARS FOR AX PERSEI 
Designation BD I Pg | Designation BD | I Pg 


derived a photographic light-curve.2 Mr. Masaaki Huruhata made estimates of 
AX Persei, and Miss Frances Wright made 229 more measures on recent plates. An 
earlier discussion of part of the Harvard material was made by Lindsay.* The photo- 
graphic magnitudes for comparison stars have been redetermined by Miss Frances 
Wright. They are given in Tables 1 and 2. 

The light-curves of Z Andromedae are given in Figure 1, that of AX Persei in Figure 2. 
In Figure 1 are given visual magnitudes, as published by Mr. Leon Campbell from the 


1 Harvard Bull., No. 906, 1937. 
2 Pop. Astr., 49, 445, 1941. 3 Harvard Bull., No. 888, 1932. 
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Fic. 1.—Photographic and visual light-curves of Z Andromedae. Harvard spectra are marked by 
symbol X. 
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observations made by the A.A.V.S.O., combined in 50-day means. The observations 
are not published here, but are available to scientific workers at the Harvard Observa- 


tory. 


It is evident from the figures that both stars display a definite periodicity. Maxima 
of Z Andromedae are tabulated in Table 3; those of AX Persei in Table 4. The lengths 
of cycles for the two stars are compared in Table 5. The similarity is striking. 

The visual and photographic light-curves for Z Andromedae permit the separation of 
the light-curves of the two stars by the method described in an earlier paper.‘ The results 
are given in Figure 3. Unfortunately, no visual light-curve is available for AX Persei. 


x | | 

10+ 4 


Y2 Y3 Y4 Y6 


Fic. 2.—Photographic light-curve of AX Persei. Harvard spectra are marked by symbol X 


TABLE 3 
MAXIMA AND MINIMA OF Z ANDROMEDAE 
JD Maxima | JD Minima ||} JD Maxima | JD Minima |} JD Maxima | JD Minima | JD Maxima | JD Minima 
15000 15350 19600: 20060 23440 23820 27840 28180 
15740 16110 20460 20800 24200 24550 28600: 28830 
16410 16810 20980 21270 24900 25400 29120 29370 
17200 17540 21600 22100 25630 26000 29540 29970 
18100 18400 22240 22670 26370 26700 
18750 19200: || 22820 23380 27200 27460 
TABLE 4 
MAXIMA AND MINIMA OF AX PERSEI 
JD Maxima | JD Minima || JD Maxima | JD Minima || JD Maxima | JD Minima |} JD Maxima | JD Minima 
4 C. Payne-Gaposchkin and Constance Boyd, Ap. J., 104, 357, 1946. 
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The spectral variations of Z Andromedae and AX Persei are no less remarkable than 
their light-curves, and a concurrent study of the two types of data is necessary for an 
understanding of the systems. 

Z Andromedae has been studied, notably by H. H. Plaskett, by P. W. Merrill, and at 
the Yerkes Observatory. A summary of the salient features of the spectra is given in 
Table 6. The times at which the various spectra were obtained are indicated in Figures 
1 and 3, by the symbols given (in parentheses) in the bibliography that follows Table 6. 


TABLE 5 
DETERMINATION OF LENGTH OF CYCLE 
Z ANDROMEDAE AX PERsEI 
Length of Cycle No. of Length of Cycle No. of 
(Days) Cycles (Days) Cycles 
Median cycle......... 714+82 (p.e.) 42 685+79 (p.e.) 28 
JD 15000-29540 | ...... JD 11900-16600 |........ 
Consecutive cycles. . 692 (mean) 21 671 (mean) 7 
WHOS. JD 15350-29970 |........ JD 25150-2580 
Consecutive cycles. . 696 (mean) 21 670 (mean) 4 
Average mean cy- 
24000 25000 26000 27000 
T qT | T T qT T | T T | T T T | i i 
\ 
28000 29000 30000 31000 
or T T T | T T T T T 
e 
| \ 
2 “J 
j 
4b / \ 


Fic. 3.—Resolution{of light-curve of Z Andromedae into red (circles) and blue (dots) components 
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SUMMARY OF SPECTRA OF Z ANDROMEDAE 


TABLE 6 


JD or Year Observer* Description 
23642-23678(10)..| H.H. Plaskett | Emission spectrum representing a combination of (1) a stellar part 
with lines due principally to H, He1, Mg ut and Feu, and (2) a 
nebular part with lines due to He 1, N and [nebulium]”’ 


24761-24873(4)... 
25099-25105(2). . . 


23662, 23674. .... 
1932-1933....... 


29523-29603... .. 


29867-29995... .. 


30000-30215..... 


30360-30400... .. 
30556-30664... .. 


30664-30725..... 


P. W. Merrill 
F. S. Hogg, I 
F. S. Hogg, II 


P. Wellmann 
Yerkes I 


P. W. Merrill 
P. W. Merrill 
Yerkes Ia 


P. W. Merrill 
P. W. Merrill 
Yerkes II 


Yerkes IIT 


Yerkes IV 
Yerkes V 


Yerkes VI 


Yerkes VlIa 
Yerkes VII 


These two spectra were used for a spectrophotometric study, which 
led to a color temperature of 5200°+900° K 

Spectrum as described by Plaskett 

7i0 bands detected on Plaskett’s spectra of these dates 

“Twelve spectrograms obtained during the last two years show 
that the 770 bands strengthened as the star faded. The emission 
lines have weakened relative to the continuous spectrum... .” 

Heu, N ut present, with weak violet absorption; 770 still present 

Cygni spectrum with Plaskett’s spectrum missing; 770, 
He, N 1m, O11, have disappeared; emissions of H, Feu, Ti 
remain with violet absorption borders; forbidden lines weak; dur- 
ing the interval, absorption increases, emission decreases slightly 

Shell-type spectrum; 4686 absent 

4686 barely visible 

Spectrum like Plaskett’s; also a P Cygni spectrum ascribed to ex- 
panding shell; 4363 relatively strong; some evidence of late com- 
ponent 

Shell spectrum disappearing 

Only traces of shell spectrum remain 

Increased intensity ratio, nebular/auroral lines of [O m1]; P Cygni 
absorption disappears; Balmer continuum appears in emission: 


spectrum has post-nova character 
Feit O 111 fluorescent 


Increased intensity ratios: [Vv] 
Het Feu Siu, Mgu 
[Nem] ’ Het, mt)’ Continuum 

No variation in structure of nuclear features, but weaker; decrease 


of [O 111], [Ne 111], [Ne v];slow decrease of 111] since 1940; 


strong Balmer continuum in emission; ian hall still near photo- 
spheric surface 

[Ne 11] increased; now as in August, 1940 

Strong [Fe vu, also [Fev]; peak excitation on JD 30500; hardly 
any absorption, but evidence of late spectrum in visual region 

Slow decline in excitation; [Fe vir] and [Fe v] still present (30725) 
but weaker than on 30541 or even than on 30664 

Increase of excitation 

Increase of [Fe vit] and [Nev]. Fert very weak; Fe decreasing; 
[Fe 1] increasing 


* Authors and symbols: 


H. H. Plaskett, Pub. Dom. Ap. Obs., 


Vol. 4, No. 10, 1928 (P). 


P. W. Merrill, Ap. J., 99, 15, 1944 (M). 


F. S. Hogg, I, Pub. A 


.S.P., 44, 328, 1932 (H). 


F. S. Hogg, II, Pub. A.A.S., 8, 14, 1934 (H). 

P. Wellmann, Beob. zirk., 21, 105, 1939; V.J.S., 75, 53, 1940 (W). 
Yerkes I: P. Swings and O. Struve, Ap. J., 91, 601, 1940 (Y1). 
Yerkes Ia: P. Swings and O. Struve, Jbid., 93, 356, 1941 (Y1a). 
Yerkes II: P. Swings and O. Struve, Jbid., 94, 296, 1941 (Y2). 
Yerkes III: P. Swings and O. Struve, Jbid., 95, 152, 1942 (Y3). 
Yerkes IV: P. Swings and O. Struve, Ibid., 96, 257, 1942 (Y4). 
Yerkes V: P. Swings and O. Struve, Jbid., 97, 194, 1943 (Y5). 
Yerkes VI: P. Swings and O. Struve, Jbid., 98, 93, — (Y6). 
Yerkes VIa: O. Struve, Ap. J., 99, 209, 1944 (Y6a) 

Yerkes VII: P. Swings and O. ‘Struve, Ap. J. 230, 1945 (Y7). 
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Nineteen short-dispersion spectra of Z Andromedae, taken at Harvard with the 8-inch 
Draper telescope and one or two objective prisms provide some additional information, 
which is given in Table 7. Successive columns give the (hitherto unpublished) comments 
of Miss Cannon, Mrs. Mayall, and the writer. The two last spectra in the list were 
taken by Mrs. Mayall, who has kindly made them available. 

The Harvard spectra are somewhat scattered in time, and the first six fall in an inter- 
val in which the course of the light-curve is uncertain. The remainder are collected in 
Table 8, which relates position on the light-curve to general spectroscopic behavior. 

The information pieced together for Z Andromedae f:om light-curves and spectra 
may be summarized as follows: 

Both ‘‘components” vary together with the 690-day cycle. The blue “component” 
has an observed range of about 6 mag.; the red “‘component”’ varies by 2} mag. at most. 

At the highest maxima of the blue component, the spectrum is that of a minor nova 
outburst, which is followed, as the star declines, by the typical spectral development of 
a slow nova. 


TABLE 7 
HARVARD SPECTRA OF Z ANDROMEDAE 


JD Miss Cannon Mrs. Mayall Mrs. Payne-Gaposchkin 
$3014... 2s Poor; bright lines seen 4686=3; Hy=1 4686; Hy weaker 
16). Spectrum seen, no bright | Very poor, but bright lines | Too poor 
lines distinctly seen in red 
or bands 
HB=2, 4686=3, Hy=3, | HB, Hy, Hé, decreasing in- 
Hé=1 tensity, 4686 
HB, 4686; Hy, Hé, He, | HB=4, 4686=5, Hy=5, | HB, Hy, Hi, He 
Hy =4686; H fainter Hé=1, He=1 
56142... ... Bright lines seen; too near ae 4686=4, Hy=4, | HB, 4686, Hy, H4; very poor 
edge Hé= 
16462)... 3s: HB, Hy, Hé, 4686 distinctly | HB=4, 4686=5, Hy=5, | HB, 4686, Hy, Hé; weak con- 
seen; He ? Hb=3 tinuum; good 
16467 HB, 4686, Hy, He seen; | H8=3, 4686=4, Hy=5, | HB, 4686, Hy, Hd; Hy is very 
all bright, Hy perhaps Hé=2 strong; medium continuum 
the strongest 
16485. HB, Hy, H6; all are Hy=3, | HB, 4686, Hy, 4686= Hy 
right 6= 
Too faint, barely seen 4686=1 Continuum just seen; if 
bright lines, they are very 
faint 
bright lines they are very 
faint 
lines 
HB=2, 4686=1; no back- | Two bright lines 
ground 
Very bright; spectrum Continuum very strong, 
burned out, but Hf ap- bright Hg, Hy? 
pears to be bright=1 
He=1 
Hé=1 medium continuum 
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Titanium oxide is recorded during the intervals when the red component is relatively 
bright (JD 23600 on; JD 27000 [strengthened]; JD 29475; JD 30600); it disappears 
during high maxima. 

The color temperature of 5200°, determined by Plaskett, is consistent with the color 
index of +0.8 observed at that time. 

The successive appearance of He 11 4686 as the star declines and of [O 111] 4363, fol- 
lowed by the nebular lines, which increase in intensity relative to the auroral line, all 
seem (in the light of Table 6) to be part of a regularly recurrent pattern. 


TABLE 8 
SUMMARY OF INFORMATION FROM EARLY SPECTRA 


Position on Light-Curve on Description 
{ Faint continuum; H (? 4363) seen 
Faint, just before rise..| 20072 | No continuum; H, 4686 seen on one plate, 4686 alone on another; 
spectrum apparently varying rapidly 
Rising to maximum....| 14874 | Balmer lines bright; 4686 strong; ? 4363 present 
Bright maximum...... {30613} Strong continuum; Balmer lines bright; no other bright lines seen 
Decline from bright 
MARZIMUM.... 2.6.5 «5 29849 | Medium continuum; Balmer lines bright; 4686 strong; no evidence 


of 4363 


ate maximum........ 16467; | Weak continuum; Balmer lines bright; 4686 strong; evidence of 4363 


Decline from intermedi- ote 
16485 


Intermediate minimum. { e142} Balmer lines; 4686; 4363? 


Decline to faint mini- 
eg ee 16787 | Continuum barely seen; bright lines very faint; 4686 seen 


TABLE 9 


JD | References* Remarks 

29500-600.. Yerkes Weak continuum, corresponding to gM3e; very high excitation; 
strong Balmer continuum in emission; bright H, Het, 
He ut, [Fe vit], [Ne v], [Fe x]? 

29660-30000. ......... Yerkes II [Fe [Fe v1] disappeared; Ni, N2 decreased relative to 4363; 
[Ne 11]/H, 3888/[Ne 11] 

[Ne ut] weaker than 3888; 4472 strengthened relative to 4640; 
Ni, N2 decreased; no [Fe vit] 

Yerkes III | [Fe vu], [Nev] again strong; decreased Het, Nim, C1, 
{O m|, [Ne 11]; excitation has increased since JD 30000 

Yerkes IV [Ou] fairly strong; [Fe vir] weaker; auroral/nebular [O m1] 
decreased; [Fe vj, [Fe v1] strong 

Yerkes VI | Again increasing in excitation; [Fe vu] as on JD 30300 [Ne m1] 
weaker; [Fe v] still present, weakened 


* References: 
Yerkes I: P. Swings and O. Struve, Ap. J., 91, 607, 1940. 
Yerkes II: P. Swings and O. Struve, ibid., 94, 299, 1941. 
Yerkes III: P. Swings and O. Struve, ibid., 95, 152, 1942. 
Yerkes IV: P. Swings and O. Struve, ibid., 96, 254, 1942. 
Yerkes VI: P. Swings and O. Struve, ibid., 98, 93, 1940. 
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Spectroscopic material for AX Persei is summarized in Tables 9 and 10. The informa- 
tion complements that for Z Andromedae, which, for the most part, covers bright, 
disturbed intervals. For AX Persei the information begins during the faint, quiescent 
period and covers one outburst of intermediate intensity. Only two of the Harvard 
plates, taken in 1914 during a high maximum, show appreciable blue continuum. 


TABLE 10 
HARVARD SPECTRA OF AX PERSEI 
JD Mrs. Mayall Mrs. Gaposchkin 
ver HB=2, Hy=2, Hi=1 Medium continuum, runs far to violet; 
HB=2?, Hy=1?, H6? 
ye HB=1? Medium continuum; bright lines? 
Hp=3, Hy=1, HB=2, Hy=1, Hi=1 
Hp=1, Hy? 
HB=2, Hy=1, Hg=1, Hy=}3 
Hp=3, Hy=1 HB=2, Hy=1 
29522........| HB=3, Hy=1, Hg=2, Hy=1, Hi=} 
} Hp=4, 4686=1, Hy=3, Hi=1 HB=2, 4686=1, Hy=2, Hi=1 
29524........| HB=5, Hy=3, Hi=1 HBp=2, Hy=1, Hi=1 
HBp=4, 4686=1, Hy=3, H5=1 HB=3, 4686=1, Hy=2, Hi=1 
29577.........| HB=5, 4686=1, Hy=3, Hi=2, He=1 HB=4, 4686?, Hy=2, Hi=1, He=} 
HB=5, 4686=1, Hy=3, Hi=2, He=1 HB=3, 4686=}3, Hy=2, Hi=2 
HBp=2, Hy=1 Hg=1, Hy=} 
ye HB=5, no 4686, Hy=4, Hi=2, He=1 HB=4, no 4686, Hy=3, Hi=2, He=1 
HBp=2, Hy=1 Hp=1, Hy=}3 
29931........| HB=2, 4686=2, 4640=}, Hy=1, Hi=} 
29934........| HB=4,4686=3,4640=1, Hy=1,Hi=1| HB=2, 4686=2, 4640?, Hy=1 
|? HB=3, 4686=3, Hy=1, Hb=1 - HB=2, 4686=2, 4640?, Hy=1 


During the low maximum at 29500, the Harvard plates show the bright Balmer lines 
strong and 4686 weak or missing. During the rise to the intermediate maximum at 
29930, the Balmer lines and 4686 were strong, the N 111 lines at 4640 were present but 
weak. This behavior is similar to that described in Table 8 for Z Andromedae, but the 
excitation is higher throughout. 
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THE ECLIPSING SYSTEM UV LEONIS 


SERGEI GAPOSCHKIN 
Harvard College Observatory 
Received August 3, 1946 


ABSTRACT 


The eclipsing system UV Leonis has been studied on 18 spectrograms made at the McDonald Observa- 
tory and on 562 photographic Harvard patrol plates. The spectrum outside eclipse shows two com- 
ponents of spectral class GO and G2. Orbital velocities of both components were obtained. The spectro- 
scopic elements are: K; = 165 km/sec; K2 = 180 km/sec; (a; + a2) sin 7 = 2.85 K 10° km = 4.10 ©; 
(m, + m2) sin?i = 2.56 ©. The photometric elements are: Min.= 24256724722 + 0460008582E; Max. = 
848; min.; = 9™28; min., = 9™22 (not rectified); k = 0.96 (a = 1.0);7 = 90°; a2 = 0.41; bp = 0.39; a, = 
0.43; b; = 0.41; s = 0.10; y = 1.06; LZ, = .049; and LZ; = 0.51. The combination of spectroscopic and 
photometric elements makes the system UV Leonis a well-determined one. The absolute dimensions are: 
R.-= 1.12 ©; Ri = 1.17 ©; m, = 1.34 ©; m, = 1.22 ©; Spe = GO; Sp: = G2; Mvis, = 4™74; Mvis; = 
4™91, the two latter being computed from R’s and adopted temperatures of 5750° and 5510° C, respec- 
tively. A striking difference between the less massive component of this system and that of RT An- 
dromedae (of the same period) is indicated. 


I, SPECTROSCOPIC OBSERVATIONS 


At the McDonald Observatory, in December, 1944, I obtained 18 spectrograms with 
the 82-inch mirror for the eclipsing system UV Leonis. All the spectrograms appeared to 
be of good quality. The spectrum on 10 of them shows double lines. Both spectra are of 
good quality and of G type; the classification based on the metallic lines is adopted as 
GO and G2. Nine metallic lines, chiefly of iron, together with Hy and H6, were measured. 
Each plate was measured once by the writer; the reduction was performed by Mrs. C. 
Payne-Gaposchkin. The distribution of the spectra over the phases is favorable. Table 1 
gives the spectroscopic observations with additional information (see also Fig. 1). 

The plotting of the radial velocities against the phase shows a satisfactory agreement 
between photometric and spectroscopic observations. Though the number of spectro- 
grams is not large, they seem to indicate a small eccentricity of the orbit. The method of 
Lehmann-Filhés has been used, and the results are given in Table 2. 


II. PHOTOMETRIC INVESTIGATION 


The 562 photographic plates used have been measured twice at different times. The 
field is not rich in stars, and the number of comparison stars is limited. Three were 
selected, as given in Table 3; their magnitudes were obtained by the author from an 
unpublished sequence at Harvard and by comparison with the neighboring Harvard 
Standard Region C6. The scale and color were also investigated on Harvard yellow 
plates and visually at the McDonald Observatory. 

The period has been discussed by means of twenty-five times of primary and thirteen 
times of secondary minima over an interval of forty years. A new period was found; the 
elements are Min.=2425672.722 + 0460008582. Since the period is short and the ex- 
posure time of the plates is long, the observed brightness on a given plate may be in error. 
The error may be large not only where the change of brightness is rapid but also where 
the light-curve has greater curvature, e.g., at second and third contact and on the 
shoulders of the light-curve. Since there is no method of evaluating rigorously the neces- 
sary correction to magnitude, a simple graphical procedure, involving the determination 
of the “‘center of gravity” of the exposure-time effect, has been used. First, a light-curve 
was obtained from the original observations, and the shift of the center of gravity of a 
number of points on this curve was determined for different exposure times. The uncor- 
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TABLE 1 
VELOCITIES OF UV LEONIS IN KM/SEC 


Exposure} Velocity 
JD 7) Time All 
0 (Minutes)| (All but H) 

0°531 50 —221 (8) —215 (11) 
+ 87 (7) 82 (8) 

142 45 + 32 (6) + 32 (6) 
—172 -(8) — 92 (11) 

i es .644 40 — 55 (6) — 56 (8) 
+ 30 (2) 30 (2) 

| ae 347 37 —135 (7) —132 (9) 
+ 79 (5) + 70 (7) 

| 015 38 +132 (6) +123 (8) 
—178 (8) — 167 (10) 

218 40 — 46 (9) — 50 (11) 
697 35 — 24 (10) — 27 (12) 
307 41 —154 (7 —145 (9) 
+ 86 (6) + 84 (7 

374 40 —188 (8) — 200 (10) 
+104 (8) + 98 (10) 

.553 40 —189 (4) —165 (7) 
+ 92 (3) + 81 (6) 

.066 40 + 81 (6) + 81( 6) 
—183 (4) —165 (7) 

.196 40 — 58 (8) — 48 (10) 
.689 40 — 48 (8) — 46 (10) 
07 7 .812 35 — 60 (7) — 54 (9) 
DOG .347 40 —134 (4) —133 (6) 
+ 84 (3) + 69 (5) 

31456-87276. 2 629 45 —120 (6) — 91 (8 
+ 26 (6) + 26 (6) 

683 45 — 47 (7) — 47 (8) 
31456. 0.735 45 — 49 (6) — 46 (10) 
Po I 2 4 5 6 7 8 x) xe) 

100k 4 


K, = 165 


Fic. 1.—Velocity-curves of UV Leonis 


km/sec 


K, = 180 km/sec 
Ratio of masses = 1.09 


TABLE 2 


SPECTROSCOPIC ELEMENTS OF UV LEONIS 


(a; + a2) sini = 2.85 K 10° km = 4.10 
+ mz) sin? i = 2.560 


A 
| | 
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-100 
-200 
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rected light-curve shows a considerable ellipticity, z = 0.20, with the two ranges 0768 
and 062. After correction for exposure-time effect, the ellipticity was reduced, z = 
0.10, and the minima were deeper, 0™80 and 0™74. As an illustration, a part of the cor- 
rection table is given in Table 4; tabulated quantities are hundredths of a magnitude. 

The rectification of the light-curve has been made graphically by plotting 7 against 
cos 26; the deduced value of z is 0.10. The application of the method of least squares to 
the formula / = a+ b cos 6+ ¢ cos? 6 has led to s = 0.0905 + 0.0203; y (coefficient 
of reflection effect) = 0.0236 + 0.0088; and x = 0.9999 + 0.0191. There is a slight indi- 
cation that the minima are separated by not exactly half the period; the phases of the 
minima being 07755 and 0?245, respectively. However, in obtaining the dimensions of 
the system, this small eccentricity has been neglected. The rectified light-curve is given 
in Table 5. 

The maximum has been adopted from observations outside the eclipse as 8™48; the 
ranges then are 0™74 and 0™68. From these values there is very little latitude for a and 


TABLE 3 
COMPARISON STARS FOR UV LEONIS 


Designation BD IPg 
14 2271 (9.1) 9.11 
Po 4 2 3 4 5 6 7 8 9 0 
T T T T T T T T 
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Fic. 2.—Light-curve of UV Leonis. Circles, observed normals; dots, computed uniform curve 


TABLE 4 
CORRECTION FOR EXPOSURE TIME 


EXposuRE 

TIME 

(MvuTes) | | oP10 | oP20 | oP2s | oP3s | oPso | oPos | oP7s | oP90 
4 16 4 


—12 | +34 —8 


\| 


ad 


ge 
mi 
pa 
2 
th. 
sul 
TI 
~ 
0 
| | 
by 
Tal 
iv 
give 
and 
on 
P 
of 
the 
I 
8 
gl 
oF Exc 
aut] 
Lac 
1 
2 
ad 


UV LEONIS 373 


k; the latter must lie between 1.00 and 0.96. In order to use the secondary minimum, to- 
gether with the primary, in investigating the elements, the magnitudes of the secondary 
minimum were reduced to the amplitude of the primary; the numbers are given in 
parentheses in Table 5. The errors of the normal points of Table 5 cannot be greater 
than 0704. 

The light-curve (the ascending and descending branches) is more consistent with as- 
sumption of a total eclipse and the ratio of radii 0.96 than with k = 1.00 and k = 0.96. 
The run of observations is, indeed, such that it favors a still smaller value of &, in viola- 
tion of the limits set by the depth of the secondary minimum. This is well illustrated 


TABLE 5 
RECTIFIED LIGHT-CURVE OF UV LEONIS 

Phase IPg Phase IPg Phase IPg Phase IPg 
0070". 9:22 2273. 8.52 8.97(9.01) || 0.7239..... 8.48 
9.15 8.45 .4869...... 9.16( .22) yk 8.48 
9.06 8.46 9.09( .15) 8.48 
8.61 8.52(. 53) 8.69( .73) 8.62 
$200. 8.52 8.58(.60) 8.51( .54) 8.91 
8.44 8.56(.58) 8.49( .50) |} 0.9874..... 9.21 

8.49 || 0.4401...... 8.74(.77) || 0.6748...... 8.45( .46) 


* Adopting the middle of the primary minimum, ¢ = 07750. 


TABLE 6 
ELEMENTS OF UV LEONIS 
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1.06 = 0.39 
0.49 z = 0.10 
0.51 b/a = 0.95 


by plotting the intensity (J) against sin* yg. The solution of the light-curve is given in 
Table 6, for uniform solution, using S. Scharbe’s method.! The computed light-curve, 
given in the figure, shows a systematic deviation at the bottom of the minimum. 


III, ABSOLUTE DIMENSIONS AND DISCUSSION 


The system of UV Leonis is one of the very few eclipsing variables of short period 
and late spectrum for which adequate data exist to obtain masses and radii. Both com- 
ponents have well-defined orbital velocities, and the minima are so deep that the ratio 
of radii is practically uniquely determined. The combination of the spectroscopic with 
the photometric elements leads to the absolute dimensions given in Table 7. 

It is of interest to compare these values with those for stars of similar character. Table 
8 gives all known G stars, other than the sun, for which masses have been determined. 
Except for those indicated otherwise, all values are taken from the compilation by the 
author.2 There is a striking similarity in masses for UV Leonis and systems like AR 
Lacertae, WZ Ophiuchi, and RT Andromedae, the first two of which are determined 


1 Poulkovo Bull., 94, 31, 1923. 
* Harvard Reprints, No. 201, 1940. 
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Fic. 3.—Relative dimensions of UV Leonis 


TABLE 8 
MASSES OF G STARS 

Star Spectrum Mass Period Star Spectrum Mass Period 
GO 1.500} 0463 G2 1.410} 1490 
GO 1.41 4.16 G2 2.03 3.38 
OS GO 1.34 0.60 RS CVn G8 1.74 4.80 
G1 0.50 0.41 RT And* G8 0.99 0.63 


* C. Payne-Gaposchkin, Ap. J., 103, 291, 1946. 
t S. Gaposchkin, Pub. A.S.P., 56, 199, 1944. 
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with very high accuracy. Even more interesting is a sharp difference between the masses 
of short-period systems, like W Ursae Majoris and ER Orionis, and those of the rest of 
the group. This may be partly due to the fact that the former have periods of less than 
half a day, that their components are close to each other, and that the reflection effect 
on the radial velocities (making the amplitude smaller) is considerable, which will lead 
to a falsification of the masses. 

Another point of interest is the comparison between the system of RT Andromedae 
and the present one. They have similar periods (0.6 days), about the same ellipticity, 
and the primaries are of about the same spectral class; but the secondary spectrum of 
the system RT Andromedae is of definitely later type than that of UV Leonis. There is, 
however, a great difference in photometric elements, exemplified in the different light- 
curves. RT Andromedae represents a combination of a main-sequence star and a sub- 
giant instead of a pair of similar stars, as with UV Leonis. To understand why they are so 
different in one respect and so similar in another involves a knowledge of the origin and 
evolution of double stars. 


The writer wishes to express his appreciation of the hospitality of Dr. Struve during 
his stay at the McDonald Observatory, which made possible the securing of the spectro- 
grams used in the present investigation, and of the generous permission of Dr. Shapley 
to make this visit. 
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THE ECLIPSING SYSTEM RX GEMINORUM 


SERGEI GAPOSCHKIN 
Harvard College Observatory 
Received August 3, 1946 


ABSTRACT 


The system has been studied on 19 spectrograms of the McDonald Observatory and on 1376 photo- 
graphic patrol plates of the Harvard Observatory, supplemented by 239 visual estimates by Enebo. 
Outside the eclipse the spectrum shows but one component, of spectral class A4. Orbital velocities were 
obtained only for this component. At minimum—a constant brightness of short duration—the spectrum 
is practically unchanged. No trace of a secondary spectrum was observed. The spectroscopic elements 
are Vo = +33.5 km/sec; K = 25.5 km/sec; e = 0.295; w = 21°7; f = 0.018 ©; a2sini = 4.07 X 10°km 
= 5.85 ©. The photographic light-curve exhibits a small reflection effect, slight asymmetry of the 
primary minimum (fall is slower than rise), and a constant phase of small duration at the bottom of the 
minimum. The photographic elements for (D) solution are: Max. = 9™20, A; = 2™02, A2 = 0702, k= 0.40, 
r, = 0.188, r. = 0.075, i = 83°9, L, = 0.84, and y = 30.9; the computed spectrum of the invisible compo- 
nent is KO. The visual light-curve is less well determined than the photographic one. Its relative elements 
for (U) solution are: Max. = 9™25, A; = 1™76, k = 0.33, nr, = 0.198, re = 0.075, i = 83°0, L. = 0.82, 
y = 31.5; the color index during the total minimum is 0"21, at maximum, 0™05. The spectroscopic and 
photographic observations indicate that RX Geminorum is a typical system, consisting of a main-se- 
quence star and a subgiant. The approximate dimensions of the components are: masses, 3.1 © and 
0.62 ©; radii, 2.20 © and 5.5 ©; absolute visual] magnitudes, 1.1 and 2.7; and spectra A4 and KOp. 


The observed spectrum at primary minimum is, however, of class A. 


I, SPECTROSCOPIC WORK 


During my stay at McDonald Observatory I obtained 19 spectrograms of RX Gemi- 
norum with the Cassegrain spectrograph of the 82-inch reflector. The spectrum is simple, 
with well-developed hydrogen lines. Some fainter lines at \\ 4172-4178 and Ad 3900- 
3913 indicate that the star is, if not a supergiant, at least not a normal main-sequence 
star. Only one spectrum (A4) was measurable. All measurements were made by the au- 
thor and all reductions by Mrs. Cecilia Payne-Gaposchkin. The radial velocities from 
hydrogen, calcium, and iron lines are given in Table 1. The distribution of the plates is 
favorable, though their number is too small for detailed discussion of such spectroscopic 
phenomena as the rotation effect, which is probably present in the spectrum of this 


system. 
II. SPECTROSCOPIC ORBIT 


Though the plates are not numerous, the general trend of the radial velocity-curve of 
the brighter component (A4) is unmistakable; the orbit is eccentric. Using the Lehmann- 
Filhés method, we obtain the spectroscopic elements for all lines as shown in Table 2. 

I have only one plate on the descending branch of the minimum close to totality. 
There is definite indication that the minima came earlier in 1944-1945 than predicted by 
the ephemeris, so that the plate is taken, if not entirely during the totality, at least dur- 
ing a part of it. A simultaneous visual estimation at that phase was easily made, for there 
is a very close companion star at 11.0 mag. The visual magnitude of RX Geminorum was 
only 0.1 brighter than this companion—that is, 10.9—at JD 30435.959. Our spectrum 
shows at this phase the K line unchanged and sharp. On this basis, therefore, the spectral 
class would be about A5. Wyse, who classified the star at minimum on spectrograms of 
small dispersion, has given it as A3p. But on our spectrum there is a definite indication 
of the presence of \ 4227, which is not visible on any of the other spectra. The hydrogen 
lines are definitely present. The primary presents no difficulties for classification; it is 
about A4 in the Morgan system. The distribution of energy in the spectrum suggests 4 
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TABLE 1 
VELOCITIES OF RX GEMINORUM IN KM/SEC 
JD Phase H Feu Cam 
0°306 +17: (4)| +50: —22: — il: 
pa eee .554 10.1 (5) 6.5 + 8.5 + il 
.631 30.8 (6) +27 + 34 
.720 51.8 (5) 59 +44 + 57 
0, .886* 27.2 (4) 29 +29 + 22 
‘re .040 18 (5) 15 +25 + 10 
im .205 25 (4) + 13 
its . 280 4 (4) + 13 
.358 20.2 (5) 18.5 +21 + 22 
he .446 25.8 (5) 24 +24.5 + 32 
he 14.8 (5) 10 +11 + 32 
0, .616 37.2 (5) 38.5 +22.5 + 44 
0- .694 61 (5) +65 + 68 
its .780 65 (5) 65 +48 +100 
32, .862 66 (5) 63 +66 + 72 
nd .944 50.8 (4) 38.5 +51 + 75 
.596 45.8 (4) 48 +40 + 47 
nd 0.745 +66 +70.5 +44.5 +100: 
p. 
* In minimun, faint. 
‘i TABLE 2 
% SPECTROSCOPIC ELEMENTS OF RX GEMINORUM 
“e K, = 25.5 km/sec Vo = 33.5 km/sec 
f = 0.018 a, sini = 4.07 = 5.860 
e = 0.295 P = 124209* 
w = 21°7 
* From the photometric work. 
is 
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n Fic. 1.—Velocity-curve for RX Geminorum; ordinates are radial velocities (km/sec), abscissze are 


n phases in units of the period. 


‘ 
IS 
a 


378 SERGEI GAPOSCHKIN 


lower temperature, in conformity with the surface brightness. The velocity-curve shows 
a possible rotation effect. No doubt the star deserves more attention spectroscopically, 


III. PHOTOMETRIC OBSERVATIONS 


The photographic observations (1376) will be discussed first. The comparison stars 
were determined with reference to Selected Area 50, with the stars of 8.55, 10.34, 11.99, 
12.00, and 12.08 mag. Eight plates with different exposure time from 60 to 120 minutes 
were used, and the magnitudes are given in Table 3. The star « is the one very close to the 


TABLE 3 
COMPARISON STARS FOR RX GEMINORUM 
Star Scale BD I Pg Vis 
+8 34°1458 (9.1) 8.62 +0.03 8.95 
+4 33°1412 (8.8) 9.00:+ .12 9.26 
— 6 33°1413 (9.3) 9.96 + .02 9.51 
—15 33°1416 (9.3) 10.63 + .02 9.95 


variable. Some other comparison stars of intermediate brightness were used, which, how- 
ever, have no BD numbers. The visual magnitudes were determined in reference to the 
photovisual magnitude given by C. Wirtanen and A. N. Vyssotsky' for their field 35°11, 
together with other sequences. 

I have used over twenty minima to investigate the possible correction of the period, 
over an interval of forty years. The period seems to have been constant until JD 2428000, 
when it shortened slightly (see Fig. 2). 


IV. THE PHOTOGRAPHIC LIGHT-CURVE 


All photographic estimations were combined into fifty-one normal points; the light- 
curve is given in Table 4. The phases go and ¢ are the original (referred to JD 2400000), 


10 20000 


25000 30000 
bd a 


00 


Fic. 2.—Phase of minimum of RX Geminorum. (referred to primary minimum of mean light-curve) 
as a function of Julian Day. 

Fic. 3.—Reflection effect for RX Geminorum; ordinates are photographic magnitudes; abscissae are 
phases. 


and the one reduced to the middle of the primary minimum (¢) = 0.902) as the zero 
point. 

There is a definite indication of a reflection effect (see Fig. 3); the correction for the 
reflection effect has been made with the formula Am = 0™01 cos ¢. The minimum seems 
to be slightly asymmetrical, the descending branch the steeper. 


1Ap. J., 101, 141, 1945. 
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Solutions based on the method of Scharbe? for uniform and total darkening were car- 
ried out. Both represent the observations well within the observed uncertainty. The 
probable (or mean) errors can be inferred from the fact that the probable error of a 
single estimate is between 0709 and 0™12. 


V. VISUAL OBSERVATIONS 


The visual observations of eclipsing binaries have until recently been progressively 
neglected, owing to the efficiency of the photographic observations. But an investigation 
of thick atmospheres requires independent work in several wave lengths. The system of 


TABLE 4 
PHOTOGRAPHIC LIGHT-CURVE FOR RX GEMINORUM 
¥0 IPg #0 n IPg 
38 9.20 .128 30 9.21 . 786 
29 9.20 . 187 42 9.19 .831 
30 9.18 .310 17 9.23 .951 
24 9.18 .332 14 9.25 .959 
32 9.18 . 384 12 9.73 .977 
44 9.18 .406 12 10.50 .981 
54 9.19 .458 8 11.10 .993 
32 9.20 .482 10 11.22 .000 
29 9.21 .535 7 11.07 .009 
32 9.18 . 560 7 10.68 .014 
59 9.18 . 588 7 10.27 .021 
30 9.16 .638 16 9.35 .032 
28 9.18 .659 13 9.21 .043 
41 9.19 . 706 24 9.16 0.086 
24 9.18 0.738 


RX Geminorum consists of an A star and a component of later spectral class. It has a 
well-defined totality, and the period is 12 days; we might expect a thick atmosphere ef- 
fect, though not a large one. It is unfortunate that there is no well-observed visual light- 
curve to be compared with the photographic curve. All that I have are 239 visual obser- 
vations of Enebo, made more than thirty years ago.* The distribution of the observations 
is rather unfavorable, for the observer had the object of determining the period of the 
star rather than a reliable light-curve. Since the magnitude scale of comparison stars used 
by Enebo was not on the international scale, I have determined magnitudes for the re- 
duction of his observations, The reduction has been made with the help of the following 
relations: 


S.G. = 9"°10+ 0.83 (Enebo — 8"5) 


? Poulkovo Bull., 94, 31, 1923. 
§ Ark. f. Math. och Naturvidensk. B, Vol. 33, No. 8, 1913. 
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to m = 10™18 and after m = 10.18 


S.G.=10750+ 1.11 (Enebo— 1018) . 


The visual light-curve of RX Geminorum is given in Table 5. 


TABLE 5 
VISUAL LIGHT-CURVE OF RX GEMINORUM 

n IPg n IPg ¢ I Pg 

.007... 5 10.78 ee 6 9.21 .882.. 6 9.20 
. 5 10.55 .409... 6 9.25 .909... 7 9.30 
5 10.19 6 9.29 .945... 9.32 
5 9.32 .476.. 6 9.23 5 9.26 
5 9.24 7 9.23 .965.. 9.34 
.046.. 5 9.24 7 9.31 .969... 5 9.50 
.079... 10 9.25 6 9.23 5 9.75 
2 9.32 .634.. 6 9.25 10.09 
i ee 9 9.21 .673... 9 9.27 .979. 5 10.16 
9 9.24 8 9.2i .982. 5 10.24 
9 9.28 7 9.26 .986. 5 10.51 
0.270.. 10 9.22 0.806... 8 9.23 SS s 10.71 
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Fic. 4.—Primary minima of RX Geminorum 


(D) curve; circled crosses, computed (U) curves. 


O sun 


Fic. 5.—Relative dimensions for the system of RX Geminorum 


. Large dots represent normals, small dots, computed 
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A uniform solution was carried out for this light-curve, exactly as for the photographic 
light-curve. The results for the several solutions are summarized in Table 6. 


VI. DISCUSSION 


The system RX Geminorum is a normal one in respect to period, range, and photo- 
metric and spectroscopic elements. In order to see the whole interplay of the components 
we must consider the absolute dimensions. By applying the mass-luminosity relation, 
together with the relative elements for (D) solution, we find a mass for the A4 compon- 
ent equal to 3.1©; radius, 2.90; and the visual absolute magnitude, +1.1. We find the 
spectral class of the larger component by using the ratio of the surface brightness (y = 
30.9) and the arbitrarily adopted photographic surface brightness 34 for the brighter 
component. The deduced spectrum is KO. With this spectral class and with the relative 
radius we can find the mass and radius of this component. But similar dimensions can 
be found from the A4 component and relative elements and from the mass function, by 


TABLE 6 
PHOTOGRAPHIC AND VISUAL ELEMENTS OF RX GEMINORUM 
PHOTOGRAPHIC 
(1376 OssERv.) VISUAL 
(239 OpsERv.) 
U 
U D 
Radius-of A ater... 0.065 0.075 0.075 
OF 0.192 0.188 0.198 
Ratio of surface brightness. . .. . 46.1 30.9 30.9 
Spectrum of the smaller star....| A4 A4 A4 
Spectrum of the larger star..... K4* K4* K4* 
Maximum brightness.......... 9™20 9=20 9™25 
Minimum brightness........... 14.22 11.22 11.01 


* From surface brightness. 


adopting a mass of 3.1© for the A4 component. Table 7 gives the intercomparison of the 
results. 

According to Table 7, the system is a combination of two stars, of which one is on the 
main sequence in every respect—spectrum, A4; mass, 3.1© ; radius, 2.5© ; and luminos- 
ity, +1™1; the other is a typical subgiant of spectral class KO with a radius of 5.50 
and luminosity +2™7 but with mass probably less than that of the sun. The large red 
tar totally eclipses the smaller one for 3.4 hours; but the spectrum remains practically 
the same as at maximum, A4. The color index 0™21 at the minimum indicates also that 
the spectral class is near A. This system is not the first one for which such a discrepancy 
\AS observed and KO predicted) has been found, but it is one of the most conspicuous. 
Wyse,‘ who first called attention to the fact that the spectrum during primary minimum 
sof class A, suggested tentatively the presence of a third body. Wyse’s spectrum was 
taken entirely within totality: our spectrum, centered at phase 0.886, is still within the 
partial phase but at a time when 0.878 of the brighter component is covered. If the lumi- 
wosity found in our solution is adopted, the ‘‘“K component” is found to be half a magni- 


‘Lick Obs. Bull., No 464, 1934. 
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tude brighter than the uneclipsed portion of the A component at this time. The photo. 
graphic magnitude at the phase in question is 10.60. 

An alternative explanation, which must be explored, is that the eclipse is annular and 
not total. An annular eclipse seems, however, to be ruled out by the photographic light. 
curve; and Enebo’s visual light-curve is not good enough to carry much weight. From 
the observed duration of eclipse and totality, we find, for an annular eclipse, k = 0.81, 


TABLE 7 
ABSOLUTE DIMENSIONS OF RX GEMINORUM 


A4 Component KO Component 
0.8 Of 
2.20Ot 7.2 Ot 
5.500Tt 

Absolute visual magnitude........ +1.1* { +-2.1 * 

+2.7 


* From mass luminosity and the light-curve. 

t From mass function and &. 

¢ From mass function, radius of A4, and &. 

§ From mass of A4 (3.1) and mass function (0.018). 


From the depths of the two minima for the annular case, k = 0.93 (secondary minimum 
of depth 002) or k = 0.92 (secondary minimum inappreciable). The shape of the mini- 


mum is incompatible with such a high value of k; the largest value that can be deducedis 
0.72. This value is impossible because it leads to a negative inclination. For the darkened 
solution the matter is even worse. In other words, the system remains a puzzle. 


The writer extends his thanks to Dr. Struve for the opportunity to work at the 
McDonald Observatory in December, 1944. 
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THE ECLIPSING SYSTEM RY GEMINORUM 


SERGEI GAPOSCHKIN 
Harvard College Observatory 
Received August 3, 1946 


ABSTRACT 


The eclipsing system RY Geminorum has been studied from 13 spectrograms of the McDonald Ob- 
srvatory, 1667 photographic patrol plates of the Harvard College Observatory, and 73 visual observa- 
tions of Beyer. Outside eclipse the spectrum shows only one component, of class A2. During the total 
minimum the spectrum of the fainter, larger, and less massive is K2. Orbital velocities were obtained 
only for the brighter component (A2); the spectroscopic elements are: Vo = +9.0 km/sec; K = 30.0 
km/sec; a sint = 4.0 X 10°km = 5.75 ©; and f = 0.022 ©. The photographic elements for the (D) 
glution are: Max. = 8™52; A; = 2™80; A. = ?; k = 0.42; r4 = 0.086; rx = 0.205; 7 = 85°0; y = 64.6; 
Ly = 0.92; Lx= 0.08; Spx = K2 (predicted); and P = 9430090. The visual light-curve is less well 
determined, but the minmum is fully covered with observations. Its elements are: Max. = 8™45; 
4;= 200; Ao = ?; k = 0.35; ra = 0.079; re = 0.225; i = 80°4; La = 0.84; Le = 0.16; y = 42.5 
for (U); and the color index during the totality is 087. The comparison between the spectroscopic and 
the photometric observations indicates a satisfactory agreement. There is a definite eccentricity of 
the spectroscopic orbit. The comparison between the visual and the photographic dimensions (U) indi- 
cates the possibility of a thick atmosphere. Like RX Geminorum, this system consists of a main-sequence 
and a subgiant star. 


I, SPECTROSCOPIC OBSERVATIONS 


At the McDonald Observatory in December, 1944, I obtained 13 spectrograms of the 
eclipsing system RY Geminorum, of period 9 days. 

The spectrum of the brighter component presents a beautiful series of hydrogen lines 
and metallic lines, which are sharp and indicate giant characteristics. The spectrograms 


TABLE 1 
RADIAL VELOCITIES OF RY GEMINORUM IN KM/SEC 
JD Phase Mean H Feu Tiu Can 
(Km/Sec) 

0°349 —24.4 — 3.5 —34 —30 —30 
Serer 310 —12.1 —12.5 —32 + 3.5 —i1 
ence 425 — 3.5 — 5.5 + 1 —12 +1 
Ce ear 522 +10.0 — 3 +11 + 3 +22 
625 +23.0 +20.0 +12 +18 +41 
+50.4 +38 +43 +43 +60 
OS re .849 +32.7 +40 +37 +32 +66 
.. .058 + 0.5 — 1.5 +11 —17 +16 
. 266 — 9.8 — 6 —14 —32 
.022 + 6.4 — 3.5 4 — 6.5 +13 
123 — 0.6 + 3.0 —16.5 +1 +10 
0.329 —24.1 —16.5 —38 —32 —14 


* At minimum. 


were measured by the author and reduced by Mrs. Payne-Gaposchkin. The velocities are 
given in Table 1. The number of the plates is not large, but they are well distributed over 
the whole period of the variation. One plate is taken at the time of totality. The blue end 
of the spectrum has weakened greatly, and some new metallic lines appear in the red 


383 


d 
1S 
ned 
Bias 
th 
} 
ae 


384 SERGEI GAPOSCHKIN 


end: The spectrum is definitely of late type; my estimate is K2. The spectroscopic ele. 
ments are almost a copy of those of the system RX Geminorum, with which RY Geni. 
norum shows great similarity. They were derived by the method of Lehmann-Filhés (see 
Table 2). 


II. PHOTOGRAPHIC OBSERVATIONS 


The 1667 photographic estimations were obtained on the Harvard patrol plates, cover. 
ing an interval of forty-five years. The estimations of minima have been used to check 


TABLE 2 
SPECTROSCOPIC ELEMENTS FOR RY GEMINORUM 
asini = 4.00- 10°km Vo = + 9.0 km/sec 
= 0.022 K.= 30.0 km/sec 
Spectrum = A2 + K2 P = 9930090 


“40 1 


Fic. 1.—Velocities of RY Geminorum. Ordinates are radial velocities in km/sec, zbscissae are phases 
in units of period. 


15000 20000 25000 30000 
1.00 
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Fic. 2.—Phase of minimum of RY Geminorum as a function of Julian Day 


the period of 9430090. As is seen from Figure 2, no correction is needed. The photo- 
graphic light-curve reduced with the above period is given in Table 3. 

The magnitudes are on the International scale and were obtained by comparison with 
Harvard Standard Region C4. The comparison stars are given in Table 4. 

The photographic light-curve exhibits one of the largest ranges known. No secondary 
minimum is indicated by observations. Primary minimum is total. The solutions for this 
light-curve, using Scharbe’s method,! are given in Table 6. At the bottom of the minimum 
the (D) solution is better than the (U) (see Fig. 3). 


VISUAL OBSERVATIONS 


Asystem like RY Geminorum, with two widely different components and a rather long 
period, may furnish a good case for the study of a thick atmosphere, by comparing the 
relative dimensions obtained in the photographic and the visual wave lengths. Unfortu- 
nately, as with RX Geminorum, there is no well-observed visual light-curve. 


1S. Scharbe, Powlkovo Bull., 94, 31, 1923. 
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We have 73 visual observations by Beyer,” of which 63 are within a minimum. I have 
investigated the visual magnitudes of Beyer and found that they are close to the Inter- 
national scale: 8"35 of Beyer corresponding to 8"45 of Gaposchkin, and 10.35 of Beyer 
corresponding to 10"45 of Gaposchkin. The visual solutions were carried out with the 
same method as the photographic solutions. 


TABLE 3 
PHOTOGRAPHIC LIGHT-CURVE OF RY GEMINORUM 
Phase a IPg Phase be IPg Phase a I Ps 
40 8™52 27 8253 11 8™68 
060... | 63 8.52 47 8.52 8 9.53 
086... | 31 8.52 8.52 7 9.94 
* Computed with go = 1/P (JD —2400000). 
TABLE 4 
COMPARISON STARS (PHOTOGRAPHIC) 
Star IPg P.E. Star IPg P.E. 
8.68 + .03 10.65 + .02 
TABLE 5 
VISUAL LIGHT-CURVE OF RY GEMINORUM 
Phase Mag. No. Phase Mag. No. Phase Mag. No. 
8™45 1 10™50 6 9m56 4 
8.58 2 10.42 | 12 8.99 3 
8.89 6 10.43 8 8.54 2 
9.76 3 10.37 5 


?A.N., 234, 88, 1928. 
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IV. THE RELATIVE DIMENSIONS 


Comparing, now, the photographic dimensions with those from the visual light-curve, 
we find a small difference, as in the case of the very similar system RX Geminorum. The 
difference is of the same amount (4 or 5 per cent) for the two systems and in the same 
direction. It is small, and future photometric observations will be required to confirm its 
reality (see Table 6). 


20 26 24 26 28 20 22 24 26 98 
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Fic. 3.—Primary minima of RY Geminorum. Large dots represent normals; small dots, D solution; 
circled crosses, uniform solution. Ordinates are I Pg and I Pv; abscissae are phases. 
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Fic. 4.—Relative dimensions in the system RY Geminorum 


V. DISCUSSION 


The system RY Geminorum is a twin of the system of RX Geminorum. Both consist 
of widely different components: a main-sequence star and a subgiant, revolving around 
each other in about the same period. The similarity of the two systems can be seen from 
Table 7. Both systems show large ranges of light-variation: 202 for RX Geminorumand 
2™80 for RY Geminorum, with total eclipses. Also both have shown some indications 
that point to the existence of thick atmosphere. But there is one marked difference— 
the fact that, while the observed spectrum during the totality of RY Geminorum is in full 
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agreement with the color class, predicted from the ratio of surface brightness and the 
color index (k2 and K2), this is not the case for RX Geminorum (a5 and KO). The ob- 
served radial velocities of the centers of gravity of the systems are +9.0 for RY Gemi- 
norum and +33.5 for RX Geminorum. 

If we compare these two systems with that of AS Eridani* period of 2466, we see a 
great similarity in the spectrum, mass, radius, and luminosity of the brighter component 
(A); but, in respect to the fainter one, a conspicuous difference exists. Both fainter com- 
ponents in RX and RY Geminorum are of intermediate luminosity between giant and 
main sequence, while that of AS Eridani is on the main sequence. 


TABLE 6 
RELATIVE DIMENSIONS OF RY GEMINORUM 
PHOTOGRAPHIC PHOTOGRAPHIC 
VISUAL VISUAL 
U | U 
U D U D 
er 0.35 0.42 0.35 84.2 85.0 84.4 
93.9 64.6 42.9 | Spectrum;...... K5 
TABLE 7 
COMPARISON OF RX AND RY GEMINORUM 
RX Gem RY Gem RX Gem RY Gem 
12421 9430 4.07-10® km | 4.00- 10° km 
A4, KO A2, K2 3.1,0.620 | 3.5,0.700 
Mass function....... 0.0180 0.0220 2.2,5.50 | 2.7, G40 


The writer extends his thanks to Dr. Struve for the privilege of working for a 
month at the McDonald observatory. 


NotEe.—Since the paper was sent to the printer I have received a letter from Dr. McKellar 
at the Dominion Observatory (Victoria, B.C., Canada) in which he gives the following radial 
velocities of the star: 


JD Phase Velocity 
28199. 877 0.934 +23 km/sec 
209. 844 .010 —2 
486.944 .802 +56 
533.885 . 846 +32.7 
815.029 .077 + 7.0 
966.787 . 393 — 1.6 
982.775 .112 — 3.0 
29550.015 0.099 — 6.5 


Though the distribution of the velocities over the period is less favorable than my own, they 
confirm my results. An interesting plate is at phase 0.934 (during totality), 1°20" before the 
middle of the eclipse. My plate at the phase 0.951 was taken 2°20™ after the middle of the 
eclipse. The respective velocities are +23 and —20.3 km/sec. 


*S. Gaposchkin, Harvard Bull. (in press). 
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SOME RADIAL-VELOCITY AND LINE-INTENSITY MEASURES IN 
THE SPECTRUM OF 6 CANIS MAJORIS* 


ANNE B. UNDERHILL 
Yerkes and McDonald Observatories 
Received August 8, 1946 


ABSTRACT 


Coudé plates of 8 Canis Majoris taken in February and March, 1942, have been measured for radial 
velocity. The observations of February 13 and March 4 cover the 6-hour period nearly completely. The 
amplitude of the velocity-curve and y velocity vary, being 20 and 26.8 km/sec, respectively, on Febn- 
ary 13, and 7 and about 31 km/sec on March 4. A variable amplitude is in accord with previous observa- 
tions. A variable y velocity is consistent with Henroteau’s observations of 1922. A period of 0.252809 
day connects the present velocity-curves and that obtained by Struve and Swings in 1941. The times of 
maximum and minimum radial velocity agree with those predicted by Meyer’s formula. The line-quality 
changes follow the previously established period of 0.2513015 day. The line-sharpness period and the 


radial-velocity period are distinct. 
The central absorptions and equivalent widths of the O 1 and @ lines vary slightly over the 6-hour 


cycle, being largest at times of sharpest lines. The variation in the He 1 lines is less conspicuous, being 
of the order of the probable errors of the measurements. Curves of growth show that the O17 ions and 
Het atoms have zero or very small turbulent velocity. The populations of the Het levels and the On 
levels investigated are close to the values found by Unsdld for the BO dwarf, r Scorpii. 


The star 8 Canis Majoris is known! to show radial-velocity changes of varying ampli- 
tude with a period of about 6 hours. The quality of the spectral lines changes periodically 
also, the lines being sharp and deep at times and quite diffuse at other times. A varia- 
tion of 0.03 mag. in light in about 6 hours has been found by Fath.? 

The purpose of this investigation is to ascertain whether the period of the line sharp- 
ness and the period of the velocity changes are identical and to obtain an estimate of the 
turbulent velocity from curves of growth. The lines show considerable turbulent broad- 
ening. 

This investigation is a sequel to the work of Struve and Swings,’ which was based ona 
series of coudé spectrograms taken at the McDonald Observatory in January, 1941. 
Further plates of 8 Canis Majoris were taken by Swings and Struve at McDonald with 
the coudé spectrograph, series of plates, each extending over practically the full 6-hour 
period, being taken on February 13 and March 4, 1942. These plates form the basis of 


the following discussion. 
I. RADIAL-VELOCITY MEASUREMENTS 


All the plates were measured for radial velocity. Since the plates are too long to be 
accommodated all at once in the measuring engine, the plates were taken in sections, 
three regions on each plate being measured and reduced separately, and then a mean 
velocity for the plate was derived from the mean velocity of each region. To obtain 
the mean velocity of the plate, the mean velocity of each section was weighted according 
to the number of lines measured in that section. The velocities from each section agreed 
within their probable error. 

Table 1 gives details regarding the regions measured. Most of the lines used were 
O rand Het lines. The hydrogen lines were not used for determining the mean velocity 


* Contributions from the McDonald Observatory, University of Texas, No. 128. 
1 Henroteau, Lick Obs. Bull., 9, 157, 1918; Pub. Dom. Obs., 8, 31, 1922; Meyer, Pub. A.S.P., 46, 202, 
34. a 


2 Lick Obs. Bull., 17, 116, 1935. 3 Ap. J., 94, 99, 1941. 
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TABLE 1 
REGIONS MEASURED FOR RADIAL VELOCITY 
. Wave-Length Dispersion Source of Stellar 
Region Range (A/Mm) Wave Lengths 
3950-4036 At He, 2.01 Kiihlborn, Veréff. Univ. Sternw. Berlin-Babels- 
berg, Vol. 12, No. 1, 1938 
Re ey 4068-4188 At Hé, 2.51 Struve and Swings, A p. J., 94, 99, 1941 
4337-4529 AtHy,3.52 | Ibid. 
TABLE 2 
RADIAL-VELOCITY MEASURES 
Plate Date U.F. No. of Velocity Phase Line 
Cd 1942 (Days) Lines (Km/Sec) (P =07252809) | Quality 
|. Seale February 13.049 43 +25.2 0.000 3 
13.088 47 32.8 0.154 
|. ee 13.123 51 35.6 0.292 5 
13.150 56 36.6 0.399 2 
197 13.169 44 30.6 0.474 P 
13.190 45 0.557 1 
199 13.207 41 23.4 0.624 1 
13.228 42 17.0 0.707 a 
Macchia 13.257 50 17.0 0.855 2 
March 4.097 54 33.7 75.345 3 
220 4.123 47 75.448 4 
221 4.149 56 28.9 75.551 4 
ee 4.168 56 28.0 75.626 5 
223 4.201 53 28.9 75.756 4 
ee 4.239 37 27.8 75.906 1 
5.059 53 31.7 79.151 2 
5.077 49 32.5 79.222 3 
po ee 6.056 50 33.5 83.094 3 
pee 6.077 45 +32.5 83.177 2 
Struve and Swings’s Measures 
1941 
January 24.104 68 +38.17 


ERRATA 


Owing to an error in the reductions, the velocities obtained from plates Cd 193-201, inclusive, 
should be increased by 2.8 km/sec; those from plates Cd 219-228, inclusive, should be increased 
by 1.4 km/sec; and those from Cd 233 and 234 should be increased by 1.3 km/sec. The y 
velocities are now +29.6 km/sec on February 13, 1942, and +32.4 km/sec on March 4, 1942. 
Agreement is still not obtained with the Lick observations of a constant y velocity of +34 
km/sec for the system 8 Canis Majoris. 
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of each section. An attempt was made to measure all the lines listed by Kiihlborn‘ or by 
Struve and Swings? for each section; but the number of lines measured varies from plate 
to plate quite considerably, as the line quality changes over the 6-hour cycle. This varia- 
tion of line quality was noted by the previous investigators and is quite conspicuous on 
the present plates, since it is possible to measure many more lines at epochs of sharpest 
lines than at epochs of very diffuse lines. 

Table 2 gives the radial-velocity measures reduced to the sun. The y velocities ob- 
tained from these short-period radial-velocity changes are +26.8 km/sec on February 
13, 1942; +31 km/sec on March 4, 1942; and +31 km/sec on January 24, 1941. Struve 
and Swings’s measures are reproduced here, as they will be used to supplement the new 
material. The observations of January 24, 1941, and March 4, 1942, do not extend over 
quite a whole cycle, hence the y velocity is estimated from the appearance of the section 
of velocity-curve available. 

A varying ¥ is in accord with Henroteau’s! observations of 1922, which indicate that 
there is a short-period velocity variation of changing amplitude superimposed over a 
systemic velocity, which changes with a 42-day period. If the y velocity from Struve and 
Swings’s measures is put on the descending branch of Henroteau’s curve for y velocity, 
the more recent measures fall on the curve pretty well. There is not sufficient evidence 
to define the curve any better than Henroteau has done. The present observations are 
not consistent with a constant y velocity of 34 km/sec, as reported by Meyer! from a 
study of the plates obtained at the Lick Observatory. 

It was found that a period of 0.252809 day satisfactorily connects the short-period 
velocity-curves over about 1600 cycles, which is the scope of the observations discussed 
here. Time, fo, is taken arbitrarily at 1942, February 13.049, U.T., and the phase of each 
plate for a period of 0.252809 day is given in the sixth column of Table 2. Figure 1 
shows a plot of the three most nearly complete series of measures. Maximum velocity 
occurs at phase 0.335. The long-period change in amplitude and y velocity of these 
short-period velocity-curves is evident. According to Meyer, the amplitude changes with 
a period of 49.1 days. 

In the seventh column of Table 2 are given visual estimates of the line quality. 
These estimates are based on a scale of 5, the value 1 being given for broad, diffuse lines, 
which show no central core; 5 to lines which have a sharp, narrow core and very small 
wings. Value 3 is halfway between 1 and 5, the lines showing some core, but also fairly 
diffuse wings. Values 2 and 4 are intermediate between 1 and 3, and 3 and 5, respective- 
ly. No definite asymmetries were noticed in the lines. 

Small arrows in Figure 1 indicate the epochs of sharpest lines. These times are con- 
sistent with the period of 0.2513015 day, whict was established by the Lick! observa- 
tions. The times of sharpest lines are not in phase with the epochs of maximum velocity. 
Plate XIV shows the group of O 0 lines near Hy and illustrates the change in appearance 
of the lines. On February 13, sharpest lines occurred at phase 0.292 near maximum 
velocity, whereas the lines are very diffuse at phase 0.624, intermediate to low velocity. 
On March 4—19 days later—the situation is reversed, sharpest lines occurring at a time 
of low velocity and quite diffuse lines occurring at the highest velocity measured. This 
evidence shows that the short-velocity period and the line-width period are distinct, 
thus confirming the work of Henroteau and that of Meyer. 

Figure 2 shows the observed short-period velocity-curves and the curves computed 
from the formula developed by Meyer' to represent the short-period variations. In order 
to make the comparison as convenient as possible, the calculated velocities have been 
plotted around the observed y-axis in each case rather than around the constant y-axis 
of 34 km/sec found by Meyer. It is seen that Meyer’s formula predicts time of maximum 
and minimum very well in these three cases, though the fit in the intermediate part of the 


4 Veriff. Univ. Sternw. Berlin-Babelsberg, Vol. 12, No. 1, 1938. 
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curves is not quite so good and the amplitudes are not quite right. Neither of the periods 
of 0.2513015 day and 0.2500222 day, which are combined by Meyer, will fit the radial- 
velocity observations when used singly. 

The radial-velocity measurements suggest that 8 Canis Majoris is a pulsating star 
of short period. Upon these variations is superimposed a velocity variation of period 
of about 42 days, which could be interpreted as binary motion around an invisible com- 
panion but also might be similar to the long-period radial-velocity variation shown in 
some RR Lyrae-type stars to be superimposed over the short-period fluctuations. It 
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Fic. 1.—Velocity-curves of 8 Canis Majoris 


seems likely that the difference in period between the radial-velocity changes and the 
changes in line shape, as well as the varying amplitude of the radial-velocity-curves, is 
connected with the long-period variation. 


II. EQUIVALENT-WIDTH MEASUREMENTS 


Equivalent widths of the He 1 lines, the H lines, and the lines of the 3s — 3p transition 
array of O 11 were derived from the most suitable plates. The equivalent widths of the 
0 and the H lines are greatest at the times of sharpest lines. The change in equivalent 
widths of the He lines is indeterminate, being of the order of the probable error of the 
measures. The measured central absorptions of the lines were found to vary over the 
cycle, being greatest at epochs of sharpest lines. The effect is most noticeable for the 
0 11 lines and less conspicuous for the He 1 and H lines. 

The individual measures of equivalent width and central absorption are given in 
Table 3. The equivalent widths are given in milliangstroms and the central absorptions 
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in percentages. No attempt has been made to correct these central absorptions for in- 
strumental profile; however, for any given line they are indicative of the true change in 
central intensity, though direct comparison of lines of widely different wave length and 
total absorption would be hazardous. The equivalent width measure is W; the central ab- 
sorption is A,. The O 11 lines are arranged in multiplets of the 3s — 3p transition array; 
and the Het lines are arranged according to series. The measures for plates Cd 221 and 
195 are put side by side, since the lines are sharp in these cases. The lines are visually 
quite diffuse on the other plates, Cd 196, 197, 200. The last two columns of the table 
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25- 


U-T 
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Fic. 2.—Comparison of observed velocity-curves and those predicted by Meyer’s formula. The 
calculated values have been plotted around the observed + axis, instead of around y = 34 km/sec, as 
given by Meyer. 


give, respectively, the mean of the Cd 195 and 221 measures, and the mean of the 
measures from plates Cd 196, 197, 200. 

The measures for the O 11 line \ 4367 are very uncertain, as this line was very weak, 
and are not used in constructing the curve of growth. The measures for the He 1 line 
3965 and the O 11 line \ 3973 are also rather uncertain, as these lines lie in the wings 
of He and most of the plates were rather thin in this region. When a line fell in the wing 
of a H line, the wing was taken as the continuum for that line. A correction for the blend- 
ing effect of hydrogen was not attempted. The lines 3965 and \ 3973 were used in 
constructing curves of growth; but, in placing the theoretical curve, little weight was 
given to them. 


5Zs.f. Ap., 21, 22, 1942; Revised Multiplet Table, p. 8, 1945. 


R. 
sion: 
ines 
lin 
as 
the 
in B 
HD 
— 
) 
3983 
30 432¢ 
4337 
é P 4345 
4367 
a 
4347 
20 4351 
441: 
4417 
oN MARCH 4 1942 
30 ° 400! 
414 
416' 
Bays 443 
3071 
410 
434 
i 
4 lin 
cal 
the 
cul 
: ber 
fitt 
He 


r in- 
ge in 

and 
] ab- 
ray; 

and 
lally 
able 


8 CANIS MAJORIS 393 


These line-intensity observations are directly comparable with those obtained by 
R. M. Petrie® in a study of HD 199140, a 6 Canis Majoris star, and confirm his conclu- 
sions about the variation in central intensity and slight change in equivalent width of the 
lines over the short period. No conclusion can be drawn from the present observations 
as to the reality of a change in spectral class. Furthermore, the present data show that 
the period of line-quality change is distinct from that of the short-period velocity changes 
in 8 Canis Majoris, although a similar difference does not appear to be reported for 
HD 199140. 


TABLE 3 
EQUIVALENT-WIDTH MEASURES 
INDIVIDUAL VALUES MEAN VALUES 
N Ete-| cq 221 Cd 195 Cd 196 Cd 197 Cd 200 Sharp Diffuse 
MENT Lines Lines 
3973 Ou 85) 17 88) 17 117| 16 88} 16 102) 16.5 
3983 63) 14 54! 10 53] 10 57| 12 10 
4317 210} 28 181} 29 122} 19 20 142} 22 195} 28.5} 130} 20 
4320 178) 30 176} 29 146} 21.5} 162) 21 139} 22 177| 29.5} 149) 21.5 
4326 13 109} 15 56} 10.5} 96) 12 67; 10.5} 84) 14 73] 11 
4337 65} 10 93) 15 67| 11 8 9 79| 12.5} 72) 9 
4346 166} 25.5} 160) 25.5) 108) 21 179} 28 123} 15 163} 25.5} 137] 21 
4349 277| 38 290) 40 258} 32 290} 33 227) 29 284) 39 258) 31 
4367 444 7 65) 8 27| 4 36) 7 46) 6 
4347 128) 16 174) 23.5) 101] 14 124) 16 149} 18 151} 20 141| 18 
4351 169} 25 174) 31 19 141; 19 156} 21 28 156} 30 
4415 36 39.5} 224! 27.5) 304! 33 237| 28 260} 38 255} 29.5 
4417 208) 34 190} 34.5) 182} 25.5} 241) 27 200} 22 199) 34 208) 25 
4452.... 124) 15 127| 17 84) 14 93) 11 72} 11 126) 16 83) 12 
4026....| Het} 546) 33.5) 610) 36 662) 41 738} 37.5) 493) 31.5} 578) 34.5) 631) 36.5 
4471.... 57 56 718) 44 770} 52 689! 46 842) 56.5) 726) 47 
4009 265) 19 260} 21.5} 200} 20 216) 18 184} 16 262} 20 200} 18 
4144 28 333) 31.5) 398) 30 358} 28 366} 29.5} 315} 30 374| 29 
543} 42 470} 34 544) 39.5) 487) 36 543) 42 500} 36.5 
3965 113) 21 112) 14.5) 133} 21 113} 19.5} 122) 18 
4169 97| 14 116) 15 114) 14 94) 12 141} 14.5} 107) 14.5} 116) 13.5 
4437... 137; 23.5) 112) 17 102) 14 97) 14 116} 15 124) 20 105} 14 
3970....| H 1649) 40 { 1725) 1822} 43 | 1471) 38 1687) 40.5) 1646) 40.5 
2610) 52 | 2310) 48 | 2124) 45 | 1925) 43 | 2610) 52 | 2120) 45 
4340 3057; 54 | 2980) 57 | 2522) 48 | 3434) 50 | 2198) 48 | 3018) 56 | 2718) 49 


Curves of growth were constructed from the O 1 lines (Fig. 3, a) and the Het 
lines (Fig. 3, 6) by plotting log W/) against log gif — (5040/7) xi. Here g: is the statisti- 
cal weight of the lower level of the line, Xi the excitation potential of that level, and T 
the temperature. A temperature of 25,200° is assumed. The solid curve is the theoretical 
curve of growth computed in each case from the approximate equations quoted by Gold- 
berg.” Mean equivalent widths from plates when the lines are sharp and from plates 
when the lines are diffuse have been plotted separately. The theoretical curves can be 
fitted to the observations without any vertical shift, thus indicating that the Om and 
Het atoms have small turbulent velocities. 


® Pub, A.S.P., 9, 53, 1939; J.R.A.S. Canada, 40, 149, 1946. 
7 Ap. J., 89, 623, 1939. 
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The total width in angstroms of the O lines was measured on the profiles, at a depth 
equal to half the central absorption. The mean width at half-absorption of all the lines 
corresponds to a velocity range of 45 km/sec. The width is slightly more when the lines 
are diffuse and slightly less when the lines are sharp; but the accuracy of the measure- 
ments is not sufficient to warrant treating the two cases separately. Greenstein® has 
shown that the effective resolving-power of the coudé spectrograph is about 35,000, 
corresponding to a Ad of 9 km/sec at 4300 A. The thermal velocity of the O11 ions is 
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Fic. 3.—Curves of growth: a, 3s—3p transition array of O 0; }, all available He 1 lines 


5.1 km/sec at 25,200°. Thus some other broadening agent beyond thermal motion and 
finite resolving-power of the spectrograph must be found to account for the observed 
width of the lines. It is suggested that motion of optically thick sections of the atmos- 
phere causes the broadening which is observed. The intensity of the line is not affected 
by these mass motions, since, within the range of depth in the atmosphere in which the 
line is formed, the atoms have only thermal motions. However, when the integrated 
starlight from all over the disk is observed, the parts of the line formed in each element 
of atmosphere are seen displaced to the violet or to the red according to whether the ele- 
ment in which the line was formed was expanding or contracting, and hence a broadened 
profile is obtained. The curve of growth is not affected, since the absorption coefficient of 


8 Ap. J., 95, 194, 1942. 
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the atoms is not affected by these superimposed motions, any more than it is when the 
star is in rotation. 
Horizontal displacements of the curve of growth give 


No we? 1_ 


and thus enable us to estimate the population of the lower levels of the lines. The qual- 
ity NV. is the number of atoms in the requisite state of ionization; b(7) is the partition 
function; and v is the thermal velocity of the atoms. The results are given in Table 4. The 


TABLE 4 
POPULATION OF LEVELS 
Het Ou 
Level log NV Level log Nf3s—sp 


total f-value of the 3s — 3p transition array of O 1 is not known, hence only log Nf,,_;, 
can be determined. It seems likely that f,,_;, is close to 1. 

These populations are comparable to those found by Unsdld> for r Scorpii. He found 
log N = 14.95 for the 2°P level of He1, and log N/f,,-;, = 14.26 for O 11. 


This investigation was undertaken at the suggestion of Dr. O. Struve and owes much 
to his guidance. The author is grateful for many stimulating discussions with Dr. Struve 
and with Dr. J. L. Greenstein of problems arising in the course of the investigation 
and expresses her thanks for their help. 
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SPECTROGRAPHIC OBSERVATIONS OF 
RY PERSEI AND RZ OPHIUCHI* 


W. A. HILTNER 
Yerkes and McDonald Observatories 
Received August 5, 1946 


ABSTRACT 


The eclipsing variables RY Per and RZ Oph have been investigated spectrographically. The velocity- 
curves for hydrogen and helium in RY Per are not identical. Helium gives a rotational velocity approxi- 
mately twice that of hydrogen. Satellite lines of hydrogen before first contact suggest that a stream 
of hydrogen gas flows from the FS star to the smaller B4 star. RZ Oph shows hydrogen emission at all 
phases. The emission has its origin in a gaseous ring around the smaller F8 primary star. A significant 
er was not obtained; however, an approximate upper limit for the semi-amplitude has been 
established. 


RY PERSEI 


The variable star RY Persei is an eclipsing binary with a period of approximately 1 
week. Its variability was discovered by Mme Ceraski in 1906. A number of investigators 
have observed the light-curve; Nijland! established a period of 6.863571 days, and Shap- 


TABLE 1 
WAVE LENGTHS OF STAR LINES FOR RY PERSEI 


B4 Star FS Star 
3819.61 Feu+Cri+Mgt....| 4351.84 
4471.48 


ley? and Fetlaar* have made solutions based on Nijland’s observations. Recently Wal- 
bach and Gaposchkin‘ published a solution based on a light-curve obtained from esti- 
mates on patrol plates. The photographic range observed by them is 2.01 mag. at pri- 
mary minimum, and a secondary minimum was reported with a depth of 0.01 mag. These 
observers also noted a periodic change in the period; a sine curve with a semi-amplitude 


* Contributions from the McDonald Observatory, University of Texas, No. 129. 
14.N., 176, 172, 1907; 183, 281, 1910; 231, 207, 1927; B.A.N., 2, 231, 1924. 
2 Contr. Princeton U. Obs., No. 3, 1913. 

3 B_A.N., 3, 195, 1926. 

4 Pub. A.A.S., 10, 72, 1940; see also Harvard Bull., No. 915, 1941. 
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SPECTROGRAPHIC OBSERVATIONS 397 
TABLE 2 
RADIAL VELOCITIES OF RY PERSEI 
RapiAL VELOCITIES 
(In Ka /SEc) 
PHASE 
PLATE DaTE U.T. (In Days) 

H Het F Star Cau K 
am 6710. . 16.468 0.169 —39.4 $90.0 

ant 16.505 0.206 —43.8 + 0.8 

een Feb. 25.156 0.218 —51.3 —27.1 

>... 4.107 0.305 —28.8 — 7.9 

5468..... 11.070 0.405 — 1.0 — 6.7 

Ors 5469..... 11.082 0.417 +12.3 —25.0 

5470..... 11.097 0.432 — 3.9 2.9 

5471..... 11.111 0.446 | + 2.0 67 

CQ 4037..... 1944 Dec. 5.089 0.514 —12.9 —16.5 

Gf/2 5374..... 1945 Mar. 5.071 1.269 +13.8 —15.4 

5.090 1.288 — 6.1 —19.3 

12.070 1.405 — 8.3 —10.3 

5489..... 12.092 1.427 —19.3 —10.3 

5490..... 12.102 1.437 —20.2 —14.1 

CQ 4043..... 1944 Dec. 6.384 1.809 — 0.8 > | —20.1 

Gf/2 5502..... 1945 Mar. 13.095 2.430 10.9 —20.8 

CQ 4047..... 1944 Dec. 7.033 2.458, —20.3 ow —11.3 

4051..... 7.184 2.609 —17.1 —17.4 

Gf/2 5330.....| 1945 Feb. 28.135 3.197 —15.1 ca ee —12.2 

28.156 3.218 —14.0 —12.2 

28.174 3.236 — 9.9 — 4.6 

5396..... Mar. 7.106 3.304 —16.4 + 0.2 

7.136 3.334 — 8.3 SRB + 0.2 

. Feb. 22.072 3.998 — 1.7 — 9.0 

Gf/2 5334..... Mar. 1.166 4.228 — 3.2 —12.1 

5406..... 8.069 4.267 — 5.2 — 7.2 

S408... 8.090 4.288 +15.2 — 3.5 

CQ 3403..... 1944 July 4.462 4.749 +13.3 oe OMT. ck —13.3 

1945 Feb. 23.068 4.994 +36.2 — 9.0 

4360. .:.. 23.092 5.018 +32.4 — 6.0 

Gf/2 5424..... Mar. 9.071 5.269 +39.0 | —10.9 

| COM 1944 Dec. 3.204 5.492 +26.2 — 3.5 

1945 Feb. 24.069 5.995 +20.9 — 8.9 

e 4370; 24.094 6.020 +32.0 — 5.9 

3.113 6.175 +46.3 — 4.3 

10.095 6.293 +56.0 — 6.9 
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398 W. A. HILTNER 
TABLE 2—Continued 
| | RapIAL VELOCITIES 
(In Km /SeEc) 
| | PHASE | 
PLATE | DATE | jai (In Days) | 
| | Het F Star | CanK 
CQ 4033..... | 1944 Dec. 4.099 | 6.387 | +36.5 | — 24.3].......... — 7.0 
4034..... 4.124 6.412 +31.5 =) ~ 0.9 
Gf/2 6617..... | 1945 Sept. | 25.434 6.589 +62.5 oS ee + 6.7 
6619..... | 25.472 6.627 +57.2 + 6.7 
| 


of 0.5 hour and a period of 19 years agreed with the observations. The following relative 
elements of the system were obtained: 


L,=0.157 ,5 r2=0.122A, 
r,=0.197 A, i= 89°0. 
L. = 0.843 


RY Per was placed on the observing program of the 82-inch reflector of the McDonald 
Observatory in July, 1944, and 62 spectrograms were obtained from July 4, 1944, to 
October 16, 1945. Two arrangements of the Cassegrain spectrograph were employed: 
two quartz prisms and a 500-mm camera (designated ‘“‘CQ’’), giving a dispersion of 
55 A/mm at Hy, and two glass prisms and a 180-mm Schmidt camera (designated 
“G f/2’’), giving a dispersion of 75 A/mm at Hy (see Table 2). 

The spectrum of RY Per at maximum light is B4 with moderately diffuse lines (see 
P|. XV). The helium lines probably appear more diffuse than the hydrogen lines, but the 
comparison is difficult. On the spectrograms taken with the quartz prisms the hydrogen 
lines were usually measured for radial velocity to H13, beginning with Hy, while the 
spectrograms obtained with the glass prisms were frequently measured to 10, begin- 
ning with Hf. The line He was always excluded because of a strong interstellar Ca 11 
H line. Four helium lines were normally measured. The wave lengths employed for de- 
termining radial velocities are listed in Table 1. 

Table 2 gives the observed radial velocities. The phases reckoned from primary mini- 
mum were computed by the formula 


Primary minimum = JD2418216.701+6.863571E, 


- given by Nijland.® The velocities for hydrogen and helium are listed in separate columns. 
The observed radial velocities are plotted in Figures 1 and 2. It is clearly shown that the 
two velocity-curves are not in agreement. Graphical solutions of the two velocity-curves 
yielded the results shown in the tabulation on facing page. 

The radial velocity of the secondary component, observable at primary minimum, has 
been measured on seven spectrograms. Its spectral class is F5, with giant characteristics; 
the \ 4077 line of Sr 1 is rather strong relative to Fe 1 4063. The assumed wave lengths 
are also listed in Table 1 and the observed radial velocities are recorded in Table 2 and 
plotted in Figures 1 and 2. 


5 The subscript 1 refers to the star in front at primary minimum. 
6 A.N., 231, 207, 1927. 
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There are several interesting features in the spectrum and the velocity-curves of 
RY Per. First, the rotational effect of hydrogen and helium is different, with the latter 
element showing nearly twice the rotational velocity. A similar phenomenon was ob- 
served in RZ Scuti by Neubauer and Struve.’? In RZ Sct the range for helium was 
250 km/sec and that for hydrogen, 100 km/sec; and in RY Per helium and hydrogen 
show rotational velocities of approximately 200 km/sec and 110 km/sec, respectively. 
In RY Per, as in RZ Sct, it is probable that hydrogen forms a slowly rotating envelope 
about the reversing layer of the star. Second, the hydrogen velocity-curve from phase 
6.15 to phase 6.5 is anomalous. The observed radial velocities are systematically large. A 
similar phenomenon was also observed in RZ Sct, but in that case hydrogen and helium 


| 
| Hydrogen Helium 
6.863571 (assumed) 6.863571 (assumed) 
+4.5 km/sec — 26 km/sec 
27 km/sec 36 km/sec 
2.610% km 3.4X 10° km 


RY Persei— Hydrogen lines 


F/2G 


| co ° 


° 2 4 6 8 10 Phase in days 


Fic. 1.—Velocity-curve of RY Persei from the H lines 


were similar. In RY Per it is clear that the two velocity-curves are different at this 
phase. Third, the velocity-curves for the two elements are different in amplitude, and 
still more pronounced is the difference in y velocities. It is well established that the 
hydrogen-curve normally falls above the helium-curve. This same phenomenon was ob- 
served by Deutsch® in BD + 55°616. He found that the y velocity for hydrogen was 
14 km/sec more positive than that of He 1 and that the semi-amplitude for hydrogen was 
38.4 km/sec, in contrast to 65.6 km/sec for that of Het. It is also observed that in 
RY Perand BD + 55°616 the velocity-curves for the two elements in question are simi- 
lar when the brighter components are receding from the observer. To interpret the dis- 
crepancies in the observed velocity-curves it appears necessary to assume that the 
brighter, smaller component is surrounded by an asymmetrical envelope and is approach- 


7 Ap. J., 101, 240, 1945. 8 Ap. J., 102, 496, 1945. 
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ing the star with respect to the He 1 atoms. That some such mechanism is present in the 
case of RY Per is substantiated by three spectrograms obtained on March 10, 1945, 
which suggest that a stream of hydrogen flows from the F5 star to the B4 star. On the 
above date, between phases 6.27 and 6.29 (0.10 days before first contact), the hydrogen 
lines showed a second component displaced to the red, resembling the satellite lines ob- 
served in 6 Lyr,® which were interpreted in that case in terms of a flow of material from 
the smaller F star to the larger B9 star. The interstellar K line has a velocity of —9.4 
km/sec. 


RY Persei— Helium lines 
km/sec T T T T T T T T T T 
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Fic. 2.—Velocity-curve of RY Persei from the Het lines 


RZ OPHIUCHI 


The variability of RZ Oph was also discovered by Mme Ceraski, in 1905. A period of 
261.9 days was established by Beyer'® and others. The depth of the total primary mini- 
mum is 0.80 mag. The eclipse lasts 16 days, or 0.061 P, and totality lasts 8 days, or 
0.031 P. Outside of eclipse, mag. 10.0, the spectrum is that of an F8 supergiant, similar 
to y Cyg. At primary eclipse the spectrum closely resembles that of  Cyg, a K5 super- 
giant. In full light hydrogen emission at H is present as a double line, with a separation 
of the components of 284 km/sec. This emission undergoes an eclipse similar to that in 
some other eclipsing binaries with hydrogen emission. The behavior of 8 emission near 
and during primary eclipse is given in Table 3. The spectrograms were obtained at the 
McDonald Observatory with the Cassegrain spectrograph and have a dispersion of 
120 A/mm at H8. This dispersion is sufficient easily to resolve the components on East- 
man 103a-O emulsion (see Pl. XVI). 

The observations of the eclipse of the ring suggest that the inner radius of the ring 


® See O. Struve, Ap. J., 93, 104, 1941, 10 4.N., 234, 94, 1928. 
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is approximately equal to the radius of the secondary star and that the outer radius is 
about 50 per cent greater. 

A total of thirty-seven spectrograms of RZ Oph were obtained during one cycle, of 
which all but seven were outside of eclipse. In addition to estimating the character of HB 
emission, the spectrograms were measured for radial velocities. The pertinent data con- 
cerning the spectrograms are given in Table 4, along with the measured radial velocities. 
The phases were computed by the formula given by Beyer.’ 


Primary minimum = JD2418630.65 + 261.943E . 


TABLE 3 
Hg EMISSION IN RZ OPHIUCHI 
Computed Phase Character of H@8 Emission 
Both components present 
Only red component present 
Red component much stronger 
Red component much stronger 
Components about equal 
Only violet component present 
Only violet component present 
Violet component present and 
possibly the red component 
TABLE 4 
RADIAL VELOCITIES OF RZ OPHIUCHI 
| | Radial || Radial 
Plate Date | UT Phase Velocity |} Plate Date UT Phase | Velocity 
(G f/2) (1945) (In Days)| (In (G f/2) (1945) (In Days)| (In 
| | | Km/Sec) Km/Sec) 
SSO ion | Mar. 3.500 | 52.14 | + 9.6 || 6518....... Sept. | 20.101 | 252.74 | + 1.7 
§.489 57.13 | + 2.6 21.091 | 253.73 | + 0.1 
SS: | eee 17.494 | 66.14 | — 5.7 || 6562..... 22.252 | 254.90 | — 7.5 
SOO 26.469) 75:12.) 23.155 | 255.80 | + 4.8 
30.462 | 79.10 | + 0.9 || 6588....... 24.122 | 256.76 | —13.2 
ST acca ss May | 30.414 | 140.06 | — 6.1 || 6599....... 25.103 | 257.75 | — 6.6 
SIO is was June | 12.394 | 153.04 | — 2.5 || 6624....... 27.094 | 259.74 | + 4.9 
17.279 | 157.92 | — 1.7 ... 28.097 | 260.74 | + 2.6 
Cy Ce 21.408 | 162.05 | + 2.2 || 6645....... Oct. 1.122 1.82 | + 6.9 
LS Seer July 18.356 | 189.00 | — 3.7 || 6666....... 7.169 7.87 | — 4.3 
Oy as ators 21.207 | 191.85 | + 0.3 || 6673....... 12.088 | 12.7 — 5.4 
31.272 | 208.92 | =. 13.065 13.77 | +13.6 
ere Aug. | 13.262 | 214.90 | + 5.9 || 6689....... 14.107 14.81 | + 8.8 
14.246 | 215.89 0.0: || 15.093 | 15.79 | + 0.4 
16.292 | 217.94 | — 1.7 6700....... 16.075 | 16.77 | — 0.5 
| eee 30.208 | 231.85 | — 4.8 || 6736...... 18.076 18.77 | — 6.3 
Ee) ree Sept. 7.215 | 239.86 | + 6.7 || 6789....... 23.075 | 23.77 | +14.5 
1 15.154 | 247.80 | + 1.6 || 6824....... 26.055 | 26.75 | — 0.8 
: 18.213 | 250.86 | +10.6 
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Although the spectrum has exceptionally sharp lines, the resulting measurements are dis- 
appointing as can be seen in Figure 3. It is obvious that we have no information concern- 
ing the velocity-curve except that an approximate upper limit can be established for the 
semi-amplitude. However, if no orbital motion is assumed, the resulting probable error of 
a single observation is + 3.8 km/sec, which is very similar to the value +3.3 km/sec that 
Struve" obtained from a series of spectrograms of a Per with the same spectrograph. 
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Fic. 3.—Radial velocities of RZ Ophiuchi 


This suggests that the semi-amplitude is certainly not greater than 5 km/sec. The upper 
limits of a sin i and f(m) then take on the values 18 X 10° km and 0.0030, respectively. 

A discussion of eclipsing binary systems with gaseous hydrogen rings was recently 
published by Struve,'® and RZ Oph was among the stars discussed. 


It is a pleasure to thank Drs. Struve, Sahade, and Edmondson for obtaining some of 
the spectrograms. 


1 Ap. J., 102, 125, 1945. 
12 Ap. J., 103, 95, 1946; Observatory, 66, 208, 1946. 
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THE RATIO OF INTERSTELLAR ABSORPTION TO REDDENING 


Jesse L. GREENSTEIN 
Yerkes Observatory 
Received July 31, 1946 


ABSTRACT 


The wide-base-line photoelectric six-color photometry of Stebbins and Whitford provides a new 
determination of the variation of space reddening with wave length; the derived variation is near 1/°-5*, 
The ratio of photographic absorption to short-base-line photoelectric colors (C; scale) lies between 10 
and 15. The curvature of the energy-curves of the reddened B stars with respect to unreddened B stars 
may produce an underestimation of the reddening if measured in color classes; longer-base-line color 
scales, such as the Greenwich gradients, and the V—I scale show larger reddenings. An observational 
technique for the discovery of reddened stars is based on the apparent discrepancy between long- and 
short-base-line colors. The present absorption law is interpreted with reference to the theory of absorp- 
tion and scattering by small particles; a direct determination of the ratio of absorption to reddening may 
permit a decision as to whether the particles are dominantly metals or dielectrics—the latter now seem- 
ing more probable. 


The photoelectric six-color photometry!” by Stebbins and Whitford provides a major 
forward step in the determination of the ratio of interstellar absorption to reddening. 
The large wave-length range covered provides, in itself, a definite minimum value. 
Extrapolating observed absorptions for estimation of the absorption at infinite wave 
length now involves smaller uncertainty. Certain difficulties arising in the use of short- 
base-line color excesses are also clarified. In their first paper,! Stebbins and Whitford 
derive the variation of interstellar absorption with wave length. Since, to a first approxi- 
mation, this variation is nearly proportional to 1/A, they compute the deviation of the 


TABLE 1 
POINTS MEASURED IN SIX-COLOR PHOTOMETRY 
U Vv | B | G | R I 
Ridin eke 3530 4220 4880 5700 7190 10300 A 
2.83 2.37 2.05 1.75 1.39 0.97 


observed energy-curves of reddened stars from those expected on the basis*of the 1/A 
relation. They compare a reddened and an unreddened B star by forming the differences 
between the observed magnitudes at six colors. Their six effective wave lengths are given 
in Table 1. The stars observed are of spectral types O5—B3; many of the reddened stars 
are supergiants. Their reductions involve the use of colors measured at the ultraviolet 
point, U. Since U may be affected by the Balmer absorption (or emission) lines and 
continuum, I have reinvestigated the reddening as derived from the points V to I, inclu- 
sive. The Paschen lines and continuum should be substantially weaker than the Balmer 
lines, and there is little reason to expect effects at point I. There are no other lines 
strong enough to affect materially the measured energy of the star within the wide filter 
bands employed. Absolute-magnitude effects should be small, although not absolutely 
excluded, in the range V — I. 

The reduction of the data on the V — I base line follows directly the method used 
by Stebbins and Whitford! in their Table 5. The magnitude differences between the red- 


1Ap. J., 98, 20, 1943. 2 Thid., 102, 318, 1945. 


403 


; 
= 
= 
a 
: 
: 


404 JESSE L. GREENSTEIN 


dened and the unreddened star are reduced by a scale factor to the values that they 
would have had if the reddening (and possible temperature difference) had been such as 
to give a V — I color difference of 1.00 mag. (This assumes a proportionality of the ab- 
sorptions, or the existence of a single functional form for the reddening law.) The ex- 
pected differences, if the reddening varied as 1/A, are then subtracted and the residuals, 
D, tabulated. The mean residuals, grouped by spectral type, are given in Table 2, which 
gives the number (NV) of reddened stars; the mean observed color excesses E(C;) on 
the short-base-line photoelectric scale ;’ the actual mean observed 6(V — I) for each group; 
and the values of D, the deviation from the 1/) law, as reduced to the case 6(V — I) = 1.00 


TABLE 2 
DEVIATIONS OF THE REDDENING FROM THE 1/A RELATION 


Tyre N | EG) | 

| | V | B | G R I 
| 7 | +0719 | +103 | 02000 | +07035 | +0"033 | +02035 | 0=000 
8 | 40.34) +42.06 | 0.000 | +0.044 | +0.055 | +0.041 | 0.000 
6 | +0.28 41.61 | 0.000 | +0.046 | +0.066 | +0.012*| 0.000 
6 | +0.30 +1.78 | 0.000 | +0.061 | +0.054 | +0.032 | 0.000 
Mean Dj............ | 0000 | +0046 | +0049 | +0031*] 02000 
Reddening for 1.000 0.771} 0.557 | 0.300 | 0.000 
Observed reddening for, | | 
V—I=1.00........ 1.000 0.817 0.606} 0.331* 0.000 


* Two B1 stars have discrepant measures at point R, which reduce the value of D far below the mean. If these stars are 
omitted, D for B1 stars at point R becomes +0™032, and the mean D for all stars becomes +0™036. This change would make the 
curve for the variation of D with 1/\ considerably smoother. 


+ The method of choosing pairs of stars is somewhat arbitraryand may be criticized because it involves a comparison of 
stars of high and low luminosity. Since it is theoretically possible that the energy distribution in a supergiant may show intrinsic 
differences (and curvature) when compared to that in a dwarf, a part of the observed curvature might not be due to space red- 
dening. To investigate this possibility, six relatively slightly reddened supergiants of types B1 and B2 were selected as well as 
four highly reddened supergiants. The luminosity classifications were from unpublished data by Dr. W. W. Morgan. The mean 
V—I color of the first group was —1.78 mag. and of the second group +0.94 mag. The reduction to determine the law of red- 
dening was then carried through as before, and the deviations, D, from the 1/A law were reduced to a reddening 6(V — I) = 
1.00 mag. The results were as follows: 


Point V B G R I 
D +0m054 +0™066 +0m055 0™000 


The new values of D are similar to those found in Table 2, being larger by an amount slightly in excess of their probable error. 
Any intrinsic curvature of the energy-curves of the supergiants should cancel out in this comparison. We may therefore safely 
exclude absolute-magnitude effects as the source of the observed curvature. 


mag. The ratio of (V — I)/E(C)) is found to be 5.8; this is an underestimate because of 
the systematic errors introduced by the method of selection of the pairs of stars. 

The observed deviations in Table 2 are similar to those found by Stebbins and Whit- 
ford,! but on a smaller scale because of the reduced wave-length range covered. The 
agreement shows that the Balmer absorption has little systematic effect on the six-color 
photometry of the early B stars. The reddening law is one which varies less rapidly 
than 1/d. The observed mean absorptions for a star whose V — I reddening is 1 mag. 
are plotted in Figure 1, together with extrapolations to infinite wave length. An approxi- 
mate fit to the observed values of 6m can be obtained. Assume that the variation of 
5m() is to be represented by an interpolation formula of the form 


5m (n, = (1) 


3 Stebbins, Huffer, and Whitford, Ap. J., 91, 20, 1940. 
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where }, is the effective wave length of the infrared photoelectric measure. Then form 


the difference 
D (n, = 6m — bm (1, d), (2) 


in order to obtain the predicted deviations from the 1/) relation, 


D = — 1/2) — (1 (3) 
| | 
-15 
Am Absorption for 
(\-1) Star with Reddening 
V-I= 100 

-10 — 

Extrapolation 

-05 


7 
Extrapolation 


»* 
20 
+05 on 
+10 4 
| | | | 
25 20 15 10 05 
in po 


Fic. 1.—The observed mean values of the reddening for a star with V — I = 1.0 mag. are plotted. 
The dotted line represents the expected values if the 1/d relation is valid. The solid curve represents the 
values of the reddening if the \~** relation is valid. The zero point of the magnitude scale of reddenings 
should be the value for 1/A = 0.0, i.e., if the 1/ relation were valid, 5m(A) would be +1.69 mag. for the 
point V; if the A~*® extrapolation is adopted, 5m(A) = +2.54 mag. for point V. 


We can then compute D(n, \) for the desired values of and X. The results of such com- 
putations are given in Table 3 and compared with the observed values of D. Auxiliary 
quantities tabulated are the constants cn; dm (I — ©), the extrapolation of our inter- 
polation formula to infinite wave length to obtain the absorption at point J on the as- ; a 
sumption that there is zero nonselective absorption; E(C) is the predicted reddening on es: 
the short-base-line photoelectric scale (1332 B stars) ; A (pg) is the absorption at the pho- 
tographic effective wave length, taken for reddened stars to be 4400 A. 

It can be seen that 6m(I — ~) increases strongly with decreasing nm. In Table 3 a 
minimum value of A(pg)/E(C)) is given in the last line, based on the following assump- 
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tions. Beyond the point I, let us assume that the absorption decreases as 1/); let us also 
assume that there is no nonselective absorption. Then, 


A (pg) = en (1/Xpg— 1/01) + (4) 


A (pg) = 1/27) + 0°69 (5) 


Comparison of the observed and predicted D in Table 3 shows that the index in the 
law of reddening is close to 0.56; a value of 0.75 seems definitely too large. Then the ratio 
of absorption to C; reddening is 14:1, by extrapolating m = 0.56 to infinite wave length; 
the ‘“‘minimum” value, using the 1/\ extrapolation is 9:1. The extrapolation of the ob- 
served absorptions over the V — I region to infinite wave length is sensitive to , par- 
ticularly when m < 1. This can be seen in Figure 1, where the upward curvature ob- 


TABLE 3 
COMPUTED VALUES OF THE DEVIATIONS FROM 1/A 


CompuTep D(n, d) 
CoLor 
D 

n=1.0 0.75 0.62 0.56 0.50 
0.046 0.000 0.024 0.039 0.043 0.050 
0.049 0.000 0.026 0.041 0.046 0.052 
0.031 0.000 0.017 0.028 0.031 0.034 
0.000 0.000 0.000 0.000 0.000 0.000 
0.69 1.05 1.35 1.54 1.77 
4.83 5.24 5.49 5.59 5.71 
py 1.63 1.99 2.30 2.49 2.72 
7.87 | 10.4 12.6 13.9 15.5 
Minimum A(pg)/E(Ci).|......... 7.87 8.54 8.98 9.13 9.35 


tained with m = 0.56 results in an absorption difference between I and © of 1.54 mag. 
This may possibly be an overestimation. If we extrapolate, using the 1/A rela- 
tion, 6m(I — ©) is only 0.69 mag. The small curvature of 0.049 mag. in the V — I 
region results in 0.85 mag. increase in 6m(I — ~), The true value of dm(I — ~) may 
lie between the two estimates; we shall later discuss the possibility that the extrapolation 
might yield less than 0.69 mag. Nonselective absorption, free-free absorption in inter- 
stellar hydrogen, and electron scattering produce some small absorption in the infrared 
and are all neglected in this investigation. The ‘““minimum” value of A(pg)/E(Ci) = 9:1 
may be taken as a relatively conservative one at present. 

These conclusions confirm early spectrophotometric determinations, such as that of 
Greenstein,’ in which an index m = 0.8 was suggested by an upward curvature of energy- 
curves of reddened B stars. The investigation by Greenstein and Henyey,° based on the 
C; scale of photoelectric color excesses, involved completely different assumptions than 
those made in spectrophotometric work and seemed to indicate that the nonselective 
component of the reddening is small; a ratio of about 9:1 was also found for A (pg)/E(Ci). 


4 Ap. J., 87, 151, 1938. 
5 Ap. J., 93, 327, 1941. 
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THE RELATION BETWEEN COLOR-INDEX SCALES 


The new photoelectric material on wide-base-line colors helps to clarify a complex and 
neglected problem in the photometry of reddened stars. There is now little doubt that 
the energy-curves of the reddened B stars deviate from the black-body-curves (assum- 
ing that the unreddened B stars radiate like black bodies). The stars of later spectral 
type may also have small deviations from a black body, caused by spectral lines or by 
different types of variation of opacity with wave length. A single linear reduction from 
one short-base-line color-index scale to another may not be possible. The existence of 
different reduction-curves for space-reddened B stars and for late-type stars (which I 
will call “‘temperature-reddened stars’’) is well shown in a detailed comparison by the 
Greenwich observers® of Greenwich gradients’ with various color systems. In their 
Figures 1, 3, and 4, gradients are compared with three photoelectric scales. Gradients 
are measured from A 4200 to A 6300, while the photoelectric scales have short base lines 
within the region AX 3800-4800. Photoelectric colors (for example, the C; scale), in a cer- 
tain sense, underestimate the reddening; the deviations from black-body radiation now 
found in the B stars will explain, at least qualitatively, the fact that the reddened B stars 
show the C; scale to be relatively compressed in comparison with gradients. Another 
demonstration of this phenomenon is found in a comparison by Sherman and Morgan® 


TABLE 4 
COLOR CLASSES OF REDDENED B STARS 

Star Grad. V-I C1 Becker Bottlinger 
{7 {5 a8 a3 a2 


of color classes of reddened stars, as determined on the spectrophotometric broad base 
line and on three photoelectric systems. Their Table 4 is reproduced, in part, in my 
Table 4, with the addition of color classes determined on the V — I base line and with 
the use of the color classes based on the Greenwich gradients (‘‘Grad.”’). The colors of 
main-sequence, unreddened, later-type stars were used to establish the relation, in each 
color system, between intrinsic color and spectral type; the color class of a reddened star 
was then determined from its observed color. (W. W. Morgan has also demonstrated the 
existence of this effect in several unpublished investigations of the relation between 
spectral type and color.) The Greenwich gradient scale and the V — I scale are both suf- 
ficiently long to give accordant measures of the reddening suffered by these stars; the 
short-base-line colors underestimate the reddening by roughly one-half. The C; scale of 
color excesses is extremely useful in all investigations of space reddening; nevertheless, 
the new V — I colors are the most sensitive means of discovery of reddened stars. 
The V — I scale also permits an accurate calibration of the absorption corresponding 
to a given C; reddening. A similar calibration remains to be made for the International 
scale; it is probable that the relation between V — I and International colors will not 
be the same for late-type stars as for very highly reddened stars. The ratio of absorption 
to International color excess is not accurately known at present and may be best deter- 
mined, not from effective wave lengths, but from direct observations of International 
color excesses of stars with known reddenings on the V — I scale, using highly reddened 
B stars. 


§ M.N., 100, 196, 1940. 7 Tbid., 100, 189, 1940. 8 Ap. J., 89, 516, 1939. 
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The usefulness of broad-base-line colors is clearly demonstrated in Figures 2 and 3. 
In Figure 2 are plotted the observed values of V — I and C; for the North Polar Se- 
quence stars, and for stars of types B3 and earlier (space reddened but with substantially 
similar intrinsic colors). The two groups define clearly separated relations between the 
color systems; a certain nonlinearity is present in both groups. For simplicity, a linear 
relation was fitted to each by least squares (star NPS 2r was omitted from the NPS 


group): 
Temperature-reddened NPS stars , C,=0740+07269 (V-D), (6) 
Space-reddened B stars , C,=0716+07159(V—-I). (7) 


o NPS Stars 7 
+10F— Reddened B Stars os 


+08 
+06 
+04 
+02 


-04 
30 +10 +20 +30 


| 
Fic. 2.—The plotted points represent stars for which C; and V — I are known. The least-squares 
solutions are plotted. 


There is a large difference in slope between the two groups. In equation (7) the scale co- 
efficient is also definitely lower than the value 0.173 derived in Table 2 and also lower 
than any predicted value of E(Ci)/V — Lin Table 3, even with n = 0.56. This difference 
of slope may be affected by the method of selection for the data in Table 2 or by the non- 
linearity of the relation between the scales of colors. It does not seem profitable to refine 
the analysis further; the uncertainty and possible variability with color of the effective 
wave lengths of the C; scale will affect the predicted value of E(C;). The value of (V — I) / 
C; = 1/0.159 = 6.3, in agreement with Stebbins and Whitford,!? will be adopted. The 
result is an increase of the ratio A(pg)/E(C), for n = 0.56, to 15.6, and of the ‘‘mini- 
mum”’ value to 10.3. 

In Figure 3 the Greenwich gradients are plotted against V — I for temperature- 
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reddened stars (types B6 or later) and for early B stars separately. The least-squares 
solutions give 


Temperature-reddened stars , Grad. = 1.30+0.677 (V-I), (8) 
Space-reddened B stars , Grad. = 1.45+ 0.682 (V-—I). (9) 


In these two sets of long-base-line colors there is negligible difference of scale and a small 
zero-point shift. The wave-length range of the gradients, about 0.9 u~', seems sufficient 
to give a measure of the reddening (in color classes) similar to the longer V — I base 
line. This supposition is confirmed by the rather scanty data in Table 4; here the gradi- 


Grad. 

+16-—— @ Reddened B Stars < B85 ° 
© Unreddened Stars > BS 

+14-— — 
+12-— 
_ 
+08-— 
+06-— 
+04-— 
+02-— 
-02-—- 

-250 -200 —150 -100 - 050 000 +050 


V-I Photoelectric 


Fic. 3.—The plotted points represent stars common to the lists of V — I and the Greenwich gradients, 
The least-squares solutions are plotted. 


ents and V — I measures give accordant reddenings, in color classes, which are larger 
than those determined by the short-base-line photoelectric measures. ' 

A similar phenomenon has been discovered observationally in several unpublished 
investigations by Dr. W. W. Morgan at the Yerkes Observatory. Photographs have been 
taken in the galactic plane, using the Tikhoff method; the color-curve of the refracting 
telescope provides a separation of radiation in the different colors. An illustration of a 
plate obtained in such a way is given by Miinch and Terrazas; red and blue light ap- 
pear as separate extrafocal rings. When panchromatic plates are used, the base line for 
colors is approximately 4500-6000 A. Faint reddened B stars have been found with color 
classes in the late K’s; the photoelectric excesses correspond to considerably earlier color 
classes. It is obvious that our knowledge of B stars fainter than apparent magnitude 
+10 is very incomplete; if such stars are in the galactic plane, they will be highly red- 
dened. Consider, now, a B star with a large reddening, for example, HD 169034, a cBO 
star. According to Stebbins and Whitford,' this star has a V — I color of +2.01 mag. 
and a C; color of +0.45 mag. A near-by unreddened B star like HD 204172 has V — I = 
—2.18 mag. and C; = —0.19 mag. The color difference on the V — I scale is +4.19 


° 4p. J., 103, 371, 1946. 
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mag., and on the C; scale +0.64 mag. On the other hand, the star NPS 3r (HD 114282), 
type gK2, has V — I = +2.00 mag. and C; = +0.98 mag. When compared with 
HD 204172 the color difference on the V — I scale is +4.18 mag. and on the C; scale 
+1.17 mag. If we form the differences of colors between the highly reddened B star and 
the temperature-reddened gK2 star, we obtain: 


6(V-I), HD 169034—HD 114282 = +0.01 mag. 
6C:, HD 169034—HD 114282 = — 0.53 mag. 


The B star has nearly the same V — I color as the gK2 star but is — 0.53 mag. bluer on 
the C; scale. We can derive a statistical mean value of this expected color difference on 
the C, scale between a space-reddened and a temperature-reddened star for which 
V — I colors are the same by subtracting equation (6) from equation (7): 


6C, = — (10) 


The particular case of HD 169034, referred to above, would give a predicted value of 
6C; = —0.46 mag., in fair agreement with the observed —0.53 mag. 

Consequently, it may be suggested that new faint reddened stars can be selected as 
follows: A wide-field telescope is used to obtain short-dispersion objective-prism (or -grat- 
ing) spectra in the galactic plane. An infrared emulsion of type I-N (or further to the 
infrared) should be used. The spectra can be broken into three wave-length regions by 
means of suitable filters. Interference-type filters or gelatin filters can be used. Two suc- 
cessive exposures could be made, one with a Wratten 45 filter, which would isolate bands 
near 8500 and 4800 A, and then with a Wratten 35, which gives a band at 4200 A. Pos- 
sible better sets of filters may be developed to obtain three regions simultaneously. The 
three images would possess an extrafocal quality, because of their appreciable wave- 
length range and would permit accurate visual estimates of color classes on a short base 
line in the blue-violet and on a long base line from violet to infrared. Such a classification 
would reveal reddened B stars as objects of large, broad-base-line color and relatively 
small, short-base-line color. For a star like HD 169034, the point at 4800 A would be 
about 0.5 mag. brighter than in normal late-type stars. A spectral classification of sus- 
pected stars might be required to eliminate late-type stars with distorted energy-curves 
from the reddened B stars. A three-ring Tikhoff exposure method or a three-color direct 
photograph could also be used in the discovery of highly reddened stars. 

It is even more striking to compare the long-base-line colors of two stars with the 
same C; color. From equations (6) and (7) we obtain the difference of V — I color of 
space and temperature-reddened stars of the same C}: 


6(V—I) =2.54C,+ 0.49 mag. (11) 


The coefficients in equation (11) are not accurately determined by the present data. 
However, we can compare HD 169034, the cBO star, with the stars of the North Polar 
Sequence. We find that Polaris, HD 8890 (cF7), has C; = +0.49, and V — I = —0.30. 
Consequently, 6(V — I) is +2.31 mag. for approximately equal C;, between the space- 
reddened and the temperature-reddened star. The value given by equation (11) is 
+1.63 mag., which is probably within the accuracy to be expected. (In all these dis- 
cussions no consideration has been taken of statistical effects to be expected when the 
actual regression-curves are replaced by simple linear least-squares relations.) In any 
case, stars with as much space reddening as HD 169034 would show excess brightnesses 
of more than 2 mag. in the infrared, when compared with late-type stars of the same 
C, color. 


es 
] 
: 
] 
r 
( 
( 


~~ 


INTERSTELLAR ABSORPTION 411 


THE NATURE OF THE INTERSTELLAR PARTICLES 


Of the 28 stars now observed by Stebbins and Whitford, only @ Ori, the Trapezium 
group, shows an obviously different type of reddening law. The upward curvature of the 
energy in the blue and violet regions in 6! Ori is much greater than is shown by the mean 
D’s tabulated. This observation agrees with that of Baade and Minkowski.’ The excep- 
tional type of reddening may be connected with the fact that most of the reddening is 
produced close to hot O stars in the dark and bright nebulosity of the Orion nebula. It 
may be suggested that (1) radiation pressure alters the frequency distribution of particle 
sizes" in the dust cloud and thus changes the integrated absorption coefficient ; (2) high- 
temperature ultraviolet radiation changes the optical properties of the material; or (3) 
the high-energy density of radiation in the nebula causes particles of certain ranges of size 
and composition to evaporate. Since many stars have large parts of their reddening pro- 
duced in dense dark nebulae in which the stars are immersed, an observational discyimi- 
nation may be possible. It should be remarked that there exists another peculiarity of the 
absorption in the Orion nebula. The diffuse interstellar band at 4430 A is known to be 
unexpectedly weak in 6! Ori, according to W. W. Morgan; on the other hand, I find 
the band very strong in HD 147889, a B2 star deeply immersed and reddened in a dark 
nebula, in which there is also a strong reflection bright nebula. The presence of the hotter 
stars or the emission lines in the nebula seems to be associated with the difference be- 
tween \ 4430 in @' Ori and HD 147889. (An investigation of \ 4430 in bright nebulae is 
now in progress.) Another suggestive observation is that by Joy,'* of T Tauri variables 
in dense dark nebulosities; here the presence of emission lines suggests possible evapora- 
tion of gaseous material from the interstellar particles close to the stars. 

Excluding 6' Ori, there is a remarkable uniformity in the reddening law; possible 
differences of particle sizes or frequency distributions between thick, dark nebulosities 
and general interstellar space seem to have little effect. The observed reddening law 
covers a very large range of wave lengths; certain types of particles may be searched for 
which would give the observed law. It has been shown" that, with a properly chosen 
frequency distribution of particle sizes, either metals or dielectrics would give almost any 
desired law of interstellar reddening. If, however, we assume that particles of only one 
size are present, the observed reddening law and the ratio A(pg)/E(C,) (determined 
independently) would decide between metallic and dielectric particles. 

Let us assume that particles of one size are present. Let the radius of the particle be 
a, and let us define a = 27a/X. Define the efficiency factor Q(a), such that 


Extinction = Q (a) ra’. (12) 


The quantity Q(a) is an effective cross-section of the particle and is given by-the Mie 
theory for spherical particles; it is a function of the refractive index of the particle V 
and of a. The absorption in magnitudes is proportional to Q; Greenstein" gives values of 
Q(a) for various metals and dielectrics; new computations and a résumé of older data 
are given by van de Hulst."‘ In all cases, Q(a) has a maximum value" near a(N — 1) = 
2; for small a it drops, roughly proportional to a‘ for dielectrics and to a for metals; for 
larger a, Q(a) rises and falls, asymptotically approaching a limit, Q(@) = 2. If we are 
to explain the absorption by a single particle size, we must find a section of these curves 
sufficiently large to cover a change in a by a factor of about 3 (since 1/d varies from 1 to 
2.8 in the photoelectric measures) and such that the slope of the curves is compatible with 
the observed law of interstellar reddening. Particles with a greater than that at which 


10 Ap. J., 86, 123, 1937. 
Harvard Circ., No. 422, 1937. 13 Ap. J., 102, 168, 1945. 
12 4.J., 51, 21, 1944. 14 Optics of Spherical Particles, Thesis, Amsterdam, 1946. 
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the first maximum of Q(a) occurs are excluded, since they would make the stars appear 
bluer rather than redder, and even at point U the interstellar absorption is still observed 
to increase with increasing 1/A. The interpretation of these curves in terms of the inter- 
stellar absorption is further complicated by the fact that most real substances have a re- 
fractive index which varies with wave length; in general, the refractive index increases 
toward the ultraviolet and causes the actual Q(a) to increase more rapidly with decreas- 
ing wave length than is shown in any of the plotted extinction-curves. For example, it is 
known that small metallic particles (a << 1) have an extinction coefficient of the form 


(X) =constant X 1/\ X imaginary part (13) 


where the imaginary part of the complex refractive index yields a coefficient’ to be 
multiplied into 1/A to obtain the final total extinction. For Fe particles the net result is a 
variation near 1/’; for Ni particles, near 1/A*. Consequently, the observed interstellar 
absorption-curve cannot be explained by very small metallic particles and, of course, 
also cannot be explained by very small dielectric particles. The increase of the refractive 
index in both metals and dielectrics toward the ultraviolet is a quite general phenomenon, 
connected with the existence of absorption bands in the far ultraviolet. 

We must therefore determine whether particles of one size, near the maximum of the 
Q(a) curves, can explain the observed interstellar absorption-curve. The total extinction 
per gram of dust particles of this size is about the same" for metals and for dielectrics 
and is about 10*-10°. For purposes of simplicity, let us neglect the variation of refractive 
index with wave length. Assume that the particle radius, a, is such that Q(a) has its 
maximum value somewhere in the near ultraviolet, say at point U (1/A = 2.83). Then 
at point I, (1/A = 0.97), a(I) = 0.342 a(U). Ani investigation of the form of the curves!""!4 
for the behavior of dielectric particles was then made; the results are attended with con- 
siderable uncertainty, since the form of the Q(a) curves is rather poorly defined by the 
available computations near a = 1. If we measure the slope, y, of the curve log Q(a) = 
f(log a), the interstellar reddening will vary as 


6m (A) =constant X (1/A)7 (14) 


in the range of wave length at which the slope y is measured. For a dielectric with 
N = 4, if we take a(U) = 5, the slopes y are 0.59 at point U, 1.2 at point V, 1.9 at point 
B, 2.7 at point R, and 3.1 at point I. This very rapid change of the shape of the dm(A) 
curve is not compatible with the observations. If the index N is greater than #, the varia- 
tion of y is even steeper. In the case of metallic particles, Schalén has shown” that a 
rather more satisfactory fit is obtained; the better fit arises because the limiting value of 
y is unity, if the variation of refractive index with wave length is neglected, when 
a << 1. The recent work of van de Hulst,'* however, indicates that the rate of varia- 
tion of the slope with wave length will be slower if we use a lower index of refraction for 
dielectric particles. Curves for dielectrics with different indices, N, are similar if Q(a) 
is plotted against (V — 1)a. It may be estimated that if V were between 1.1 and 1.2, the 
variation of y with wave length would be compatible with the observations. 

While the interpretation of the absorption as due to dielectric particles of a certain 
size meets with certain difficulties, there is a possibility of distinguishing between metals 
and dielectrics which may eventually require us to assume that the particles are dielec- 
tric. For very long wave lengths, a — 0; for metals the absorption will vary roughly as 
1/d in this range, and for dielectrics as 1/A*. Consequently, the absorption between 
point I and infinite wave length will be smaller for dielectrics than for‘thetals. The re- 


15 Schalén, Uppsala Ann., Vol. 1, No. 2, 1939. 
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sult will be that A(pg)/E(C;) will be smaller for dielectrics. We have shown that if the 
absorption for 1/\ less than 0.97 (i.e., in the far infrared) varies as 1/A, A(pg)/E(Ci) 
will have the “minimum” value of 10.3. It is quite possible that this value is too large 
to be compatible with other direct observational determinations of the ratio of absorp- 
tion to reddening. If this should prove to be true, the absorption beyond point I must 
vary more rapidly than 1/); the dielectric particles would then be a more satisfactory 
source of such absorption. The observed interpolation formula, 1/A°** must break down 
in the fairly near infrared, since it gives a value of A(pg)/(C;) of 15.6, which is almost 
certainly incompatible with direct observations of this ratio. Observations of reddened 
B stars at wave lengths near 20,000 A could be decisive. 

The entire set of arguments can be carried through similarly in case a mixture of 
particles of different sizes are present in space; the effects would be less pronounced, 
and there would be no difficulty fitting the observed law of reddening by the Q(a) curves 
of dielectric particles even with high index. If an upper limit to the particle size present 
in space should exist, there would still remain a possibility of distinguishing between 
metals and dielectrics by the behavior of the absorption-curve in the far infrared. 
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CONTINUOUS EMISSION IN THE ORION NEBULA* 


Jesse L. GREENSTEIN 
Yerkes and McDonald Observatories 
Received August 17, 1946 


ABSTRACT 


Short-dispersion slit spectra of the central part of the Orion nebula have been measured spectrophoto- 
metrically. The energy distribution in the continuum is determined with respect to 6! Ori C; the adopted 
stellar color temperature is based on gradients and photoelectric colors. The nebular continuum has a 
color temperature of 12,000° on both sides of the Balmer limit; the Balmer emission discontinuity js 
1.6 mag. If the energy distribution in the Balmer continuum is interpreted as recombination emission, 
its electron temperature is about 65,000°. This excessively high value can be replaced by a more reason- 
able value near 12,000° if we assume that the optical thickness of the nebula in the Balmer continuum 
is appreciable. Suggested observational tests of this hypothesis include the predicted absence of fila- 
mentary structure in direct photographs taken in ultraviolet light. The complex effects of radiation 
scattered by dust particles in the nebula and of particle reddening are considered. The intensity of the 
hydrogen lines indicates an electron temperature near 9,000°; lines near the series limit reveal a peculiarly 
rapid rise in intensity for high quantum numbers, as they approach the intensity at the series limit. The 
density and size of the nebula are estimated to be sufficient to give appreciable self-absorption. 


THE OBSERVATIONS 


The Orion nebula, NGC 1976, has high surface brightness and a large angular diam- 
eter. The excitation and illumination have been ascribed dominantly to the stars of the 
Trapezium group, 6! Ori, of spectral types O7—-B1. While a large part of the nebular 
radiation is in the emission lines, there exists a strong continuous spectrum. Spectro- 
photometric investigations of the continuum have been made by Greenstein and Henyey! 
and by Barbier,? who measured both the emission discontinuity at the Balmer series limit 
and the color temperature, T., of the continuous emission. Greenstein and Henyey used 
the McDonald nebular spectrograph, and their results refer to the outer parts of the 
nebula, between 3’ and 8’ from the Trapezium; they found an emission discontinuity of ~ 
—0.4 mag. and a gradient of the nebular continuum in the range AA 3920-6060 of +0.85, 
making the nebula slightly bluer than 6' Ori. Barbier? used a conventional slit spectro- 
graph of higher resolution and studied the inner and brighter part of the nebula, about 1’ 
distant from 6! Ori. He finds an emission discontinuity of —1.6 mag. and a gradient of 
+2.14, in the range \\ 4180-4670, making the nebula very much redder than the star. 

The resolution of the McDonald nebular spectrograph is about 30 A; Barbier’s spectra 
seem to resolve about 8 A. The great strength of the emission lines and their confluence 
near the series limit require the highest possible resolution. While Swings and Struve® 
obtained a spectrum with the f{/2 camera of the McDonald Cassegrain spectrograph, the 
long exposures required make photometry difficult. Consequently, I have re-observed 
the central part of the nebula with the two-prism (quartz) f/1 camera at McDonald; the 
exposures are 20 minutes on 103a-O film with a slit width projected on the plate of 0.02 
mm. The dispersion is 200 A/mm at \ 3727, and the observed resolution is about 4A. 
The line 19 is the last hydrogen line resolved on the plate. Microphotometer tracings of 
high and low magnification are shown in Figure 1; the resolution seems to be considerably 
higher than in the tracing reproduced by Barbier.” A reproduction of a typical spectro- 
gram is given in Plate XVII. There was no star brighter than 13 mag. on the slit; a long slit 
was used, running in the east-west direction, and the nebula was kept in a fixed position 


* Contributions from the McDonald Observatory, University of Texas, No. 130. 
1A p. J., 89, 647, 1939. ? Ann. d’ap., 7, 80, 1944. 3 Ap. J., 96, 310, 1942. 
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during the exposure. Spectrophotometric comparison was made with the star 6! Ori C; 
the star was reduced about 5 mag. in brightness by means of a rapidly rotating sector. 
Calibration was provided by a tube photometer with ultraviolet, blue, and green filters. 
One exposure of the nebula was obtained through a rotating sector to reduce the bright- 
ness by 2.5 mag. and permit comparison of the nebular continuum and the strong emission 
lines. Three normal exposures of the nebula and five of the star were also obtained. 
The mean energy distribution in the nebula was determined on three plates; the con- 
tinuous emission can be located on the tracings with a fair degree of certainty between 
d 4000 and \ 4750, as can be seen in the high-magnification tracing in Figure 1; between 
d 3640 and dA 4000 the overlapping of the hydrogen lines is significant. Points to be meas- 
ured were selected to be as free as possible from weak nebular lines listed by Wyse;‘ the 
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Fic. 2.—The differences of magnitude between the Orion nebula and 6 Ori C are plotted against 
1/X in yw. In this figure the measures near the Balmer limit are smoothed, and the emission lines 
not shown. Open circles are uncertain points in the far ultraviolet, where less than three plates were 
measured. The dotted lines represent the adopted relative energy-curves. Please note that the point at 
1/A=2.90 has been plotted 0.10 mag. too low; when plotted correctly, it falls exactly on the dotted line. 


location of the continuum near strong lines could be fairly well established by extra- 
polation from undisturbed regions of the spectrum. Near H8 and N1, N2, the location of 
the continuum is uncertain; and in the far ultraviolet the ozone bands and the weakness 
of the spectrum introduce larger accidental errors. The internal agreement of measures 
in the three nebular spectra and in the five stellar spectra was excellent. Spectrophotom- 
etry was particularly simple because the nebular continuum, in the mean, had the same 
density as the star (differences of magnitude ranged from +0.8 to —0.8) and almost the 
same color. The mean observed differences of magnitude, dm(A), nebula minus star, are 
plotted in Figure 2. 

It is now established that the energy distribution in a normal reddened B star differs 
slightly from that of a black body.®*7 In the single case of 6! Ori, it has been shown®*8 
that the reddening is of an abnormal character. I have plotted the measured differences 


4 Ap. J., 95, 356, 1942. 
5 Stebbins and Whitford, Ap. J., 98, 20, 1943. 7 Greenstein, Ap. J., 104, 403, 1946. 
6 Tbid., 102, 318, 1945. 8 Baade and Minkowski, Ap. J., 86, 123, 1937. 
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of brightness given by Stebbins and Whitford® between the Trapezium group and the 
mean of 6, e, and ¢ Ori, over the range from their point U (1/A = 2.83 w-') to point 
I (1/A = 0.97 uw), in Figure 3. The reddening of 6! Ori is +1.29 mag. on the U — I scale 
of colors. I have also plotted the values of 6m(\) expected if the reddening were like that 
of the normal reddened B stars as well as those expected if the reddening followed the 1/A 
relation (which is the expected form if the radiation is that of a black body). The devia- 
tion from the 1/) relation is large for 6! Ori. The present range of wave length (2.06 < 
1/X < 3.06) is sufficiently small, however, that we can approximate the observed energy 
distribution by a black body. If we use the last three measured points (B, V, U), we can 
approximate by a straight line of slope +0.40, with deviations of only +0.02 mag. The 
relative spectrophotometric gradient, d¢, of 6! Ori with respect to 6, e, and ¢ Oriis +0.37. 


o 6! Ori— BO Stars 
-0.6 -— @ Mean Reddened B Stors 


Reddened 


0.0 


28 26 24 22 20 IW 6 LO Wy 


Fic. 3.—The photoelectric observations of the difference in magnitude between 6 Ori and 4, e, and ¢ 
Ori are shown as open circles. The solid dots represent the mean values observed for normal reddened 
B stars, and the dotted line, the black-body curve (or 1/A reddening). 


The Greenwich color temperature measures’ for 6 and ¢ Ori give their mean gradient as 
—0.16 with respect to the mean AO star; if we adopt 14,000° for the color temperature 
of the mean AO star, we obtain ¢ (6, e Ori) = 0.87, or T, = 23,000°. Then ¢(@! Ori) = 
1.24, or T. = 12,500°. The value of 7. for 6 and ¢ Ori is approximately confirmed by 
Stebbins and Whitford’s measures® of the V — I long-base-line colors. They adopt T. = 
5500° for the standard dG6 stars; then 7, = 20,400° for 4, e, and ¢ Ori, or g = 0.94. On 
their scale, however, the mean AO star has 7, = 11,300°, which may be rather low. A 
slight upward revision of T, for the standard dG6 stars would bring both the AO- and 
the B-star temperatures into agreement with the Greenwich values. The uncertainty of 
T.. for 6 Ori is probably less than 2000°, and the adopted value of 12,500° may possibly be 


considered slightly high.!° 


M.N., 100, 189, 1940. 


10 In all these discussions it is assumed that the intrinsic colors of the four stars of the Trapezium are 
the same; it is also assumed that there is no appreciable difference in the intrinsic colors of 6, e, and ¢ Ori. 
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PROPERTIES OF THE NEBULAR CONTINUUM 


By inspection of Figure 2 we see that the nebula may be slightly redder than 6! Ori C; 
the dotted line at the right represents the observations well and has a slope correspond- 
ing to dg = +0.07, for 1/X < 2.5. For 1/A > 2.75, dg = 0 fits the observations. Conse- 
quently, g (neb) = 1.31 and T, = 11,500° in the blue. In the violet, y (neb) = 1.24 and 
T. = 12,000°. The value of the discontinuity at the series limit, D, is 1.60 in mag. or 0.64 
in logio J, in good agreement with Barbier,’ who finds 0.56. However, Barbier gives 
T. = 6950° in the region AA 4190-4685, and 7, = 10,300° to the violet of the series lim- 
it. The differences between the results of the present investigation and those of Barbier 
are significant in the blue region. His resolving-power is relatively low, apparently, since 
he chooses only three points, at AA 4685, 4567, 4190, as substantially unaffected by emis- 
sion lines and determines the gradient in the nebula from those points. I believe that the 
results of the present investigation, using a larger number of points, determine the energy 
distribution more accurately. Inspection of the present plates shows that line blending 
is not serious for \ > 4000 A. In the ultraviolet region Barbier uses the energy beyond 
the Balmer limit in 6 Ori and a CMa as a standard. He adopts g2 = 1.48, which gives 
T. = 9800° in a CMa; however, the mean AO star in the French color-temperature work 
has ¢g2 = 1.39. The use of this higher 7, for his standard star would bring the color tem- 
perature in the Balmer emission into close agreement with the new value of 12,000°. 
However, his value of 7, = 6950° in the blue is, I believe, definitely too low. 

The bright O9 star, 6? Ori is only 2’ distant from 6! Ori; however, it is bluer by —0.14 
mag. on the C; scale of photoelectric colors.'' We may conclude that 4? Ori is well in front 
of the nebula and also less luminous than 6! Ori. It probably contributes relatively little 
to the excitation and illumination of the nebula. We shall assume that 6” Ori can be neg- 
lected in further discussions of the source of the nebular continuum. It is not probable 
that 6' Ori is completely behind the dark and bright nebulosity that makes up the Orion 
nebula, since the nebulosity seems substantially opaque for radiation of fainter and 
more distant stars. 

The known possible sources of the continuum are: (1) recombination of ionized hydro- 
gen atoms, (2) scattering and reflection by dust particles, (3) free-free emission, and (4) 
electron-scattering. It is improbable that process (2) would show a discontinuity at the 
Balmer-series limit; therefore, for \ < 3646, it is apparent that bound-free transitions in 
hydrogen constitute the main source of the intensity. In the region between A 4000 and 
d 4700 the emission may be of complex origin. Process (2) may be operative, and Paschen 
emission (Pac) similar to the Balmer continuum (Bac) may exist. The ratio of Paschen to 
Balmer emission can be computed theoretically. Since 6' Ori is reddened (and the redden- 
ing seems to be substantial even in the short distance between 6! and 6? Ori), it is possible 
that there is scattering and reflection of starlight by the reddening particles. We shall 
call this diffused continuum C and assume that it is similar to that observed in reflection 
nebulae. The discontinuity in intensity at the Balmer limit will be measured by D, 
where 


I (X}<3646) C+ Bac+2Z Pac 


D = 3646) Pac 


(1) 


We will include the Paschen and higher-series continua and the free-free transitions in 
> Pac. Barbier” has used the observed value of D to find the electron temperature, T., 
of the hydrogen gas; as T, increases, D decreases. If C ¥ 0, T, will always be overesti- 
mated. Electron-scattering also reduces D; however, the optical thickness of the nebula 
due to free electrons is probably small. Barbier gives D as a function of 7, from computa- 


1 Stebbins, Huffer, and Whitford, Ap. J., 91, 20, 1940. 
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tions based on the theoretical work of Cillié.*1* After enlarging his computations, I 
give in Table 1 the theoretical values of D for various values of C/Bac. The value of D 


when C = 0 is given by 


foe] 
2 
~~ etm 
D = logis ‘ (2) 
2 
0° etm 
TABLE 1 
THE BALMER DISCONTINUITY, D 
| 
C/Bac C/Bac 
Te | 
0 0.1 0.2 0.3 0 0.1 0.2 0.3 
5. 2:34 1.02 0.77 0.63 20,000°...| 0.72 0.58 0.49 0.42 
7,500. ..<| 1.68 .96 25,000... .60 49 .42 
10,600.....{ 1.31 .68 30,000. ... 51 41 .36 
15,000....} 0.93 0.70 0.57 0.49 40,000....) 0.38 0.33 0.29 0.26 


where x, is the energy required to ionize hydrogen from the nth level. At the Balmer 
limit, 4v = x2; and we obtain 


1 h? kT 
1 e 
(3) 
h? kT, 


eo 
1 


D= logio 


n=3 


The last term in numerator and denominator of equation (3) measures the free-free emis- 
sion, negligible at low T.; when T, < 7500°, terms with 2 > 4 are also negligible. 

The observed value of D = 0.64 (Fig. 2) gives T, = 22,000° if we assume C = 0. On 
the other hand, if C ¥ 0, the value of T, can be quite low. For example, if C/Bac = 0.2, 
then T, = 12,000°, in agreement with the color temperature. Such a value of C/Bac cor- 
responds to C/E Pac = 2.7; about 70 per cent of the nebular continuum for \ > 3646 
would be diffused stellar radiation rather than Paschen emission. 

Beyond the Balmer limit, both = Pac and C are small compared to Bac. If we neglect 
C for \ < 3646, we get another estimate of T,—one which we should theoretically expect 
to be the most direct and reliable. Essentially as a consequence of the assumed Maxwell- 
ian velocity distribution in the electron gas, the sum of the Balmer and higher continua 
gives an intensity J, (neb): 


I, (neb) dv . (4) 


12 M.N., 92, 820, 1932. 
13 [bid., 96, 771, 1936. 
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Let us compare this intensity with that in a black body at temperature 7,. Then the 


ratio is 
I, (neb) dv (5) 
If we express this ratio in stellar-magnitude differences 6m, nebula minus star, and differ- 
entiate with respect to frequency, we can write equation (5) in a form similar to that 
used in the determination of spectrophotometric gradients. If we use as units 1/\ in yo, 


we obtain 


where ¢ is the usual stellar gradient. Our observations (Fig. 2) show that there is a zero 
slope of the 6m-curve in the region beyond the Balmer limit. Thus, 


14,300 
T. 


The mean wave length is about 0.34 uw, and g = 1.24 for @' Ori; then equation (7) gives 
T, = 65,000°. If we assume that the total uncertainty of the gradient of 6! Ori and of the 
comparison of the nebula and the star gives a possible error of about +0.15 in the quan- 
tities to be inserted in equation (6), we can estimate a rough uncertainty of 7, near 
+20,000°. Barbier obtained 7, = 33,000°, using the same method. While the two de- 
terminations differ, they are essentially in accord in indicating an enormously high elec- 
tron temperature; at these high temperatures a small difference in the assumed color 
temperature of the comparison star makes a very large difference in the electron tempera- 
ture of the nebula. (The same difficulty is encountered in determining accurate color 
temperatures of hot stars.) What seems definite is that the energy in the Bac decreases 
only slowly toward the ultraviolet. If interpreted as a simple recombination spectrum, 
this strong ultraviolet emission gives an electron temperature which can be conservative- 
ly estimated as larger than 40,000°. Such a value seems incompatible with all other de- 
terminations of 7, in planetary or diffuse nebulae and with other methods which can be 
applied to the Orion nebula. 

This surprising phenomenon in the Bac region may be related to the anomalous con- 
tinuous emission observed in the Crab nebula and in planetaries. Let us discard the as- 
sumption made as to the transparency of the nebula to its own radiations and consider 
the effects of self-reversal in the Balmer continuous emission. It should be remembered 
that the 2s’S level of hydrogen is highly metastable and that the hydrogen absorption 
lines are observed in various astronomical sources under conditions of high dilution. We 
consider the effect of self-absorption within the Bac as follows. Let the only source of 
nebular emission in the Bac region be the recombination of ions; neglect the stimulated 
emission because of the high dilution. The ionization is governed by the exciting star and 
by the diffuse nebular radiation in the far ultraviolet. The electron temperature, 7., is 
constant through the nebula; the nebula is considered to be a homogeneous plane-parallel 
slab, with x the depth measured from the front surface and «» the thickness. Let the only 
source of absorption be the bound-free transitions of hydrogen (although the free-free 
transitions are formally included). The emergent intensity of nebular radiation is given 


by 
= dy f cexp(— f ) dx, (8) 


where ¢, is the emission per unit volume, and &, the absorption per unit length. In the 
case of hydrogen we have 
(9) 


k,=C,Ny. (10) 


= 1.086 | +3n], (6) 


=o-—-3). (7) 
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The absorption coefficient involves the number of atoms of hydrogen in the second level, 
N2; and we cannot compute it explicitly if the metastability of the 2s?S level is appreci- 
able. It is also not directly possible to give 7, in terms of the temperature of the exciting 
star without considering the problem of collisional excitation of the forbidden lines of 
other elements. From equation (8) we obtain 


0 
Introducing the optical depth 
Tr=k,x , (12) 
and integrating, we obtain 
/kT —r,(0) 
I,dv= “Gn, edv(1—e (13) 
Equation (13) includes two correct formulae for limiting cases: 
I,dv=N;N dy, as 7, (0) (14) 
I,dy= dy as T, (0) (15) 


In equation (14) we have the energy distribution and intensity predicted by the recom- 
bination formula (9) in a thin layer; and in equation (15) we have essentially the black- 
body formula in the absence of stimulated emission. The form of equation (15) is to be 
expected, since in the limit an opaque mass radiates like a black body at its boundary 
temperature, which is 7,. If we use equation (15) to compute the electron temperature 
from the color temperature of the observed Bac, we of course obtain 7, = T-, instead 
of the value given by equation (7). If7,(0) = 3 at the Balmer limit, the energy distribu- 
tion given by equation (13) differs by only 10 per cent from the black- body distribution 
in equation (15), within the short range of wave lengths observed. In the far ultraviolet 
the energy in Bac should eventually correspond to a lower 7, than near the Balmer limit. 
However, if 7,(0) is larger than 3, the color and electron temperatures will agree well. 

If strong self- -absorption exists within the Bac, it is not necessarily true that similar 
self-absorption exists within the Paschen continuum, since there is no metastable level 
above 2s*S. Consequently, the ratio of continuous absorption coefficients will have a 
strong discontinuity near the Balmer limit. The emission lines may also show self-re- 
versal, depending on their Doppler width; the theory of the Balmer decrement might re- 
quire some revision. The discontinuity in absorption coefficient near the series limit would 
be partly masked by the increase of the mean line-absorption coefficient.'* A general de- 
pression in the ultraviolet energy distribution of stars seen through the nebula might be 
expected. (Such a depression is not visible in the spectra of the Trapezium stars.) An- 
other observational test involves the filamentary structure of emission nebulae. In photo- 
graphs taken at A < 3646, sharp features, such as bright filaments, would be weakened 
and diffuse in appearance if they are behind most of the nebular gaseous absorption. On 
the other hand, photographs taken by the light of forbidden lines (in which no line self- 
absorption exists) at \ > 3646 could show sharp filamentary structure, diffused only by 
the dust-particle absorption. Remarkable photographs by Baade" of the Crab nebula 
show in general that the filamentary structure is very pronounced in the [NV u] lines 
and much weaker in the blue and violet where continuous emission dominates; (the Crab 
nebula is not directly comparable with the Orion nebula because of the absence of hydro- 


14 By continuity, the large absorption coefficient in the Bac requires that the line-absorption coefficient 
in the higher series members also be large. 


16 Ap. J., 96, 188, 1942. 
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gen lines in the former). It would be interesting to determine whether photographs of the 
Orion nebula, taken at \ < 3646, reveal a more amorphous structure than those taken 
in the blue, green, or red. In many nebulae the emission lines are influenced by general 
expansion. If the line He 1, \ 3188 (23S—4*P°) could be observed with sufficiently high 
dispersion to be seen as double, then the longer-wave component should be relatively 
weakened in comparison to the longer-wave component of \ 3888 (2°S—3*P°). 

The hypothesis of strong self-absorption within the Bac is novel and untested; it is ap- 
parently required by the anomalously high 7, given by the energy distribution in the 
Bac, using equation (14). A value of JT, much higher than 10,000° is probably excluded 
by our knowledge of the spectra of planetary nebulae. In spite of their higher excitation 
and ionization and their hotter exciting stars, electron temperatures'® seldom exceed 
10,000.° Aller'’ has pointed out that the low intensity of \ 4363, the more highly excited 
forbidden line of [O m1], corresponds to 7, = 10,500°. In spectra of stars seen partly 
through the Orion nebula, O. C. Wilson'* found the absorption line \ 3888 of He 1, direct 
evidence for the population of a highly excited metastable state. We shall temporarily 
adopt 7, = T. = 12,000° for the Bac emission; then we obtain a consistent group of 
measures of T, in the nebula. The Bac yields 12,000°; the small Balmer discontinuity 
can be explained by the existence of some reflection spectrum; the forbidden-line emis- 
sion spectrum is about that expected at this 7,. If no Balmer continuous absorption ex- 
ists, the energy distribution in Bac gives T, between 65,000° and 33,000° (Barbier’s 
value); the Balmer discontinuity, without any reflection spectrum, gives 22,000°; the 
forbidden line spectrum, 10,500°. 


THE REFLECTION CONTINUUM AND THE EFFECT OF REDDENING 


An observational distinction between the reflection continuum and the recombination 
hydrogen continuum in the blue is difficult. Polarization measurements in the spectrum 
may be possible, but polarization is not always present in reflection nebulae. The Paschen 
discontinuity is relatively small and badly distorted by terrestrial absorption bands. (If 
we define a quantity like D in eq. [1] for the Paschen limit, it is 0.48 for 7, = 10,000° and 
0.20 for T, = 25,000°.) A reflection continuum may be expected from the general theory 
of reflection nebulae, unless the dust particles in the Orion nebula are very different from 
those near later-type stars in the Orion region. A rough value of the intensity of the con- 
tinuum, C, may be computed. Consider the illumination by the Trapezium stars at a 
distance r of a part of the nebula lying in the celestial sphere; let R be the distance of the 
stars from the earth and m their apparent magnitude as seen at the earth. An observer 
at r will see the Trapezium as of magnitude m’, where 


2 
= m—t+ 2.5 () (16) 


The absorption in the line of sight to the earth is ¢;. Express the angular separation of the 
point in the nebula from the star, in radians, as a; then 


m'’=m—t+5 logioa . (17) 


If interstellar dust exists at that point, with total optical thickness &, diffusing starlight 
uniformly in all directions with albedo y, the first-order diffused intensity is, per unit 
solid angle, 


(18) 


16 Menzel, Aller, and Hebb, Ap. J., 93, 230, 1941. 
17 Pub. A.S.P., 58, 165, 1946. 18 Pub. A.A.S., 9, 274, 1939. 


| 
| 
E's 
og 


ORION NEBULA 423 


in units of the incident intensity. Convert to surface brightness, M(m), of the nebula in 
magnitudes per square degree at a distance a’ in minutes of arc: 


M(n) =m+t—t+5 logwa’— 2.5 logis (yt) — 6.1. (19) 


This formula is similar to Hubble’s relation. Assume that the optical depth is about uni- 
ty, 4 = tf = 1, and that y = 0.5. Then at one minute of arc from the Trapezium, 
M(n) = +0.05 mag. per square degree in the visual range of the spectrum (since m = 5.4 
for the Trapezium refers to the visual range). This is an extremely high surface bright- 
ness, greater than that of most planetary nebulae, but evidence is not available to prove 
it sufficient to account for the observed continuum. 

There exists some reddening material in front of 4 Ori; if this material produces the 
continuous spectrum, C, and is also located in the region where hydrogen emission occurs, 
the problem of color effects in the nebula is very complex. Certain simple models de- 
scribing the effect of the dust on the hydrogen emission may be considered; we wish to 
consider the effect of the reddening on the color temperature of the continuum. 

1. If the dark material lies completely in front of both @ Ori and the nebular emission, 
the nebular emission is reddened as much as the light of 6! Ori. Then we correct ¢ (6! Ori) 
to give a color temperature equal to its effective temperature. The value of T, for 
 < 3646 is then about 20,000°. There is no expected continuum, C. 

2. Let the hydrogen and dust particles be well mixed throughout the nebula, and let 
the nebular particles have zero albedo. There will be no expected continuum, C, caused by 
diffused starlight. Model reflection nebulae have been investigated by Henyey and 
Greenstein,!* but in these the fluorescence process was neglected. Let the rate of emission 
of Bac be proportional to € per unit volume and let the absorption coefficient of light by 
the dust be k’ (omitting the effect of hydrogen absorption within Bac). The equation of 
transfer can be derived as follows: 


dIl=—k'Idst+eds. (20) 
The solution is of the form 
€ 
Ten= f erdr. (21) 


The optical thickness of the dust is 7}. From the reddening-curve of 6! Ori in Figure 3 
it is apparent that 7 is very large; therefore let r; + ~. The surface brightness of the 
nebula, then, is 


(22) 
Consider two wave lengths, with A: < Au, 
(23) 
content (24) 
€2 ky 


The reddening of the Bac by dust absorption is measured by the color excess, E, as fol- 
lows: 


E= —2.5 (1ogio “3 — =), (25) 
€2 €1 
ks 
E= 2.5 logio - (26) 
ky 


19 4p. J., 88, 580, 1938. 
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If k’ varied as i/X", then 


E=2.5n logio (27) 


While equation (27) could be used with normally reddened B stars, in the case of 4! Ori 
it seems best to use the observed reddening law. From Figure 3,  < 1. The absorption 
at point J (1/ = 0.97) is unknown but may be estimated at 1 mag. Then at A, = 4000 A, 
ky = 2.18; and at \; = 4800 A, &; = 1.99. Inserting these values in equation (26), we 
find the reddening of the nebular emission to be E = 0.10 mag. If we correct the nebular 
continuum in the region 4000 < \ < 4800 for this reddening, we obtain a correction, 
6g = —0.17. Then the true nebular continuum has g(neb) = 1.07, or T, = 16,000° in 
the blue. In the ultraviolet the nature of the reddening-curve is not known accurately, 
but for \ < 3646 we may estimate JT, = 14,500°. 

3. In models (1) and (2) we have neglected the diffused stellar radiation completely. 
Since the albedo of the dust may be high, one limiting model is that with y = 1. The 
particles scatter radiation without absorption. Since energy is conserved, all radiation 
emitted by the nebula escapes at the same rate as it is emitted. The observed and initial 
intensity distributions of the nebular radiation are the same, and no correction is re- 
quired to the color temperatures, except for the reddening in front of the nebula. In model 
(3), T, is about 12,000° in the Bac region, if we assume, as before, that the nebula is opti- 
cally thick in hydrogen absorption beyond the series limit. In case y = 1, the diffused 
stellar radiation in the blue would be large and the Pac emission relatively unimportant. 
The energy distribution in C then could be studied according to the theory of reflection 
nebulae.'* However, since the variation of absorption with wave length is abnormally 
slow in @' Ori, it is improbable that the color of the nebula and that of the star would 
differ appreciably.”° 

In summary, we may estimate that the dust absorption may change the color tempera- 
ture of the nebular emission and that 12,000° < T. < 20,000°. Part of the continuum for 
\ > 4000 probably originates by diffusion by small particles in the nebula. The intensity 
of this reflection continuum lies between 10 and 20 per cent of that in the Balmer con- 
tinuous emission. It is suggestive that the Balmer jump observed in the outer part of the 


Orion nebula! was smaller than that measured in the present investigation. It is possible . 


that the reflection continuum, C, is relatively more important in the outer parts of the 
nebula, thus reducing the observed D. A relative decrease of C in the inner part of the 
nebula could be caused by radiation pressure, which reduces the density of interstellar 
dust near 6! Ori, or by the alteration of the optical properties of the dust particles in a 
region of relatively high radiation density. 


THE LINE SPECTRUM 


One spectrogram was obtained with a rotating sector to reduce the intensity of the 
emission lines by a factor of 10. The exposure time was the same as for the normal nebular 
exposures. The photometry of the lines, unfortunately, is not accurate, since a narrow 
slit was used. Further, 1 is in a region where the plate sensitivity is decreasing rapidly. 
The accuracy of line intensities for H8, Hy, and Hé may be rather low. The lines H9, 
#10, and H11 are measurable on all the nebular spectra. The other hydrogen lines have 
blending features, mainly He1 and [Ne ut]. The intensity of the lines was expressed in 
units of the energy within one wave-number of the stellar spectrum; the energy at \ 3646 
(in the Bac) was measured also with respect to the star. Since we know the color tempera- 
ture of 6' Ori, we can determine the ratio of the intensity in a line to that in one wave- 
number at the series limit. The lines of He 1 at the series limit do not offer any serious 
problem. The results are given in Table 2. The assumptions involved in the computation 


20 Tbid., eqs. 12 and 13’, and Figs. 1 and 4. 
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of 7, are that the lines are produced by recombination without self-reversal; the theory 
is that of Cillié!8 and the more complete theoretical investigation by Baker and Men- 
zel.?1 On assumption B of Baker and Menzel, it is found that the weighted mean value of 
T. indicated by the observed line intensities is about 6500°; if we use the theory of Cillié, 
T. = 9000°. The effect of reddening by dust within the nebula is to make the true value 
of I,,/Bac smaller than the observed, thus lowering the value of T, deduced from Table 2. 
It is quite possible that the assumption of simple recombination in the production of the 
Balmer lines is not correct; the rather low value of 7, here derived is not of great weight. 


TABLE 2 
LINE INTENSITIES COMPARED TO CONTINUOUS EMISSION 
Line Logio Jn/Bac Line Logio In/Bac 
25 
20-— 
Emission Intensity 
Orion Nebula = 
os Balmer 
| | | 


3687 3704 3712 A 
He | Hi6,Hel HIS 


3600 3613 
Hel 


Fic. 4.—Intensity in the hydrogen emission near the series limit. The dotted line at the lower right 
is the unit of intensity, the extrapolated value of the emission for 4000 < \ < 4700. The dotted line at 
the upper left is the ultimate value of the Balmer emission at \ < 3646. The dots are the observed in- 
tensities. Note that He1 lines blend with H 16 and disturb the region near the Balmer limit, The bar 
labeled Res. shows the estimated resolving-power. 


THE BALMER DISCONTINUITY 


The intensity of the hydrogen lines may be correlated with the peculiar nature 
of the Balmer emission near the series limit. The microphotometer tracings in Fig- 
ures 1 and 2 and the detailed emission intensities in Figure 4 show clearly the presence 
of a new phenomenon. In Figure 4, I plot the intensities of the lines H15—H19, the inten- 
sity in the region of overlapping lines, and the series limit itself. The underlying emission 
of = Pac and C is subtracted, and the intensity plotted is that of the hydrogen emission 
lines and the Balmer continuum. The unit of intensity is the intensity in 2 Pac and C, 
extrapolated to \ 3646. The microphotometer tracings were measured every millimeter 
(corresponding to every 2 A). Individual lines disappear just beyond H19, the point at 
which the expected resolution, 4 A, becomes insufficient to separate the lines. From H19 


1 4 p. J., 88, 52, 1938. 
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on, over a range of 40 A to the series limit, the intensities rise, becoming constant at the 
series limit (except for the He 1 lines). According to the correspondence principle, in ther- 
modynamic equilibrium the mean intensity of the emission lines lying within a frequency 
range dv near the series limit should be continuous with the intensity in the continuum, 
per dv. This is certainly not true for lines as high as 19 in Figure 4. The line H16 is 
blended with He 1; taking the average intensity in H15, H17, H18, and H19, we find a 
mean intensity of J = 1.7. At \ < 3646, however, the Bac clearly gives J = 3.4, as does 
the region far to the ultraviolet (see Fig. 2). The appearance of the very rapid increase 
in intensity of high series members is unusual. In chromospheric and prominence spectra” 
there is a clear continuity of the higher series members with the continuum. Other nebu- 
lae have not yet been studied in detail. The ratio of mean line emission to Bac for the 
lines H9, H10, and H11 is 0.37 and for H15-H19, 0.50. Obviously, for large principal 
quantum numbers the ratio is 1. 

We must consider in more detail the behavior of the expected line emission near the 
series limit. Menzel and Aller®* give the intensity in the continuum at the head of the 


Balmer series as 
Tindy = « (28) 


In equation (28) Kh is a numerical constant, and giz is the Gaunt factor (0.88 at the Bal- 
mer limit). For a hydrogen line, they give 


Ins = ENN (29) 


where 6b, is a dimensionless factor representing the departure from thermodynamic 
equilibrium (in which 6 = 1). We obtain the ratio of line to continuum as follows: 


Ine Onn22R AR/n*kT 


The interval in frequency between successive Balmer lines approaches 2R/n* as n be- 
comes large; also as n becomes large, for reasonable values of 7., the exponent of e in 
equation (30) approaches zero. Call the mean intensity when a given line is spread out 
over the interval between successive lines, Ing. We find that equation (30) becomes 


§n2 


Tz 


The ratio of the Gaunt factors, g, becomes of the order of unity, and the variation of 
mean line intensity with 7 is merely that of 4,. In the limit, asm — ~, 6, — 1. If we ex- 
amine the values of 6, tabulated by Baker and Menzel”! in their computation of the in- 
tensities of the Balmer lines, we find that 5, lies well below unity, even when n = 30, 
unless T, is very high. Adopting their model B, the nebula optically thick in the Lyman 
region, a value of };7 near 0.5 (the observed value for the lines H15—H19) requires that 
T. be less than 10,000°. An electron temperature as high as 28,000°, suggested as the mean 
value by Barbier,” gives 6:7 = 0.75 and seems to be excluded. The accuracy of the ob- 
servational determination of 5,7 is actually quite high; only small differences of intensity 
and wave length are involved, and the effects of interstellar absorption are negligible. 
Therefore, we may view the value of b,, for high values of m, as perhaps the most reli- 
able method of determining T.. We obtain T, = 9000°. The correction of the observed in- 

22 See, e.g., microphotometer tracings by Dimitroff and Menzel in Harvard Ann., 105, 99, 1937, Pl. IT. 
In the chromosphere, 7, is lower than in the Orion nebula; if the phenomenon involved only T,, it should 
be even more pronounced in the chromosphere. However, it is the deviations from thermodynamic equi- 
librium that must be responsible for the difference. 


as n—>o , (31) 


23 Ap. J., 93, 195, 1941. 
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tensity distribution just to the red of the series limit for the effect of instrumental re- 
solving-power would provide a detailed comparison of the observations with the theory 
of the asymptotic value of b, for very large n; the present dispersion is too low for such 
an attempt. A more rigorous use of equation (30) confirms the general argument based 
on equation (31). 

Near the series limit we can also neglect the effect of self-absorption, because of the 
continuity of the absorption coefficients in the lines and in the continuum. However, it 
would be an interesting theoretical problem to obtain d,, the deviation from thermody- 
namic equilibrium, in the presence of partial self-absorption. If the nebula is opaque at 
certain wave lengths, effects due to the stratification of the elements may also be of im- 
portance. 


IONIZATION, EXCITATION, AND ABUNDANCE OF HYDROGEN 


The absolute intensities of the lines and continua can be used to obtain an estimate of 
the number of atoms of hydrogen in the nebula. In this investigation no new data on ab- 
solute surface brightnesses were obtained. Some approximate quantitative considera- 
tions may be of interest. 

In the nebula the electron density is high and the energy density of radiation also high, 
compared to conditions in interstellar space. If hydrogen, in relatively high density, ex- 
ists near the Trapezium group, there will be an inner region where the hydrogen is mainly 
ionized and an outer region in which it becomes mainly neutral. Bengt Strémgren** has 
considered the problem of the extinction of the ultraviolet ionizing radiation of the O 
stars in a hydrogen cloud in space and has computed a distance, so, at which the transi- 
tion from ionized to neutral hydrogen occurs. It is tempting to consider that the bright 
nebulosity has a size roughly measured by this So, since at greater distances the hydrogen 
is mainly neutral and the recombination process will not operate. I shall assume that the 
stellar radiation is that from a black body at T = 35,000° (color temperature) ; the elec- 
tron temperature is taken as 10,000°, and the radius of the star as four times the solar 
radius. Then 

So = 20N-*/8 parsecs , (32) 


where N is the total number of hydrogen atoms per cm! (ionized NV; plus neutral No). In 
planetary nebulae, V, = NV; = N is of the order of 10*. The value of so then is 0.043 psc, 
which seems too small. The Orion nebula is at a distance of about 400 psc, and the inner 
bright region is about 4’ in radius, or 0.47 psc. Such a value is obtained if V = 300 per 
cm’, The approximation involved is the neglect of the problem of radiative transfer of 
the Lyman continuum and lines in Strémgren’s treatment. 

At such low electron densities in most of the bright nebula and in the surrounding 
neutral hydrogen region, there will be a strong tendency for the level 2sS of hydrogen to 
become strongly overpopulated with respect to the 2p?P° levels. The metastability of the 
2s level is nearly complete. Breit and Teller” find that the mixing of the 2s and the 2p 
states by random electric fields of passing electrons gives a transition probability 2s — 2p 
of about 5 X 107? sec~", at T, = 10,000° and with NV, near 10*. The simultaneous emis- 
sion of two photons results in transitions from 2s to 2p with a probability of about 
8 sec. Since this lifetime is similar to that of the upper levels of many observed for- 
bidden lines, the de-excitation of 2s by ionization or collision is probably also rare. Con- 
sequently, we may expect that the Boltzmann distribution will give the number of excit- 
ed atoms as follows: 


N (2s) _ N (2s) 


= = —5040x,/T 33 
(is) Ne 10 (33) 


4 Ap. J., 89, 526, 1939, 25 Ap. J., 91, 215, 1940. 
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Then N(2s) is about 10-4 No, for T, = 12,000°. Now if we compute the absorption co- 
efficient per atom in state 2s, at the Balmer series limit, we find that k, (atom) is near 
10~'’. If we have a column of neutral hydrogen, xo cm thick, the optical depth will be 


Try (0) =%No10-410-" . (34) 


If we wish 7,(0) to be unity, x»V» = 107'. Now, outside the ionized region we can take 
No = 10°, so x) = 10'8, or about 0.33 psc. To obtain opacity within the Balmer continu- 
um we thus require a region of neutral hydrogen of the same density and about the same 
size as the observed bright nebula. 

A rough numerical check is obtained from Ambarzumian’s evaluation” of the number 
of excited atoms in the Orion nebula from the observed surface brightness of the nebula 
in H6 and Hy. He found that there were about 10‘ atoms in each state (n = 4 and 5) per 
cm? column through the nebula (assuming no self-absorption within the lines). Since 
there is no appreciable effect of metastability for these levels, their population is re- 
duced from the Boltzmann value by the dilution W. (We neglect a possible more com- 
plete discussion based on the work of Baker and Menzel*!). Then 


0 


If we could estimate W exactly, equation (35) would give the number of neutral atoms. 
However, the exciting radiation is variable through the nebula, and the excitation cannot 
be described by the Boltzmann formula. For a rough approximation we take the value of 
W = 10~" (at the edge of the ionized region) and use the electron temperature 7, = 
12,000°. Then Ambarzumian’s observed populations yield xoNo = 10” and the value of 
7,(0) is 10. While preliminary, these attempts to determine the opacity in the Balmer 
continuum suggest an appreciable order of magnitude. The total mass of hydrogen in 
the nebula is about ten times the solar mass. 


ADDENDUM 


In private correspondence Dr. Donald H. Menzel of the Harvard Observatory has pointed 
out some theoretical difficulties in the derivation of equation (15). I am deeply indebted to him 
for discussion of some of these points. I attempted to show that T. becomes equal to 7. for a 
nebula optically thick in the Balmer continuum. In such a nebula it would be necessary to con- 
sider theoretically in detail the energy transfer in the Lyman continuum as well and derive the 
variation of 7, with distance from the exciting star. Further, I assumed that the hydrogen is 
partly ionized and partly neutral in the same volume of space. This volume must be large enough 
to contain a sufficient number of ionized (recombining and emitting) atoms, as well as neutral 
(absorbing) atoms to produce opacity. In the discussion of Strémgren’s theory in the last sec- 
tion, we see that the ionized and neutral regions may be mainly separated in space. It is possible 
to compute the energy distribution in an emitting nebula surrounded by an absorbing hydrogen 
cloud. The emergent energy is then quite different from that of a black body. Instead of equation 


(15) we find that 
= exp [—hv/kT, — 7,(0)(v0/r)'] , (36) 


where v is the frequency at the Balmer limit. The last term depends on the variation of the 
hydrogen absorption with frequency. The exciting stars, if seen through this neutral hydrogen 
absorbing material, would show a Balmer absorption discontinuity. If 7,(0) = 1, T. = T. at 
the series limit as in our first model, but there would be appreciable change of color tempera- 
ture with wave length. If 7,(0) is large, 7. would be high, and the series limit would appear as 
an absorption discontinuity. If such more complex phenomena should be observed in nebulae, 
we may have to consider this second model. 


26 Zs. f. Ap., 6, 107, 1933. 
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)- The theoretical situation is incompletely investigated at present. It would be profitable to 
ir consider the emission intensities to be expected in a nebula in which the second level of hydro- 
gen is treated separately for the 2s and the 2p states. Even in an optically thin nebula, the values 
of 6, might be changed. Dr. Menzel further points out that if the nebula is opaque, 5, might 


) approach unity even for small values of x. For example, in the chromosphere and probably in 
Be stars the observations indicate that } is close to unity. The total amount of hydrogen per 

e unit area in a Be star may well be of the same order as in the Orion nebula. The dilution is very 

1- much greater in the nebula, however. 

e The full explanation of the observations must wait on a more complete theory of the recom- 


bination emission of hydrogen. If, in fact, high opacity does require 6, near unity, the observa- 
tion that b;; = 0.5 may indicate that the nebula is not opaque. Unfortunately, the low observed 


a value of b,; indicates a low T., while the energy distribution in the Balmer continuum requires 
a a high 7,. In such a new theory it would be desirable to re-evaluate also the expected Balmer 
I emission discontinuity, which may be affected by the self-absorption from the 2s level. In a 
e planetary nebula; where the complications arising from the reflection continuum do not exist, 
p- the theoretical predictions of the Balmer discontinuity, the values of 5,, and the energy distribu- 


1- tions in Pac and Bac can all be compared with observation. 
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ON THE CONTINUOUS ABSORPTION COEFFICIENT OF THE 
NEGATIVE HYDROGEN ION. III 


S. CHANDRASEKHAR AND FRANCES HERMAN BREEN 
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Received August 14, 1946 


ABSTRACT 


In this paper the contribution to the continuous absorption coefficient of the negative hydrogen ion 
by the free-free transitions is evaluated in terms of the matrix elements of the acceleration in the Hartree 


field of a hydrogen atom. The new coefficients are larger than the earlier determinations by factors ex- 


ceeding 10 over the entire range of wave lengths of astrophysical interest. 

Tables of the continuous absorption coefficient of H~, including both the bound-free and the free-free 
transitions for various temperatures (2520° < T < 10,080°) and wave lengths (A > 4000 A) are also 
provided. It is further shown that the new coefficients aré sufficient to account for the solar continu .us 
spectrum from A 4000 A to A 25,000 A. 


1. Introduction.—In two earlier papers! the continuous absorption coefficient of the 
negative hydrogen ion has been considered, and its cross-sections for the radiative proc- 
esses leading to its ionization have been determined with some degree of definitiveness. 
And recent astrophysical discussions? relating to the continuous spectrum of the sun 
have shown that the cross-sections which were derived for these ‘‘bound-free”’ transitions 
of H- are adequate to account for the continuous absorption in the solar atmosphere be- 
tween 4000 A and d 10,000 A. Beyond 10,000 A, however, there appears to be an 
additional source of absorption, which it would be natural to suppose is due to the radia- 
tive transitions of free electrons in the field of neutral hydrogen atoms.’ But the exist- 
ing evaluations‘ of these “‘free-free” transitions make them insufficient to account for 
the observed amount of absorption beyond } 10,000 A by factors exceeding 10. Indeed, 
on the strength of this discrepancy, the existence of a hitherto “unknown source of ab- 
sorption” has been concluded.® However, as similar conclusions relating to the bound- 
free transitions have proved premature in the past, we have examined the earlier evalua- 
tions of the free-free transitions of H~ and have found, as we shall presently explain, that 
there are ample grounds for mistrusting them, even as to giving the correct orders of mag- 
nitude. We have, accordingly, made some further calculations to obtain estimates of the 
free-free transitions, on which at least some reliance could be placed. It is the object of 
this paper to present the results of such calculations and to show that these newly de- 
termined cross-sections for the free-free transitions, together with the cross-sections for 
the bound-free transitions given in Paper II, are sufficient to account for the continuous 
absorption in the solar and in the stellar atmospheres of neighboring spectral types in a 
manner which dispels any remaining belief in an “unknown source of absorption.” 

2. The inadequacy of the Born approximation for the evaluation of the free-free transi- 
tions of H~.—The essential approximation which underlies all the earlier evaluations‘ is 
that of Born’s for describing the motion of an electron in the field of a hydrogen atom. 


1. Chandrasekhar, Ap. J., 102, 223, 395, 1945. These papers will be referred to as Papers I and II, 
respectively. 
2G. Miinch, Ap. J., 102, 385, 1945; D. Chalonge and V. Kourganoff, Ann. d’ ap. (in press). 


Pa gee astrophysical importance of this process was first suggested by A. Pannekoek, M.N., 
91, 162, . 


‘L. Nedelsky, Phys. Rev., 42, 641, 1932; D. H. Menzel and C. L. Pekeris, M.N., 96, 77, 1935; J. A. 
Wheeler and R. Wildt, Ap. J., 95, 281, 1942. 


5 See particularly the discussion of Chalonge and Kourganoff (0. cit.). 
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This is apparent, for example, from the agreement of the calculations of Menzel and 
Pekeris with those of Wheeler and Wildt, in which the Born approximation is explicitly 
made. However, on consideration it appears that for the range of energies which occur 
in stellar atmospheres the Born approximation for incident s-electrons must be a very 
bad one; for, as the mean energy of the electrons in a Maxwellian distribution at tem- 


perature T is 
1k = 0.0243 (1) 


when expressed in atomic units, it is evident that it is only electrons with k? < 0.05 that 
will principally contribute to the continuous absorption. And for s-electrons with as small 
energies as these, the Born approximation must fail. This is indeed well known from the 
work of P. M. Morse and W. P. Allis.* But it may be useful to illustrate this failure in a 


- manner which will emphasize the magnitude of the errors to which the Born approxi- 


mation may lead in the cross-sections for radiative transitions. For this purpose, we have 
compared in Figure 1 the radial s-waves 


xo (Born) =sin kr (2) 


on the Born approximation with the properly normalized s-spherical waves in the 
Hartree field of a hydrogen atom, which the authors have recently tabulated.’ It is seen 
that for r < 2 the Hartree waves have amplitudes which are larger than xo (Born) by 
factors which, on the average, exceed 3. However, similar comparisons between the 
p-waves 


sin kr 
—coskr (3) 


x: (Born) = 


on the Born approximation and the corresponding Hartree waves show that equation (3) 
provides a satisfactory approximation for the energies in which we are interested. Re- 
membering that, in the evaluation of the matrix elements of the acceleration, the wave 
functions for r < 2 are all that matters (cf. Fig. 2, in which we have plotted the accelera- 
tion 7 in the static field of a hydrogen atom) and, further, that the contribution to the 
absorption coefficient arising from the s— p and p — transitions must far outweigh 
all the others, it is apparent that in the Hartree approximation we shall obtain cross- 
sections which will be larger by ten or more times the values obtained with the Born ap- 
proximation. In other words, we may expect a correct evaluation of the free-free transi- 
tions of H~ to bring about an agreement between physical theory and the demands of 
astrophysical data without the need of postulating a still unknown source of continu- 
ous absorption. 

3. Formula for the evaluation of the cross-sections for the free-free transitions of H~ on 
the Hartree approximation.—Our remarks in the preceding section have shown the in- 
adequacy of the Born approximation for the evaluation of the free-free transitions of H-. 
The use of the Hartree approximation suggests itself as the next best, though the effects 


- of exchange and polarization may very well be appreciable for the very slow s-electrons 


in which we are primarily interested. But it may be hoped that these latter effects will 
not, at any rate, affect the orders of magnitude of the derived quantities! In any case, 
to improve on the Hartree approximation would require an amount of numerical work 
which will be several fold; and the task is immense even as it is.* These considerations, to- 


6 Phys. Rev., 44, 269, 1933. 
7S. Chandrasekhar and F. H. Breen, Ap. J., 103, 41, 1946. 


8 For example, the present work has required the numerical integration of 63 radial functions and the 
evaluation of 523 infinite integrals, not to mention the computation of numerous auxiliary functions and 
tables. (All this work was done with a Marchant.) 
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Fic. 1.—A comparison of the radial wave functions of s-electrons in the field of a hydrogen atom on 
the Born (B) and the Hartree (H) approximations for various energies of astrophysical interest. (The é 
abscissa measures the distance from the center in units of the Bohr radius). 
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Fic. 2.—The acceleration # = [r—* + 2(1 + r—')]e-*” of an electron in the Hartree field of a hydrogen 
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gether with the urgency of the astrophysical needs, have prompted us to undertake in 
some detail the evaluation of the cross-sections for the free-free transitions of H—, using 
the s- and the p-spherical waves in the Hartree field of the hydrogen atom which we 
have tabulated in an earlier paper.’ 

Now, in the Hartree approximation a wave function representing a plane wave at 
infinity is given by 


i! (21+ 1) 
v= Py (u) xu (15 (4) 
where the radial functions x:(r; &) are solutions of the equation 


which have unit amplitude at infinity. (In the foregoing equations we have adopted the 
atomic system of units). For wave functions of the form (4) the standard formula in the 
quantum theory, which gives the cross-section for a transition in which a free electron 
with an initial momentum fp (in atomic units) becomes an electron with a momentum k; 
by the absorption of a quantum of appropriate energy from an incident beam of unit 
specific intensity, can be written in the form® 


2re? 


x k2| 2) | (1-1, 1, R2) |2} cm’, 
l=1 


(6) 


where the various matrix elements (/, k?|7#|/ — 1, 3) have to be evaluated in atomic 


units and 
Ak? = k?— k?. (7) 


Equation (6) gives the cross-section for a single hydrogen atom in the ground state for 
the transition in question when there is one electron with momentum &p in unit volume. 
Moreover, the wave length, X, of the radiation for which formula (6) gives the atomic 
absorption coefficient is 011.3 


(AM) (8) 
In the Hartree field of a hydrogen atom the acceleration is given by 
p= (9) 
r r 


and the matrix elements which have to be evaluated are of the type 
- 1 1 


To obtain the absorption coefficients appropriate for an assembly in which we have a 
Maxwellian distribution of electrons corresponding to a temperature 7, we must average 
a(k2; Ak®) over all initial k? (and for a fixed Ak?) with the weight function 


197.8 (11) 


*Cf. J. A. Gaunt, Phil. Trans. R. Soc., London, A, 229, 163, 1930; see particularly formula (5.24) on 
p. 194; also Wheeler and Wildt, op. cit., eq. (6) on p. 284. 
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where 5040 
(12) 


In astrophysical applications it is convenient to express the free-free atomic absorp- 
tion coefficients as per unit electron pressure. If we denote by «(Ak?) the corresponding 
atomic absorption coefficient after averaging over a Maxwell distribution of initial veloc- 
ities, we find that we can write our basic formula in the form 


(AR) 
(13) 
cm 
where 
f (R2) = 100k, (14) 


4. Details of the evaluation: tables of the necessary matrix elements.—The problem of 
numerically evaluating absorption coefficients for free-free transitions is a specially 
troublesome one, since for each wave length the coefficients must be computed for a range 
of initial velocities sufficient to allow for the averaging over Maxwell distributions for 
various temperatures. If the matrix elements for all the necessary values’ of kj and Ak? 
must be individually evaluated, then the number of wave functions which would be 
needed will be so many as to make the preblem an impracticable one on this score alone. 
It may, therefore, be useful to place on record the details of our procedure in this in- 


stance. 
First, we may observe that the infinite series in equation (13) is so rapidly convergent 


that we may ignore all terms except the first one. This may be verified in the following 
pe already indicated in § 2 that for / 2 1 we can use the Born approximation 
without any serious error. We may therefore write 

xi = (TEL) (15) 
Accordingly, for / 2 2, 


1 (16) 
Writing equation (16) in the form 


0 


and integrating by parts, we find after some further reductions that 


(18) 
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The integrals over the Bessel functions which occur in the foregoing equation can be ex- 
pressed in terms of the Legendre functions of the second kind, Q(x), and their deriva- 
tives, Q;(x).!° We find 


1 
(J, k?| = (x) — (x) Ges (x) +; (x), (19) 
where the argument for the Legendre functions is 


kitk:+4 

(20) 

Using equation (19), we can readily verify the fact that all transitions in which the inci- 
dent electron is characterized by an / greater than 1 contribute less than a fraction of 
1 per cent to the sum in equation (13). It is therefore sufficient to consider only the ma- 
trix elements (0, &3|#|1, &?) and (1, &2|#|0, &?). In evaluating these matrix elements it 
was found convenient to distinguish three cases, discussed below. 
i) k? < 0.015 and k? > 0.015.—It has been shown’ that, when k? < 0.015, the radial 
functions x0(r; k?) and x:(r; k®) can be expressed with sufficient accuracy in the form 


xo (7; = — k?Vo(r) (21) 
and 
x1 (7; R?) a(n (22) 


where Ao(k) and A,(k) are the factors which normalize the functions in brackets to unit 
amplitude at infinity and Xo, Yo, X1, and Y; are certain functions which have been 
tabulated. With this representation of the wave functions, the matrix element (0, &3|#|1, 
k?) can be expressed in the form 


1 


(23) 


The corresponding expression for (1, &2|#|0, #7) can be obtained from equation (23) by 
simply interchanging and &. 

Thus, when &? and &? are both less than 0.015, the necessary matrix elements can all 
be evaluated in terms of four definite integrals which are listed in Table 1. 


TABLE 1 


INTEGRALS FOR COMPUTING THE MATRIX ELEMENTS WHEN 3 AND k? 
ARE BOTH LESS THAN 0.015 


Wr. 0.16669 vdr......... 0.09076 


ii) k? < 0.015 and k} > 0.015.—In this case we can use the representation (20) or 
(21) only for the radial function describing the incident electron, and the required matrix 
elements can be expressed in the forms 


(0, #3) = Lf Xoxs (13 rdr— Ri Voxs (13 i) (24) 


1 
Ao (ho) 
10 Cf, G. N. Watson, Theory of Bessel Functions, p. 389, Cambridge University Press, 1944. 
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and 


(1, = FA Xixo (1; — f Vixo (9; (25) 


Ay (ko) 


The integrals which occur on the right-hand side of equations (24) and (25) have been 
evaluated for various values of k? 2 0.015 and are listed in Table 2; for these values 
of k? and for all k§ < 0.015, the required matrix elements can be found without further 
numerical integrations according to equations (24) and (25). - 

iii) k? 2 0.015 and k? > 0.015.—In this case there is no short cut to the evaluation 
of the matrix elements, and the integrals must all be evaluated individually. However, 
for this purpose we have only the radial functions tabulated in the paper by Chandrasek- 
har and Breen. 

Our procedure, then, for determining the various matrix elements with assigned &? 
and Ak? was as follows: 


TABLE 2 
INTEGRALS FOR COMPUTING THE MATRIX ELEMENTS WHEN &? <0.015 
0.64185 0.14271 0.0033763 0.0018942 
.71329 15557 .0067421 .0037565 
.74163 . 15863 .0101066 .0055872 
.75618 . 15858 .013458 .0073854 
. 76622 . 15650 .017920 .0097336 
.77079 . 15328 .022366 .012025 
.77310 . 14360 .033392 .017504 
.77164 . 13840 .038851 .020113 
. 76954 .13170 .044272 .022636 
0.76357 0.12293 0.054972 0.027421 


First, the matrix elements listed in Table 3 were computed with the known radial 
functions. (This table also includes the matrix elements computed according to eqs. 
[23], [24], or [25]). The values given in Table 3 show sufficiently smooth differences with 
both arguments to enable satisfactory interpolation for any intermediate value. Using 
these directly computed values, we next found, by interpolation, all the matrix elements 
which were needed for the evaluation of «(Ak?) according to equation (13); the matrix 
elements found in this manner are given in Table 4, which was then the basis for our fur- 
ther calculations. 

5. The continuous absorption coefficient of H~.—Even with as complete a table of ma- 
trix elements as Table 4, the evaluation of the continuous absorption coefficient for 
free-free transitions is a tiresome matter, since, for each wave length and temperature 
for which the coefficient is desired, a numerical integration must be performed. The co- 
efficients tabulated in Table 5 were determined in this manner according to equation (13). 

In Table 5 we have also entered the absorption coefficients due to the bound-free 
transitions for \ < 16,500 A. These latter coefficients were obtained by multiplying the 
cross-sections derived from those given in Paper II (Table 3) by the factor 


(0) = 4.158 K 1019957 ¢1.7266 (26) 


which gives the number of H~ ions present per neutral hydrogen atom and unit electron 
pressure. The contributions to the continuous absorption coefficient of H~ by the 
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TABLE 3 
THE MATRIX ELEMENTS (0, || 1, #3) AND (1, |7|0, 2?) 
») EVALUATED BY EXACT NUMERICAL INTEGRATIONS 
| (0, 0.0025|#| | (1, 0.0025|*#| | (0, 0.005|*| | (1, 0.005|¥| | (0, 0.0075|F| | (1, 0.0075 
es | 1, k2) 0, 2) 1, k2) 0, 1, k2) 0, 
} 
mn .001577 .0007681 .002033 .001536 0.002304 0.002304 
.002103 .0008293 .002711 .001659 .003073 .002488 
.002629 .0008738 .003390 | .001747 .003842 .002621 
.003153 .0009022 .004066 | .001804 .004608 .002707 
2 .006297 .001003 008119 | .002005 .009201 .003008 
009439 001043 .01217. ~=—.002085 .01379 .003128 
.01257 .001063 .01621 | .002126 .01837 .003189 
| 01674 .001077 .02158 | .002154 .02446 .003232 
02089 .001084 .02694 .002167 .03053 .003251 
02605 .001087 .03360 .002174 .03808 .003261 
.03119 .001087 | .04022 .002174 .04558 .003261 
= .03629 .001085 .04680 | .002170  .05304 003255 
.04135 .001082 .05333 .002164 | .06044 .003246 
0.05135 0.001074 0.06622 0.002147 | 0.07506 0.003221 
(0, 0.01 (1, 0.01|#| (0, 0.0125|*| | (1,0.0125|¥| | (0, (1, 0.015 |#| 
1, 0, k2) 1, k2) 0, 1, 0, 
ee .004975 .003609 0.005241 0.004513 0.005412 0.005412 
.009935 .004011 .010467 .005015 .010808 .006015 
.01489 .004171 .01569 005215 .01620 .006254 
01983 .004253 .02090 005317 .02158 .006377 
02641 .004309 .02783 005388 .02873 .006463 
03296 .004335 .03473 005421 .03586 .006502 
l 04112 .004348 .04332 005437 .04473 .006522 
.04922 .004349 .05186 005438 .05355 .006523 
.05727 .004341 .06035 005428 .06232 .006516 
1 06527 .004329 .06877 .005414 .07102 .006494 
0.08105 0.004296 0.08541 0.005372 0.08820 0.006445 
| 
| 
(0, 0.02 (1, 0.02 | (0, 0.025|F| | (1,0.025|F| | (0, 0.03|F| (1, 0.03 
1, k2) 0, 1, k?) 0, 2) 1, k2) 0, k?) 
0.007594 0.007593 0.007842 0.009471 | 0.008021 0.01135 
.01517 .008260 01567 | .010303 | .01603 .01237 
.02270 .008506 .02345 .010609 | .02398 | .01274 
.03023 .008619 .03123 01075 .03194 .01291 
.03773 .008672 .03898 .01082 | .03987 .01299 
4 .04707 .008699 .04862 .01085 | .04974 .01303 
.05635 .008701 .05821 .01085 .05955 .01303 
> .06557 .008685 .06774 01083 | .06930 .01301 
.07473 .008663 | .07721 01081 | .07899 .01298 
| 0.09277 | 0.008597 0.09590 0.01072 0.09811 0.01288 
| | | 
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TABLE 3—Continued 
(0, 0.0357] | (1,0.035|*] | (0, 0.04]F| (1, 0.04 |#| (0, 0.06 | (1, 0.06|7| 
1, k2) 0, 1, 0, 1, 0, k?) 
.02440 .01487 .02470 .01700 0.02544 0.02544 
.03249 .01507 .03290 .01723 .03388 02578 
.04056 .01516 .04107 01733 ° .04230 02595 
.05060 01521 .05123 .01739 .05279 02603 
.06058 01521 .06135 .01740 .06321 02605 
.07050 01519 .07139 .01737 .07357 02600 
ee .08036 01515 .08137 .01732 .08387 02594 
ee 0.09982 0.01504 0.10109 0.01720 . 10422 .02576 
(0, 0.08 (1, 0.08|7| (0, 0.1[F| (1, 0.1|¥| (0, 0.125|#| | (1, 0.125 
1, 0, &) 1, 0, 1, 0, 
.05352 .03469 .05388 .04332 0.05408 0.05408 
.06410 .03471 .06453 .04339 .06478 .05417 
.07462 .03466 .07513 .04330 .07543 .05406 
.08507 .03458 .08566 .04321 .08602 05395 
1057 .03434 . 1065 .04292 . 1070 05361 
0.1461 0.03372 . 1472 .04216 . 1480 .05270 
(0, 0.15 |#] (1, 0.15|#| (0, 0.175|#] | (1, 0.175|F| (0, 0.2{F| (1, 0.2|#| 
1, 0, 1, 0, 1, 0, 
0.06485 0.06485 0.06478 0.07552 0.06467 0.08614 
.1071 .06428 .1071 .07490 . 1070 .08547 
. 1483 .06321 . 1483 .07369 . 1482 .08414 
. 1881 .06198 0.1882 0.07230 . 1882 .08259 
(0, 0.25|#| (1, 0.257 | (0, 0.35 (1, 0.35 (0, 0.45 (1, 0.45 
1, 0, 1, 0, 1, 0, 
. 1876 . 1030 1856 1434 0.1829 0.1829 
.2071 .1021 2050 1422 . 2023 .1815 
2451 .1001 2430 1396 . 2402 .1785 
.2817 .09809 2798 1370 .2769 .1754 
0.3169 0.3124 0.1722 
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(0, 0.5|#|1, (1, 0.5|#|0, &2) - (0, 0.6|#|1, &2) (1, 0.6|#|0, 2) 

.1976 0.235 0.2351 
.2752 1942 271 .2314 

Serre 0.3107 0.1909 0.3070 0.2277 


bound-free and free-free transitions can thus be brought together in the same system 
(cm*/dyne). 

In Figure 3 we have illustrated the results of Table 5 for T = 6300° (6 = 0.8). For 
comparison we have also included the earlier evaluations. 

In Table 6 the coefficients given in Table 5 have been reduced by the stimulated emis- 
sion factor 

(1 31.320 Ak?) (27) 

which is quite appreciable in the infrared. The values given in Table 6 are such as to 
enable satisfactory interpolation for most astrophysical applications. 

Finally, in Table 7 we give the net absorption coefficient of H~, including both the 
bound-free and the free-free transitions. 


x,s10"* 4 
4 
a “ 
4 
0.4 0.6 0.8 1.0 12 1.4 16 1.8 2.0 2.2 24 26 


a» ing 


Fic. 3.—The continuous absorption coefficient of H~ per neutral hydrogen atom and per unit electron 

ressure for a temperature of 6300° K. The results of earlier determinations of the same quantities are 

included for comparison. The curve M & B is derived from the cross-sections of Massey and Bates for 
the bound-free transitions, and the curve W & W corresponds to the results of Wheeler and Wildt. 
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6. Concluding remarks.—It is not our intention here to make detailed applications of 
our new absorption coefficients to various astrophysical problems. But one comparison 
is of interest. In Figure 4 we have compared the theoretical variation of the continuous 
absorption coefficient of H~ with wave length for JT = 6300° with that derived by 
Chalonge and Kourganoff in their recent discussion of the continuous spectrum of the 
sun for the same temperature. It is seen that the agreement is very satisfactory. In any 
event, there is no basis for the conclusion that the free-free transitions of H~ are not 


2.4 


Fic. 4.—A comparison of the theoretical variation of the continuous absorption coefficient with wave 
length of H~ for T = 6300°, with the corresponding variation derived empirically from the solar con- 
tinuous spectrum by Chalonge and Kourganoff. 


adequate to account for the observed amount of absorption in the red beyond A 10,000 A. 
Indeed, it would appear that the negative ion of hydrogen by itself is able to account 
quantitatively for the entire continuous spectrum of the sun over the range of wave 
lengths \\ 4000 A to 25,000 A. 

On the physical side we should, however, emphasize that, while our new coefficients 
for the free-free transitions of H~ are probably to be trusted generally, the importance of 
exchange and polarization for the slow s-electrons may lead to further changes. It is 
likely that these effects will have a tendency to reduce our coefficients somewhat. We 
hope to return to these questions in the near future. . 
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THE CONTINUOUS SPECTRUM OF THE SUN AND THE STARS 


S. CHANDRASEKHAR AND GuIpo Mijncu! 
Yerkes Observatory 
Received September 5, 1946 


Using the recent determination of the continuous absorption coefficient of H~ by Chandrasekhar 
and Breen, we have shown that the dependence of the continuous absorption coefficient with wave 
length in the range 4000-24,000 A, which can be inferred from the intensity distribution in the con- 
tinuous spectrum of the sun, can be quantitatively accounted for as due to H~; and, further, that the 
color temperatures measured in the wave-length intervals 4100-6500 A (Greenwich) and 4000-4600 A 
(Barbier and Chalonge) for stars of the main sequence and of spectral types AO-GO can also be interpreted 
in terms of the continuous absorption of H~ and neutral hydrogen atoms. 

The problem of the discontinuities at the head of the Balmer and the Paschen series is also briefly 
considered on the revised physical theory of the continuous absorption coefficient. 


1. Introduction.—The two principal problems in the theory of the continuous spec- 
trum of the stars are, first, to identify the source of the continuous absorption in the 
solar atmosphere which will account for the intensity distribution in the continuous 
spectrum of the sun and the law of darkening in the different wave lengths and, second, 
to account for the observed relations between the color and the effective temperatures of 
the stars. In this paper we shall show that the major aspects of these two problems find 
their natural solution in terms of the continuous absorption coefficient of the negative 
hydrogen ion as recently determined by Chandrasekhar and Breen.” More particularly, 
we shall show that the dependence of the continuous absorption coefficient with wave 
length in the range 4000-24,000 A, which can be deduced from the solar data, can be 
quantitatively accounted for as due to H~; and, further, that the color temperatures 
measured in the wave-length intervals 4100-6500 A (Greenwich*) and 4000-4600 A 
(Barbier and Chalonge*) for stars of the main sequence and of spectral types AO-GO 
can also be interpreted in terms of the continuous absorption of H~ and neutral hydro- 
gen. 

In addition to the two problems we have mentioned, we shall also consider some 
related questions concerning the discontinuities at the head of the Balmer and the 
Paschen series. 

2. The mean absorption coefficients of H~ and H.—As is well known, the character of 
the emergent continuous radiation from a stellar atmosphere is determined in terms of 
the temperature distribution in the atmosphere; and, as has recently been shown,° the 
temperature distribution in a nongray atmosphere will be given approximately by a 


formula of the standard type 
(r+ q[7)), (1) 


where 7, denotes the effective temperature and g(r) a certain monotonic increasing func- 
tion of the optical depth 7, provided that the mean absorption coefficient x, in terms 
of which 7 is measured, is defined as a straight average of the monochromatic absorption 
coefficient weighted according to the net flux F{!) of radiation of frequency » in a gray 
1 Fellow of the John Simon Guggenheim Memorial Foundation at the Yerkes Observatory. 
2Ap. J., 104, 430, 1946. 


3 Sir Frank Dyson, Observations of Color-Temperatures of Stars, 1926-1932, London, 1932; also M.N., 
100, 189, 1940. 


4 Ann. d’ap., 4, 30, Table 4, 1941. 


5S. Chandrasekhar, Ap. J., 101, 328, 1945. This paper will be referred to as “Radiative Equilibrium 
VII.” 
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atmosphere and if, further, x,/x is independent of depth. On this approximation, then, 
the emergent intensity in a given frequency and in a given direction will depend not 
only on the continuous absorption coefficient at the frequency under consideration but 
also on the mean absorption coefficient x over all frequencies. 

As the discussion in this paper will establish for the stellar atmospheres considered, 
the contributions to « in the visible and the infrared regions of the spectrum are essential- 
ly from only two sources: H~ and the neutral hydrogen atoms. The cross-sections for 
the absorption by H~ for various temperatures and wave lengths have been tabulated by 
Chandrasekhar and Breen in Table 7 of their paper, while those for hydrogen can be 
found from the formulae of Kramers and Gaunt, standardized, for example, by B. 
Strémgren.® However, the evaluation of the mean absorption coefficient for wave lengths 
shorter than 4000 A is made uncertain on two accounts: First, there is the absorption by 
the metals and the excessive crowding of the absorption lines toward the violet, which is 
particularly serious for spectral types later than FO; and, second, there is the absorption 
in the Lyman continuum. On both these accounts the true values of x will be larger than 
those determined by ignoring them. But the exact amount by which they will be larger 
will be difficult to predict without a detailed theory of “‘blanketing,’’’ on the one hand, 
and without going into a more exact theory® of radiative transfer than represented by 
the approximations leading to equation (1), on the other. However, since in this paper 
our primary object is to establish only the adequacy of H~ as the source of absorption in 
the solar atmosphere over the entire visible and infrared regions of the spectrum and the 
corresponding role of H~ and H for stellar atmospheres with spectral types A2—G0, it 
appeared best to ignore the refinements indicated and simply determine x by weighting 
x, due to H~ and H (without the Lyman absorption) at the conditions prevailing at 
7 = 0.6 by the flux F{” at this level. For only in this way can we use the solution to the 
transfer problem in the form of equation (1) in a consistent manner. It should, however, 
be remembered that the effects we have ignored may easily increase x, determined in 
terms of H~ and H (without the absorption in the Lyman continuum) by factors of the 
order of 1.5 and probably not exceeding 2.!° _ 

Turning our attention, next, to the evaluation of x, we may first observe that, since 
our present method of averaging is a straight one, the contributions to x, from different 
sources are simply additive. We may, accordingly, consider the mean absorption co- 
efficient of H~ and H separately. 

Now the absorption coefficient of H~, including both the bound-free and the free-free 
transitions, is most conveniently expressed as per neutral hydrogen atom and per unit 
electron pressure in the unit cm*/dyne. The monochromatic coefficients «,, after allowing 
for the stimulated emission factor (1 — e~”’/*”), are tabulated in Chandrasekhar and 
Breen’s paper for various values of 6(= 5040/7). If we now denote by a(H~) the average 


5“Tables of Model Stellar Atmospheres,” Publ. mind. Meddel. Kobenhavns Obs., No. 138, 1944. 
7™Cf. G. Miinch, Ap. J., 104, 87, 1946. 
8 Such as, e.g., the (2, 2) approximation given in ‘Radiative Equilibrium VII,” § 6. 


® The choice of r = 0.6 for the “representative point” was made after some preliminary trials (cf. 
G. Miinch, Ap. J., 102, 385, 1945, esp. Table 3), though it is evident on general grounds that a level 
such as r = 0.6, where the local temperature is approximately the same as the effective temperature, 
would be the correct one in the scheme of approximations leading to eq. (1). 


10 Tn all earlier evaluations of x the absorption in the Lyman continuum did not, indeed, play any 
role. This was due to the manner in which x was defined in those investigations as the Rosseland mean. 
But in “Radiative Equilibrium VII” it has been shown that there is no justification for taking the Rosse- 
land means as they have been hitherto. Since the method of averaging, by which we have now replaced 
the Rosseland mean, is a straight one, it is no longer permissible simply to ignore the absorption in the 
Lyman continuum. At the same time, it is not possible to take it into account satisfactorily in the (2, 1) 
approximation leading to eq. (1). We should have to go at least to the (2, 2) approximation of “Radia- 
tive Equilibrium VII.” 
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value of the coefficients x}, weighted according to the flux F{” in a gray atmosphere at 
7 = 0.6, where the temperature is approximately the effective temperature 7., then 
the contribution x(H~) to the mass absorption coefficient x by H™ is given by 


where mz is the mass of the hydrogen atom, #, the electron pressure, and xq the degree 
of ionization of hydrogen under the physical conditions represented by 7, and #,.!! The 
values of a(H~—) found by graphical integration in accordance with the formula 
1 
)== (0.6 3 
a(H-) (0.6) dy (3) 


for various values of 6 = @, are given in Table 1. 


TABLE 1 
THE MEAN ABSORPTION COEFFICIENTS a(H~) AND a(H) 

6 a(H-) | a(H) | 6 a(H-) a(H) 
0.563 10-76 1.65X 10-2? 6.08 X 10-6 5.52 10-27 


Similarly, the contribution to x by hydrogen can also be expressed in the form 


z(H) = a (H), (4) 
where 
_ FY (0.6) 
a(H) = f ada, (5) 


where a = hv/kT, and f and D are certain functions of temperature and frequency, re- 
spectively, which have been tabulated by Strémgren.’ For reasons which we have al- 
ready explained, we do not include the Lyman absorption in evaluating a(H). The values 
of a(#7) for various temperatures are also listed in Table 1. 

In terms of a(H-) and a(H) given in Table 1, we can determine the combined mass 
absorption coefficient x according to 


[a (H-) p-+a (H)] . (6) 


K = 


H 


Values of x(H~), x(H), and x determined in accordance with the foregoing equations 
for various temperatures and electron pressures are given in Table 2. 

3. The continuous absorption in the solar atmosphere—As we have already stated in 
the introduction, one of the principal problems in the interpretation of the solar spectrum 
is the identification of the source of absorption which will predict the same dependence 


11 Tt will be noted that, in writing the mass absorption coefficient in the form (2), we have assumed 
the preponderant abundance of hydrogen in the stellar atmosphere. 


12 See the reference quoted in n. 6. 
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of the absorption coefficient with wave length in the range 4000—24,000 A, which can 
be inferred from the intensity distribution in the continuous spectrum and the law of 
darkening in the different wave lengths. While the amount and variation of the continu- 
ous absorption in the spectral region AX 4000—24,000 A can be deduced in a variety of 
ways,'* it appears that, for the purposes of the identification of the physical source of 
absorption, it is most direct to adopt the following procedure: 


TABLE 2 


THE MEAN MASS ABSORPTION COEFFICIENTS k(H~), x(H), AND k FOR VARIOUS 
TEMPERATURES AND ELECTRON PRESSURES 


be=1 be=10 be = 10" be = 10° be 
| | 3.70x10-+ | 4.9810 2.12 32.1 
x(H) 1.71X10 | 1.08 14.5 61.8 93.5 
i 1.71X107 | 1.08 14.6 63.9 126 
| 3.93x10-* | 2.60x107 | 5.82X107 6.68 68.4 
eH) | 4.18X107 | 2.76 6.19 7.10 7.27 
418X107 | 2.79 6.77 13.8 75.7 
| 9.08X10-? | -1.28K107 | 1.34 13.4 134 
eH) | 4.54X107 | 6.41X107 | 0.67 0.7 1 
463X107 | 7.69X107 | 2.01 14.1 135 
eH-) | 2.26x10-? | 2.32x107 | 2.32 23.2 232 
; 6.92X10- | 2.80X10- | 2.37 23.2 232 
| 3.63X10-2 | 3.63X107 | 3.63 36.3 363 
3.96X10-2 | | 3.63 36.3 363 
e(H-) | 5.57x10-2 | 5.57x107 | 5.57 55.7 557 
5.59x10-2 | 5.57X107 | 5.57 55.7 557 
| 1.20x107 | 1.20 12.0 120 1200 
| 2.33x107 | 2.33 23.3 233 2330 


We compare the observed intensity distribution in the emergent solar flux F (obs.) 
with the flux F{? (0) to be expected in a gray atmosphere." It is evident that the de- 
partures, 


A log Fx = log Fy (obs.) — log Fx” (0) , (7) 


must be related more or less directly with the dependence of the continuous absorption 
coefficient x, with wave length. Indeed, in the approximations leading to the tempera- 
ture distribution (1) this relation must be one-one, since, with the adopted definition 
of x, the temperature distributions in the gray and the nongray atmospheres agree." 
This suggests that, with the known value of x, due to H~, we compare the predicted de- 


13 G, Miilders, Zs. f. Ap., 11, 132, 1935; G. Miinch, Ap. J., 102, 385, 1945; D. Chalonge and V. Kour- 
ganoff, Ann. d’ap. (in press). 

4 The values of F{!) (0) can be readily derived from the entries along the line r = 0 in Table 2 of 
“Radiative Equilibrium VII.” 

15 Cf, the remarks in italics on p. 343 in “Radiative Equilibrium VII.” 
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partures from FX” (0) with those observed. The only uncertainty in these predictions will 
be of the nature of a “zero-point” correction, since a value of x different from the one 
adopted will lead to an approximately constant additive correction to A log Fy.1° 

In order, then, to make the comparison suggested in the preceding paragraph, we 
need to determine A log F) in terms of K due to H~ and an adopted x. Assuming in the 
first instance that the contribution to x is only H~, we find that 


a (H~) =5.62 X 10~-* cm*/dyne (8) 
for an adopted value of 
6, =0.8822. (9) 


The value of xx for 6 = 0.8822 can be found by, simple interpolation in Table 7 of 
Chandrasekhar and Breen’s paper. The ratios «x/a(H~) derived in this manner are 


TABLE 3 


THE PREDICTED DEPARTURES [LOG Fx (obs.)—LOG F{!)(0)] FROM GRAYNESS OF THE 
SOLAR ATMOSPHERE DUE TO THE ABSORPTION BY H~ 


LoG F) (THEO.) A toc Fy 
A Loc Ff 1)(9) 
x(H-) 1.42x(H-) - - 

«=x(H-) |x=1.42K(H-) k=x(H-) | =a1.42«(H-) 

See 0.686 0.483 14.489 14.532 14.358 +0.131 +0.174 
a re 0.881 .620 14.431 14.534 14.399 + .032 + .135 
Se 1.029 iy a 14.341 14.423 14.350 — .009 + .073 
ee 1.132 .197 14.259 14.329 14.285 — .026 + .044 
1.188 .837 14. 164 14.224 14.194 — .030 + .030 
eee 1.183 .833 14.070 14.128 14.098 — .028 + .030 
1025 .792 13.982 14.034 13.999 — .017 + .035 
LE nears 1.028 .724 13.897 13.946 13.901 — .004 + .045 
cS 0.911 .642 13.817 13.862 13.805 + .012 + .057 
| 0.788 13.742 13.784 + .030 + .072 
| Sree 0.651 .458 13.674 13.713 13.624 + .050 + .089 
Dae 0.523 .368 13.608 13.643 13.537 + .071 + .106 
0.481 .339 13.532 13.566 13.454 + .078 + .112 
Sree 0.486 . 342 13.448 13.480 13.374 + .074 + .106 
Se 0.516 . 363 13.362 13.395 13.297 + .065 + .098 
0.562 .396 13.279 13.313 133225 + .056 + .090 
0.618 13.197 13.230 13. 152 + .045 + .078 
S| | eae 0.679 .478 13.118 13.149 13.083 + .035 + .066 
0.749 13.042 13.073 13.016 + .026| + .057 
LL ee 0.820 0.577 12.971 13.003 13.952 +0.019 +0.051 


given in Table 3 for various values of \. With these values of x,/a(H~), the theoretical 
determination of A log F, is straightforward with the help of the nomogram of Burk- 
hardt’s table,!” which one of us has recently published.'* The results of the determina- 
tion are given in Table 3. In Figure 1 we have further compared the computed de- 
partures A log F, with those observed.’® It is seen that the predicted variation of the 


16 This is seen most directly in an approximation in which we expand the source function B)(7) asa 
Taylor series about a suitable point and determine the emergent flux in terms of it (see, e.g., A. Unsdld, 
Physik der Sternatmosphdren, p. 109. eq. [31.18], Berlin, 1938). 

17 Zs. f. Ap., 13, 56, 1936. 18 G. Miinch, Ap. J., 102, 385, Fig. 2, 1945. 

19 For \ > 9000 A the observed departures were obtained from a reduction of the solar data by 
M. Minnaert, B.A.N., 2, No. 51, 75, 1924; see also Unsdld, op. cit., p. 32. For \ < 9000 A the reduction 
of G. Miilders (dissertation, Utrecht, 1934) was used. 
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departures runs remarkably parallel with the observed departures over the wave- 
length range 4000—20,000 A. (The observational data do not seem specially reliable for 
\ > 20,000 A.) However, the absolute values of the predicted departures are system- 
atically less than the predicted departures by approximately a constant amount, indicat- 
ing a zero-point correction in the sense that the adopted value of x as due to H~ alone 
is somewhat too small. The calculations were accordingly repeated for other slightly 
larger values of x, and it was found that with x = 1.42 «(H-), the predicted and the 
observed departures agree entirely within the limits of the observational uncertainties 
over the whole region of the spectrum in which H~ contributes to the absorption. The 
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Fic. 1.—Comparison of the observed and the theoretically predicted departures [log F, — log F} (0)] 
from a gray atmosphere due to the absorption by H~. The circles represent the observed departures of 
the solar emergent flux from that of a gray atmosphere, while Curves I and II are the theoretically 
derived departures on the two assumptions k = k(H~) and k = 1.42 K(H~). 
agreement is, in fact, so striking that we may say that H~ reveals its presence in the solar 
atmosphere by its absorption spectrum. _ ni 

We may finally remark on the value of x = 1.42 x(H-), indicated by the comparisons 
we have just made: it can, in fact, be deduced empirically from the solar data on the 
continuous spectrum that the absorption in the violet (A < 4000 A) must increase the 
value of x derived from the visible and the infrared regions of the spectrum by a factor 
of the order of 1.5.”° 

4. The predicted color-effective temperature relations: comparison with observations.— 
All earlier attempts”! to predict the color temperatures in the region \d\ 4000-6500 A 
for stars of spectral types AO-GO in agreement with the observations have failed. This 
failure in the past has been due to the following circumstance: The observed relation 
between the color and the effective temperatures and, in particular, the fact that T, > T, 


20 Cf. G. Miinch, Ap. J., 102, 385, 1945, esp. the remarks preceding eq. (19) on p. 394. 
21R. Wildt, Ap. J., 93, 47, 1941, and Observatory, 64, 195, 1942; R. E. Williamson, Ap. J., 97, 51, 1943. 
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implies that the continuous absorption coefficient is an increasing function of A in the 
spectral region observed. But the physical theory on which the calculations were made 
placed the maximum of the absorption-curve in the region of \ 4500 A; this was incom- 
patible with the observations and, moreover, predicted color temperatures less than the 
effective temperatures, contrary to all evidence. Indeed, on the strength of this dis- 
crepancy, it was concluded that H~ as a source of absorption was inadequate even in the 
region AA 4500-6500 A, and the existence of an unknown source operative in this region 
was further inferred. However, later evaluations” of the bound-free transitions of H~ 
showed the unreliability of earlier determinations and placed the maximum of the 
absorption-curve in the neighborhood of \ 8500 A. The addition of the free-free transi- 
tions pushes this maximum only still further to the red. It is therefore evident that on 
the revised physical theory we should be able to remove the major discrepancies of the 
subject. We shall now show how complete the resolution of these past difficulties is. 

From the point of view of establishing the adequacy of the physical theory in the 
region AA 4000-6500 A, it is most instructive to consider the theoretical predictions for 
color temperatures which can be directly compared with the color determinations at 
Greenwich,’ for the Greenwich measures are based on the mean gradients in the wave- 
length interval 4100-6500 A, and it is in the prediction of these colors that the earlier 
calculations were most discordant.”* 

Now, from the Planck formula in the form 


1 
= 1? (10) 
it readily follows that 
Defining the gradient 
(12) 
in the usual manner, we can write 
1 d logit _ ) 
2-903). (13) 


If Fy, and Fy, are the emergent fluxes at two wave lengths \; and ), and if ¢ is the mean 
gradient in this wave-length interval, then we can write, in accordance with equation (13) 


$= 


where X,, denotes an appropriate mean wave length for the interval to which the gradi- 
ent ¢ refers. Equation (14) can be re-written in the following form: 


1 A logio F 
(15) 
According to equation (15), the theoretical determination of color temperatures will 


proceed by determining, first, the gradient ¢ from the values of F, at the end-points 
of the wave-length interval and then determining the temperature which will give this 


gradient. 


$= 


22S. Chandrasekhar, Ap. J., 102, 223, 395, 1945 
23 Cf. Fig. 6 in Williamson’s paper (0. cit.). 
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For the Greenwich measures \,, = 0.55 wu, and equation (15) becomes 


¢ (Greenwich) = 2.75 — 2.56 logio 


F¢s007 ' 
provided that, in determining the gradients, wave lengths are measured in microns.”4 
In Table 4 we have listed the values of «,/« for the wave lengths 4100 A and 6500 A 
for various values of 6, and p,. In terms of these values the determination of the fluxes at 
the two wave lengths 4100 A and 6500 A is straightforward with the help of Burkhardt’s 
table. The gradient ¢ then follows according to equation (16) and, from that, the color 


temperature. The reciprocal color temperatures @, = 5040/7, derived in this manner 
are given in Table 5. The resulting color—effective temperature relations are illustrated 


TABLE 4 
«,/% IN MODEL STELLAR ATMOSPHERES 
5040/Te be =10 be =10? be=108 be 
3.45 3.22 2.68 
0.161 0.181 0.332 
» 4000 0.210 0.212 0.238 0.397 
4600 0.316 0.318 0. 346 0.513 
» 6500 0.826 0.834 0.848 0.991 
1.51 1.51 1.55 
8203 0.490 0.526 0.788 
4.04 2.60 1.09 
N3647 0.178 0.427 0.685 
» 4000 0.148 0.205 0.481 0.770 
4600 0.221 0.278 0.585 0.901 
» 6500 0.579 0.637 0.935 1.25 
1.27 1.36 1.45 
8203 { 0.424 0.885 1.36 
(6.06 2.82 0.961 0.684 
3647 | 10.176 0.461 0.625 0.645 
» 4000 0.208 0.521 0.695 0.725 
4600 0.276 0.617 0.810 0.848 
6500 0.546 0.901 Lil 1.14 
0.889 1.105 1.22 
[A 8205 ‘0. 356 0.893 1.19 1.24 
1.54 0.720 
d 3647 0.518 0.601 0.614 0.611 
» 4000 0.582 0.671 0.685 0.685 
4600 0.681 0.781 0.800 0.800 
6500 0.961 1.06. 1.09 1.09 
1.178 
0.708 0.642 
3647 10,681 0.636 0.636 
» 4000 0.690 0.690 0.690 0.690 
4600 0.800 0.800 0.800 0.800 
» 6500 1.09 1.09 1.09 1.09 
{t.16 116 1.16 


24 The constant cz in eq. (12) then has the value 14,320. 
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in Figure 2. For comparison we have also plotted in this figure the Greenwich determina- 
tions for stars on the main sequence and of spectral types AO—GO (reduced, however, to 
the Morgan, Keenan, and Kellman system of spectral classification). The color tempera- 
ture of the sun for this wave-length interval is also plotted in Figure 1. It is seen from 
Figure 1 that the agreement between the observed and the theoretical color temperatures 
is entirely satisfactory, particularly when it is remembered that the earlier calculations 
failed even to predict the correct sign for 0, — @,.. It will, however, be noted that the 
observed values of @, for spectral types later than FO are somewhat larger than the pre- 


TABLE 5 


THEORETICAL RECIPROCAL COLOR TEMPERATURES AND THE PREDICTED DIs- 
CONTINUITIES AT THE HEAD OF THE BALMER AND THE PASCHEN SERIES* 


Ge 
be 

0.5 | 0.6 | 0.7 0.8 0.9 1.0 

0.68 0.80 0.88 

| 0.45 0.60 73 80 88 

[0.31] 52 65 14 80 88 

42 56 ‘64 69 ‘75 

[ .39] 59 66 14 80 88 

(Band C) [ .34] 50 58 64 69 

52 59 66 74 80 88 

H-\6{B andC)| 0.41 0.50 0.58 0.64 0.69 0.75 


* @.(G) and 6-(Band C) are the reciprocal color temperatures, 5040/7 c, appropriate for the wave-length intervals 4100-6500 A 
and 4000-4600 A, respectively; Dg and Dp, representing the logarithm of the ratio of the fluxes at the two sides of the series 
limits, are the expected Balmer and Paschen discontinuities, respectively. 


dicted values, though the agreement is as good as can be expected in the case of the sun. 
The reason for this must undoubtedly be the crowding of the absorption lines toward 
the violet in the later spectral types and the consequent depression of the continuous 
spectrum in this region. The correctness of this explanation is apparent when it is noted 
that in the case of the sun, in which allowance has been made for this effect of the lines 
on the continuum, the discordance is not present. 

Comparisons similar to those we have just made also can be made with the measure- 
ments of Barbier and Chalonge‘ on the color temperatures based on the observed 
gradients in the wave-length interval 4000-4600 A. The formula giving the theoretical 
gradient for this wave-length interval takes the form 

F 4600 


(17) 


¢(BandC) = 2.175 —7.06 logis 
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The values of «x/x at 4000 A and \ 4600 A are given in Table 4, and the reciprocal 
color temperatures derived from these in Table 5. The results are further illustrated in 
Figure 3, where the theoretical relations for various electron pressures are compared 
with the measures of Barbier and Chalonge (reduced also to the Morgan, Keenan, and 
Kellman system of spectral classification). It is seen that the general agreement is again 
good, though there are now somewhat larger differences between the computed and the 
observed color temperatures for spectral types later than FO than were encountered in 
the comparison with the Greenwich colors. This must again be due to the crowding of 
the absorption lines toward the violet in the later spectral types and the further fact that 


0.5 0.6 0.7 0.8 0.9 1.0 


Fic. 2.—Comparison of the observed and the theoretical color effective-temperature relations for the 
wave-length interval 4100-6500 A. The ordinates denote the reciprocal color temperatures and the ab- 
scissae denote the reciprocal effective temperatures (@ = 5040/7). The circles represent the Greenwich 
color determinations reduced to the Morgan, Keenan, and Kellman system of spectral classification. 


The double circle represents the sun. 


the base line for the Barbier and Chalonge colors is much shorter than that for the Green- 
wich colors. 

5. The discontinuities at the head of the Balmer and the Paschen series.—With the 
physical theory of the continuous absorption coefficient now available, we can also pre- 
dict the extent of the discontinuities which we may expect at the head of the Balmer 
and the Paschen series of hydrogen. For this purpose the values of «,/x on the two 
sides of the series limits are also given in Table 4. From these values it is a simple matter 
to estimate the discontinuities which will exist at the head of the Balmer and the 
Paschen series, and they are given in Table 5. The results for the Balmer discontinuities 
are further illustrated in Figure 4, in which the discontinuities measured by Barbier and 
Chalonge‘ for various stars are also plotted. The progressive increase of the electron 
pressure as we go from the later to the earlier spectral types is particularly apparent 
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wave-length interval 4000-4600 A. The ordinates denote the reciprocal color temperatures, and the p 
abscissae denote the reciprocal effective temperatures @ = 5040/T). The circles represent the color u 

determinations of Barbier and Chalonge for the wave-length interval 4000-4600 A, reduced to the 

Morgan, Keenan, and Kellman system of spectral classification. ; F 

Sp 

AO Al A2 A3 AS AT FO F2 FS GO 


Fic. 4.—The predicted discontinuities (D) at the head of the Balmer series for various effective tem- 
peratures and electron pressures. The circles represent the discontinuities as measured by Barbier and 
Chalonge. (The double circles represent the observations for supergiants.) 
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from Figure 4; this progression is, moreover, in agreement with what is indicated by 
the color determinations (cf. Figs. 2 and 3). In Table 6 we give the electron pressures for 
the main-sequence stars of various spectral types estimated in this manner. 

6. Concluding remarks.—While our discussion in the preceding sections has estab- 
lished the unique role which H~ plays in determining the character of the continuous 
spectrum of the sun and the stars, it should not be concluded that the various other 


TABLE 6 
ELECTRON PRESSURES FOR STARS ON THE MAIN SEQUENCE 
Type log fe Type log pe Type log pe 


astrophysical elements of the theory are equally well established. Indeed, the theory of 
model stellar atmospheres as developed by Strémgren in recent years must not only be 
revised on the basis of the new absorption coefficients of H~ but also be advanced still 
further before we can be said to have a completely satisfactory account of all the classical 
problems of the theory of stellar atmospheres. But our discussion in this paper does give 
us the confidence that the continuous absorption by H~ discovered by Wildt must 
provide the key to the solution of many of these problems. 
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NOTES 


LUMINOSITY CHARACTERISTICS ON LOW-DISPERSION 
SPECTRA OF STARS OF TYPES M0-M4 


The luminosity classification of stars of spectral classes MO-M4 has usually been 
carried out with slit spectra by means of line ratios. On slitless spectra of low dispersion, 
however, Lindblad,! Vyssotsky,? and others have found the depression of the continuum 
on the red side of \ 4227 a useful criterion for distinguishing dwarfs from giants. 

In examining the spectra of stars of classes MO—M4 taken with the 4° objective prism 
of the Warner and Swasey Observatory (dispersion at Hy = 283 A/mm) it was observed 
that it is possible to extend the criteria in this spectral region to allow the recognition of 
supergiants. 

P|. XVIII shows two sets of three spectra, the first set representing stars of class MO, 
and the second stars of class M2. The luminosity decreases from top to bottom in each 
set and is indicated by the luminosity-class symbols* at the right of each spectrum. 
RW Cep has been assigned the symbol 0 because on spectrograms of all dispersions the 
criteria of high luminosity are much more strongly developed in this star than in Ia 
stars, such as pu Cep. 


The following are the sensitive criteria: 

1. Intensity of Ca 4227.—In dwarfs this line is the most conspicuous feature in the 
spectrum. It weakens progressively with increasing luminosity until in the more brilliant 
supergiants it is actually less intense than the neighboring absorption features to the 
violet. 

2. Continuous spectrum on both sides of \ 4227.—In dwarfs the violet side is stronger 
than the red side. The two portions are of approximately equal intensity in giants, while 
in supergiants the violet side is definitely the weaker. 

3. G-band.—The break in the spectrum at about \ 4308, corresponding to the red 
edge of the G-band, becomes steadily more pronounced with increasing luminosity. In 
supergiants the strengthening of \ 4290 (due mostly to Cr) gives rise to a conspicuous 
double line, AA 4290-4308. 

4. Intensity of \ 4406.—The strength of the blended line at \ 4406 increases with re- 
spect to the features on both sides of it in passing from dwarfs to giants. The lines in this 
region become indistinct in supergiants, but on slightly underexposed plates this criterion 
is useful in checking the distinction between giants and dwarfs. 

It should be noted that the stars » Cep and RW Cep are among the most brilliant 
supergiants known. For stars similar to a Ori or a Sco (luminosity class Ib) the super- 
giant characteristics are not so pronounced; but on good spectrograms these stars, too, 
can be readily picked out. As a test, two plates covering the field of the double cluster, 
h and x Per, were examined, and the high luminosity of the four M-type supergiants 
known in this field was easily recognized. In addition, it was observed that HD 14826, 
BD+56°609 (m, = 8.5) also has a supergiant spectrum. From these two plates the 


1 Stockholms Obs. Ann., Vol. 12, No. 2, 1935. 
2Ap. J., 97, 381, 1943. 


3 Morgan, Keenan, and Kellman, An Allas of Stellar Spectra, 1943. 
458 


Sh 
} 
i 
j 
i 
t 
F 
j 
* 
\ 


PLATE XVIII 


(a). MO . 


LuMINOsITY EFFECTS ON SMALL-SCALE SPECTROGRAMS OF M-TyPE STARS 
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spectrum was classified approximately as M2 Ib. Its known small proper motion (ue = 
—0*007, us = —0*010)* is consistent with assignment of the star to the extensive group 
of supergiants which are associated with the double cluster. 


P. C KEENAN 
PERKINS OBSERVATORY J. J. Nassau 
AND 
WARNER AND SWASEY OBSERVATORY 
August 1946 


THE RADIAL VELOCITY OF 27 CANIS MAJORIS 
Table 1 contains the radial velocities of 27 Canis Majoris obtained at the Yerkes Observatory 


since October, 1943. The velocity-curve in Figure 1 covers the entire range of our observations 
since 1925. The extraordinary character of this curve defies, at present, all attempts to explain 


TABLE 1 


RADIAL VELOCITIES OF 27 CANIS MAJORIS 


Velocity Velocity 
Date (Km/Sec) Date UF: (Km/Sec) 
11:09 — 21.6 8:30 —74.8 
10:05 — 28.0 9:43 —43.4 
8:39 — 16.2 7:59 —55.3 
8:22 — 10.0 6:52 —83.3 
9:34 — 10.9 6:36 —79.0 
8:50 — 33.0 1945 Mar.. 1:24 —59.0 
8:20 — 34.4 3:02 —62.8 
7:24 — 24.2 1:31 —68.0 
1002 Jan. 7:41 + 0.8 2:52 —70.3 
6:19 — 41.2 —67.6 
1255 — 22.5 | 2:31 —72.3 
a 5:50 — 50.8 1946 Feb. 28........ aide +28.1 
4:42 — 46.7 3:00 +27.3 
5:35 — 31.4 3:20 +28.0 
5:03 — 34.4 1:38 +67.4 
Feb.. 13... 00:54. 5:82 — 30.8 2:00 +62.6 
2:25 — 70.8 1:44 +38.5 
2:26 — 72.1 1:24 +53.0 
11:09 — 97.0 1:28 +49.3 
11:28 — 78.7 1:54 +45.1 
ye 10:20 — 84.3 


4Van Maanen, Mt. W. Conitr., No. 205, 1921. 
1See Ap. J., 97, 453, 1943. 
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the variations. During March, 1946, the spectrum showed faint but unmistakable absorption 
components on the violet sides of the strong hydrogen absorption lines. Such components were 
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Fic. 1.—Radial velocities of 27 Canis Majoris 


observed on previous occasions; they were not measured. I am indebted to several persons for 
the spectrograms of 1946 and to Mrs. Martha Carlson for part of the measurements. 


O. STRUVE 
August 1, 1946 


ON THE POSSIBILITY OF TRACING POLARIZATION EFFECTS 
IN THE ROTATIONAL PROFILES OF EARLY-TYPE STARS 


In connection with the interesting polarization effects in the continuous radiation of 
early-type stars recently predicted by Chandrasekhar’ and established by Edith M. 
Janssen? for the eclipsing variable U Sagittae, the writer wants to call attention to the 
possibility of tracing effects of this kind in the rotational profiles of rapidly rotating 
early-type stars. 

According to Chandrasekhar, the polarization of the continuous radiation increases 
from 0 at the center of the stellar disk to 0.11 at the limb, in consequence of a different 
law of darkening for the radial and for the tangential components. As different parts of 
the stellar disk give different displacements of the spectrum, faint effects of polarization 
should be expected in the rotational profiles of rapidly rotating stars. 

Figure 1, for example, illustrates a favorable case. The upper curve gives the apparent 
rotational profile of the unpolarized radiation from the stellar disk. Superposed on this 
curve there will be profiles in polarized light from the regions (2) (vibration parallel to the 
axis of rotation) and from a fairly narrow central strip (1) (vibration perpendicular to 
the axis of rotation). For a spherical star the integrated polarized light will be the same 
for both vibrations, but, because of the Doppler shift, the two vibrations will appear with 
different intensity in an absorption line. In the center of the line a positive polarization 
will result (vibration parallel with the axis of rotation), whereas a negative polarization 
is to be expected in the wings of the absorption contour. 

The writer has made an estimate of the effects which might appear in the favorable 
case outlined in Figure 1. By using Chandrasekhar’s values of the intensity i(r) and 


1Ap. J., 103, 363, 1946. 
24 p. J., 103, 380, 1946, 


a 
be 
i 
de 
7 m 
di 
i 4800 3000 $500 6000 6500 7000 7300 8000 6500 9000 9600 28430000 0500 +000, 1800 2000 
Ir 
} 
| 
} 
| 
} 
+ 
| 
| 
j 
a 
| 
a 
4 
if 
€ 
Ww 
5 
B 
‘ 
te 
ex 
h 
3 
4 


NOTES 461 


d degree of polarization p(r) at different distances from the center of the stellar disk, the 
mean polarization p, and the total intensity 7, of segment (2) have been calculated for 
different values of a, that is, the distance of the corresponding chord from the center. 
In fact, we have 


cos 20d0 
=— (1) 
finde 


Ps 


Apparent profile of 
unpolarized light. 


(1) Su 


= + 
com ponents 


Fic. 1 
TABLE 1 
| te is 
a 0.04 0.030 0.0012 0.005 0.002 
.08 .022 .0017 .007 .002 
.016 .0019 .008 .003 
0.17 0.011 0.0019 0.008 0.003 


where # is the angle between the radius 7 and the normal a, and dg is the elemental area. 
By a graphical method of integration the values of Table 1 have been obtained, a being 
expressed in the stellar radius as unity and J being the intensity of the total light from 
the stellar disk. 

The quantities P, and P_ in Table 1 are the polarization effects to be expected in a 


4 
42h } 
Z S 
q i) 
I 
P i+ 
7 
| 
3 


& 


462 NOTES 


line contour according to the favorable case outlined in Figure 1. Let us first consider 
P.,, which is the polarization appearing in the line center. 

The unpolarized radiation has the intensity J — 4 p,7,. Because of the line absorption 
(which is assumed to take place above the photosphere), the unpolarized light will as- 
sume the value (1 — x)(J — 4,7.) in the line center, x being the fraction absorbed. As 
the polarized component (/) is assumed to be absent in the line center, whereas the com- 
ponents (2) are supposed to have full intensity, we obtain the degree of polarization P, 
in the line center from 


The values of Table 1 have been calculated on the assumption that x = 0.5. 

A similar calculation can be made for P_. The component (/) and one of the compo- 
nents (2) are here supposed to have full strength, whereas the other component (2) is 
supposed to be absent. We obtain 


Pale 
(3) 


Assuming as a probable value x’ = 0.25, we find the values given in Table 1. 

The calculation has shown that a positive polarization, amounting to 8 promille, and 
a negative polarization of 3 promille can be expected in rotational profiles of early-type 
stars. In the case of flattened objects, somewhat stronger effects may possibly appear. 
With sensitive methods of observation, effects of this order of magnitude might be 
established. 

From a considerable number of plates taken in 1934, the writer? was able to trace 
polarization in the Hy absorption profile of 8 Lyrae. The effects were tentatively at- 
tributed to a superposed resonance radiation, but no satisfactory interpretation of the 
rather complicated effects could be given at that time. Perhaps the observed effects 
have something to do with a phenomenon of the kind proposed above. 

In 1934 some other rapidly rotating stars were also observed by me, but no clear 
effects of polarization could be found in the absorption lines. It is to be remarked, how- 
ever, that the objects in question have spectral types later than B. 


YNGVE OHMAN 


STOCKHOLM OBSERVATORY 
September 6, 1946. 


5 Nature, 134, 534, 1934; Stockholms Obs. Ann., 12, No. 6, 16, 1936; Pop. Astr. Tidskr., 16, 132, 
1935; see also H. Zanstra, M.N., 101, 250, 1941. 
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REVIEWS 


An Introduction to the Study of Eclipsing Variables. By ZDENEK Kopat. (“Harvard Observa- 
tory Monographs,” No. 6.) Cambridge, 1946. Pp. x+220. $4.00. 


The urgent need for advanced textbooks on astronomical subjects was recently pointed out 
by Struve in his review of a Russian book! with a similar title. The present book, however, is 
not an introduction to the complete field of the study of eclipsing variables but an exhaustive 
treatise of a limited part of the subject. The photometric and geometrical properties of the stars 
are fully treated and an up-to-date survey is given of the methods for analyzing their light- 
curves. However, the wealth of both elementary and advanced information from spectrographic 
observations is completely omitted. 

Throughout the book the emphasis is laid on the mathematical analysis of the subject. As 
a result, the names of actual eclipsing binaries are rarely mentioned except in the Introduction 
and in some footnotes. Numerical examples or tables are not included. Among the 195 entries 
of the Index the name “Algol” is not found. 

The organization of the book is clear. After an introductory chapter, it contains two parts. 
The first part, chapters ii-v, deals with spherically symmetrical stars, while the second part, 
chapters vi-ix, deals with distorted stars. Each part begins with a discussion of the photo- 
metric properties of single stars, proceeds with the computations of the total light for a given 
configuration, and ends with the methods for determining the orbits from the light-curves. 
The abbreviated titles of the chapters are: (i) “Introduction,” (ii) “Nature of the Eclipses,” 
(iii) ‘“‘Determination of the Elements,” (iv) “Improvement of the Elements,” (v) “Orbital 
Eccentricity,” (vi) “Distorted Stars,” (vii) “Ellipticity and Reflection,” (viii) “Light curves 
within Minima,” (ix) ““Elements of Distroted Systems.’”’ The bibliographical notes which con- 
clude each chapter are excellent guides to the original papers. 

The greater the mathematical requirements of a subject, the more careful an author should 
be in its presentation. Each superfluous coefficient or substitution should be omitted, and a 
strong effort should be made to choose clear and consistent notations. This is not only a great 
help to the reader but also a protection against mistakes. In this connection the reviewer wishes 
to comment upon three separate points. 

1. A few minor errors have been noted. In equation (23), cos +/1 — r/r, is written instead of 
/1 — r?/ri; in equations (186), (342), and (333.1), the sign of Vrot should be reversed; the ex- 
pression (88) should be inverted; small errors in some other formulae are obvious. Wrong 
quotations of formulae, which may be traces of earlier drafts, are occasionally found. Further- 
more, the author has not always been fortunate in the choice of his notations. For instance, on 
page 114 the potential YW determining the shape of a star, is written as 


w=U+V+W; 


now the consistent notations of the component terms of W would be Wz and Wrot. The notations 
used on pages 120-22, Vr and Vrot, are confusing. Finally, the book contains only two figures, 
one of which is poor. 

2. The necessity of choosing concise notations is further illustrated by the first part of chap- 
ter viii. In the rectangular co-ordinate system X’Y’Z’ (p. 171), an arbitrary point can be given 
either by its rectangular co-ordinates X’, Y’, and Z’ or by its radius vector r and the direction 
cosines, L = X’/r, M = Y’/r, and N = Z'/r. The direction cosines satisfy the equation (271): 


Now Kopal defines by equation (266) a constant quantity 7, introduced by equation (270) 
the new variables z = Ln, y = Mn, x = Nn, and “adopts x, y, and z as new independent 


1DP. J. Martinoff, Eclipsing Variables. Reviewed in A p. J., 101, 380, 1945. 
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variables which should replace the polar codrdinates r, 6, ¢ used formerly.” This is evidently 
incorrect, since x, y, and z satisfy the identity 


a? y2+4 = 


and cannot be used as co-ordinates for an arbitrary point in space. But the author ignores this 
and infers from this identity that, “to the order of accuracy we are working,” the intersection 
of the primary star with the x~y-plane “‘can still be considered a circle.” 

The full implications of this are not at once obvious. Chapter viii results in an expression for 
the loss of light A¥ of a distorted star eclipsed by another distorted star. Distortion terms up to 
fourth-order spherical harmonics are taken into account. Since the fifth-order terms have the 
same magnitude as the squares of the second-order terms, the expression derived is the most 
general one as far as linear terms are concerned. The total number of terms in A@ is 146 (eqs. 
[282]-{289]). A preliminary examination suggests that the groups of terms denoted by aU, a?, 

Y,f2, fY, and fP are correct, whereas the groups of terms denoted by fY and f? probably need 
a revision. 

Fortunately, the following chapter (ix) is hardly affected. It starts with a description of 
Russell’s classical method of rectifiction of light-curves. The elaborate results of chapter viii 
are used only to demonstrate that Russell’s procedure will give incorrect results as soon as 
third-order harmonics are involved. 

Some remarks may be added to the discussion of the reflection effect in sections 60-63. Kopal 
denotes by /(0, u) the intensity distribution of light emerging under different angles 7 from a 
star’s atmosphere and by /*(0, «) the distribution of “reflected” light. Until Chandrasekhar’s 
recent investigations, these functions could hardly be computed. It is obvious, however, that 
they must satisfy the flux conditions: 


1 1 
f (1) 1*(0, u) cosa, (2) 
0 0 


where u = cos y, 7F = net flux radiated by 1 cm? of the star’s surface, and 7S cos a = incident 
flux on 1 cm? of the star’s surface. It should be noted that equation (1) is not satisfied by Chandra- 
sekhar’s approximations but only by the correct solution to which these approximations converge. 
The first approximation equation (230) gives a net flux which is only the fraction (2 + +/3)/4 = 
0.933 of the correct flux, whereas the fourth approximation is still too small by 0.7 per cent. 

Again the Hopf theorem, equation (225), by which J* can be expressed in terms of J, is valid 
only in the case where the atmospheric particles re-emit or scatter the radiation isotropically. 
By applying the theorem to the approximate flux, equation (230), the error is doubled so that 
the values of 7*(0, u) expressed in equation (231) give an emergent flux which is about (0.933)? = 
0.870 times the incident flux. 

It will therefore appear that, to get a useful approximation to the phase-function of a star for 
integrated light, the function f(e) (p. 153) should be increased by 15 per cent, or the function 
f!V(©) by 1.4 per cent. In most applications, Kopal’s statement that the phase-function nearly 
equals the phase-function G(e)/7, following from Lambert’s law, will be useful. 

Apart from these defects, which affect only certain parts, the book gives a reliable presenta- 
tion of the subject of the theoretical interpretation of light-curves of eclipsing variables. The 
research worker or advanced student will find it a useful presentation of the classical methods 
and a guide to the discussion of more complicated effects. However, it cannot be recommended 
as an introductory volume for students, because its treatment does not stress sufficiently the 
interaction between theory and observations, which is vital for scientific progress. 


H. C. VAN DE HULST 
Yerkes Observatory 


ERRATUM 


On page 289, Volume 104, the caption to Table 5 should read “Measured Equivalent 
Widths,” not “Measured Equipment Widths.” 
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INDEX TO SUBJECTS 
Absolute Magnitudes of Dwarf M Stars, Motions and. A. N. Vyssotsky 


Absorption Coefficient of the Negative Hydrogen Ion, On the Continuous. IIT. S. Chan- 
drasekhar and Frances Herman Breen 


Absorption Lines on the Temperature Distribution of the Solar Atmosphere, The Effect 


Absorption to Reddening, The Ratio of Interstellar. Jesse L. Greenstein 

Absorption Spectrum of Air in the Region AX 600-2000, Ultraviolet. J. J. Hopfield . 
Z Andromedae and AX Persei, The Light-Curves of. Cecilia Payne-Gaposchkin 

R Aquarii, The Light-Curves of. Cecilia Payee Gabpecdlile and Constance Boyd 

6 Aquilae, Spectroscopic Observations of. Carlos U. Cesco and Otto Struve 


a Boétis and 70 Ophiuchi A, A Comparative Study of the Spectra of. Suzanne E. A. 


8 Canis Majoris, Some Radial-Velocity and Line-Intensity Measures in the Spectrum of. 

27 Canis Majoris, The Radial Velocity of. O. Struve 

Clusters, The Diameters of Globular. Albert G. Mowbray 
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